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Statement  of  Research  Probl&m 


Research  conducted  during  the  term  of  this  Grant  was  measurement  of 
infrared  reflectance  or  transmittance  spectra  of  selected  condensed  materials, 
computational  analysis  of  the  measured  spectra  to  obtain  spectral  values 
of  the  complex  refractive  index  N(v)  =  n(v)  +  ik(v)  of  those  materials, 
and  further  analysis  of  the  spectra  in  terms  of  fundamental  intra-  and 
inter-molecular  vibrational  modes. 

II.  Summary  of  Most  Important  Results 

The  most  significant  accomplishments  of  this  project  are  determinations 
of  spectral  values  of  N ( \> )  of  bulk  natural  and  artificial  aerosol/cbscurant 
materials,  and  providing  these  values  of  N(v)  to  scientists  within  the 
aersol/cbscuration  program  at  the  Aberdeen  Proving  Ground,  Md.  We 
measured  and  sent,  to  the  Aberdeen  Proving  Ground  spectral  values  of  N(v) 
for  each  of  the  thirty-eight  solids  and  liquids  listed  in  Table  I,  Most 
of  the  materials  listed  in  Table  I  were  chosen  for  investigation  by 
Mr.  Hugh  Carlon,  CSL,  Aberdeen  Proving  Ground,  who  was  the  technical  monitor 
for  this  research  program. 

Current  knowledge  of  the  manner  in  which  electromagnetic  radiation 
propagates  through  aerosols/obscurants  is  based  primarily  on  Ml  e-scattering 
computations.  Such  computations  are  only  possible  when  there  is  prior 
knowledge  of  M(v)  for  the  aerosol/obscurant.  The  spectral  values  of  N(v) 
obtained  during  the  term  of  this  Grant  are  now  routinely  used  for  Mie- 
scattering  computations  by  scientists  at  many'  U.S.  Army  Laboratories  and 
by  many  DGD  contractors.  Tabulated  values  of  N(v)  for  these  materials 
can  be  obtained  from  the  Aercsol/Cbscuraticr.  Science  Office,  CRTAR-CI3-PS , 
Aberdeen  Proving  Ground,  Md.  21010;  or  by  consulting  publications  associated 


with  this  Grant. 


TABLE  I.  Bulk  materials  investigated  during  the  term  of 
Grant  DAAG29-76-G-0185. 


Material 


Sample 

Concentration 


Spectral  Obtained 


Range 


n  +  ik 


V0* 

5,10,20,40,50,65,75,65%, 
i  Cold  85/5  sample. 

2-33  wn 

Yes 

CL 

20,30,40,50,65,75/5 

i 

2-33  pm 

Yes 

Roe  Oil 

O 

O 

2-33  pm 

Yes 

Diesel  Fuel 

100% 

2-33  pm 

Yes 

DC-704  Silicon 
Oil 

100% 

2-33  pm 

Yes 

Plexiglas 

100% 

2-33  pm 

Yes 

hno3 

1  0.5,1,2,4,8,15.7  M 

2-33  pm 

Yes 

1-hexadecanol 

100% 

2-33  pm 

Yes 

1-cctadecar.ol 

100% 

2-33  pm 

Yes 

1-eicosanol 

100% 

2-33  pm 

Yes 

Futile 

Natural  sample 

2-33  pm 

Yes 

Limestone* 

Natural  sample 

,2-33  um 

Yes 

Water* 

purified  &  deionized 

,42-, 64  pm 

Yes 

Natural  Waters* 

six  samples 

2-20  pm 

Yes 

Graphite* 

crystalline  (EjC)  I 

L.24xl0_6-1240  i 

im  Yes 

Denotes  materials  investigated  for  which  only  partial  support 
was  provided  by  DAAG29-76-G-0185. 
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Table  II  contains  a  listing  of  powdered  materials  pressed  in  KBr 
pellets  for  which  infrared  transmittance  spectra  were  measured  during 
the  latter  months  of  the  Grant.  Determination  of  spectral  values  of 
N(v)  for  the  powdered  materials  was  not  completed  during  the  term  of 
this  Grant,  but  is  continuing  during  a  follow-on  Contract, 


TABLE  II.  Powdered  samples  pressed  in  KBr  pellets  for  which  transmittance 
spectra  were  obtained  during  the  latter  months  of  Grant 
DAAG  29-76-G-0185.  Determination  of  spectral  values  of  n+ik 
has  not  been  completed  for  these  materials. 


Material  Concentration  Spectral  Range 


Montmorillonite 

1,  o.i, 

O 

O 

■0^ 

2— ho 

um 

Cclemanite 

1,  0,1, 

0.01$ 

2-10 

vim 

Remit  e 

1,  0.1, 

0.01$ 

2-10 

pm 

Kaolin 

1,  0.1, 

O 

O 

0^ 

2-35 

pm 

Lampblack 

0.002$ 

2-35 

pm 

Activated  Charcoal 

0.002$ 

2-35 

pm 

Rood  Charcoal 

0.02$ 

2-35 

pm 

Aluminum 

2$ 

2-6 

pm 

Kaolin 

2$ 

2-6 

pm 

Montmorillonite 

^Ot 

C-tO 

2-6 

pm 

Colemarlte 

2$ 

2-6 

pm 

Illite 

2$ 

2-6 

pm 

Remit  e 

2-6 

pm 
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V.  Report  of  Research 

Research  conducted  during  the  term  of  the  Grant,  for  purposes  of  this 
report,  is  divided  into  three  parts:  Investigations  of  Infrared  Reflectance 
which  are  presented  in  Section  Y.A. ,  Investigations  of  Infrared  Transmittance 
which  are  presented  in  Section  V,B. ,  and  Other  Investigations  which  are 
presented  in  Section  V.  U 

A.  Investigations  of  Infrared  Reflectance 

Here,  in  Section  V,A, ,  we  present  the  investigations  of 
those  materials  listed  in  Table  I, 

Experimental  Methods 

A  block  diagram  of  the  reflectcmeter  and  data  acquisi¬ 
tion  system  is  shewn  in  Fig,  1.  Radiant  flux  from  a  glower 
G  was  chopped  at  C  and  was  then  focussed  by  an  f/5  optical 
system,  consisting  of  plane  mirror  and  spherical  mirror 
M^,  on  the  surface  of  the  sample  S,  The  angle  of  incidence 
0  was  6.2  degrees  for  the  central  ray  from  Mg  incident  on 
the  surface  of  S.  Radiant  flux  reflected  from  S  was  imaged 
by  an  identical  optical  system  (M^  and  M^)  on  the  entrance 
slit  of  a  Ferkin-Elmer  Ebert  double-pass  grating  monochromator. 
The  entrance  slits  were  manually  adjusted  to  assure  spectral 
resolution  v/Av  of  100  or  better.  After  passage  through  the 
monochromator,  the  radiant  flux  was  optically  filtered  to 
remove  higher  diffraction  orders,  and  was  then  focussed  on  a 
thermopile  detector  equipped  with  a  Csl  window.  The  signal 
from  the  detector  was  synchronously  processed  by  a  Princeton 
Applied  Research  (FAR)  Rode I  124  lock-in  amplifier,  The 
analog  output  signal  from  the  Model  124  was  applied  to  the 
input  of  a  FAR  Model  260  analog- to-digital  converter  (.ADC) 
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which  both  digitized  the  signal  and  inserted  an  index 

corresponding  to  the  spectral  position  of  the  mono chromat or . 

In  order  to  facilitate  data  acquisition  and  analysis 

the  system  described  above  was  interfaced  to  a  Hewlett- 

Packard  (HP)  Model  9820A  programmable  calculator  equipped 

with  an  HP  9862A  plotter  and  an  HP11223A  magnetic  cassette 

19/ 

tape  unit.  This  entire  system  was  then  interfaced-^- 
through  a  decwriter  1-36  TSO  terminal  to  an  IBM  370/168 
computer  which  was  used  for  the  subsequent  Kramers-Kronig 
analysis  cf  the  data. 

During  this  investigation  we  measured  relative  specular 
reflectance  R(v)  =  R(v) _/R(v)  .  where  s  denotes  the  sample 

O  «V 

and  w  denotes  the  water  standard.  Prior  to  use  the  water 

standard  w  was  purified,  deionized,  and  filtered  through 

a  0.2  mi  Milllpore  filter.  Aqueous  solutions  were  prepared 

using  purified  water  and  reagent  grade  chemicals.  The 

aqueous  solutions  and  purified  water  were  placed  In  separate 

petri  dishes  which  served  as  sample  holders.  Each  dish  'was 

filled  to  precisely  the  same  level  as  determined  by  use  of 

a  cathetometer.  The  levels  of  all  samples  were  monitored 

routinely  during  data  collection  and  purified  water  was 

added  to  the  samples  as  needed  to  compensate  for  evaporation, 

I 

Each  sample  -was  also  stirred  thoroughly  at  the  beginning 
cf  a  spectral  scan  to  assure  isotropic  homogeneity. 

Data  were  acquired  at  200  equally  spaced  wavelength  positions 
In  each  of  the  wavelength  regions  2-3,  3-4,  4-6,  6-8,  8-10, 
f  10-15,  15-20,  20-30,  and  30-40  urn.  For  seme  samples,  individual 

relative  reflectance  measurements  were  made  at  30.15,  32.2, 


( 

I 


34.0,  and  37.1  van,  where  atmospheric  absorption  was  a 
minimum,  rather  than  throughout  the  30-40  van  region. 

Procedures  were  slightly  different  for  handling  the 
solid  samples,  The  solid  samples  we re  placed  one-at-a-time 
on  a  small  three-legged  stand  within  the  reflectometer. 

The  top  surface  of  the  water  standard  and  the  polished 
surface  of  the  solid  sample  were  placed  at  the  same  height 
by  adjusting  the  legs  of  the  stand  while  viewing  the  surface 
edge  through  a  cathetometer,  Additionally,  the  beam  from 
a  He-Ne  laser  was  reflected  from  the  level  surface  of  the 
water  standard  to  a  fiducial  mark  on  the  ceiling  of  the 
laboratory.  The  legs  of  the  stand  were  further  adjusted  to 
bring  the  laser  beam  to  the  same  fiducial  mark  when  the 
water  was  replaced  by  a  solid  sample.  The  latter  leveling 
adjustments  were  made  so  that  the  height  of  the  solid  surface 
was  also  at  the  proper  position.  Reflectance  spectra 
were  then  acquired  for  the  solids  at  the  same  wavelength 
positions  as  for  the  liquids, 

All  spectra  were  obtained  with  the  samples  at  about  27°C, 
except  for  one  cold  sample  of  85$  phosphoric  acid  which  was 
of  about  -10°C.  The  reflectance  spectra  R(v)  were  converted 
from  wavelength  to  wave  number  and  then  plotted  either  by 
the  HP9862A,  a  Tektronix  graphics  terminal,  or  by  a  calcomp 
plotter  interfaced  to  the  IBM  370/163,  The  standard 
deviations,  based  on  three  independent  measurements,  were 
generally  +  0,005  R(\>)  but  increased  to  about  +  0,01  R(v) 
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in  relatively  narrow  spectral  regions  where  atmospheric 
water  vapor  and  carbon  dioxide  were  strongly  absorbing, 

Aa.uecus  Solutions  of  ZnClp 
a.  Objectives 

The  two  principal  objectives  for  this  investigation 
were:  (1)  to  measure  central  frequencies,  half-widths, 
and  band  strengths  of  intra-  and  inter-molecular  infrared 
active  vibrational  modes  of  liquid  water  in  aqueous 
solutions  of  ZnCl9;  and  (2)  to  obtain,  for  aqueous 
solutions  of  ZnC^,  spectral  values  of  the  complex  re¬ 
fractive  index  N(v)  =  n(v)  +  ik(v)  which  are  suitable 

o  o  £5 

for  computing  basic  Mie  scattering  parameters,  These 
two  objectives  were  accomplished  as  follows.  First, 
the  near-normal  incidence  relative  reflectance  spectra 
R(v)  =  R(v;  /R(v)  in  the  270-50C0  cm~'L  spectral  region 

o  W 

were  measured  for  20,3,  29.4,  40,2,  50,  65,  and  75  per 
cent  solutions  by  weight  of  ZnC^  in  water.  In  moles/liter 
the  concentrations  of  the  three  solutions  are  1,77  M, 

2.8  M,  4.19  M,  5.75  M,  8.82  M,  and  11,65  M,  respectively. 
Subscripts  s  and  w  refer  to  the  aqueous  solutions  and 
to  pure  water,  respectively ,  Spectra  R(v)  were  analyzed 
by  use  of  a  subtactive  Kramers-Kronig  algorithm  to  obtain 
spectral  values  for  A<j>(v)  the  differences  in  phase  shifts 
for  electromagnetic  waves  of  frequency  v  reflected  from 
the  surface  of  the  ZnC^  solutions  and  from  the  surface  of 
the  water  which  was  the  reflectance  standard.  Spectral 
values  of  N(v)s  for  the  ZnC^  solutions  were  then  computed 
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by  use  of  an  algorithm^  that  uses  R(v),  A 4> C v ) ,  a 

measured  value  for  the  angle  of  Lncider.ce  e,  and  known 
2/ 

•values—  of  N(v)v/.  The  influence  of  ZnCl-  cn  the  intra- 
and  inter-molecular  vibrational  modes  of  liquid  water 
was  determined  by  critical  examination  of  R(v)  and 
k(v)  throughout  the  spectral  region  investigated, 
b.  Background  Information  for  Discussion  of  ZnCl^  Solutions. 

A  knowledge  of  the  optical  properties  of  aqueous 
solutions  of  ZnCl2  is  of  both  theoretical  and  practical 
importance.  Prom  the  practical  point  of  view,  clouds 
composed  of  droplets  of  such  solutions  are  occasionally 
present  in  localized  geographical  regions  within  the 
earth's  atmosphere.  Analytical  prediction  of  the  scatter¬ 
ing  parameters  for  radiation  passing  through  such  clouds 
relies  primarily  cn  Me  scattering  theory.  This  theory 
requires  prior  knowledge  of  the  spectral  values  of  the 

complex  refractive  index  N(v)  =  n(v)  +  ik(v)  of  the 

So  s 

materials  comprising  the  aerosol. 

The  theoretical  importance  of  aqueous  Zr.Cl~  solutions 
lies  in  a  long-standing  interest  in  the  complicated 
structural  properties  of  liquid  water  and  in  the  forma- 

pi 

tion  of  autoccmplexes  of  the  Zr.  metal  ion.  The  present 
knowledge  of  the  structure  of  liquid  water  and  of  aqueous 
ZnCl2  solutions  is  reviewed  here,  with  emphasis  cn  those 
properties  that  provide  greater  understanding  of  the 
accurate  measurements  of  Lnfrared  reflectance  spectra 
described  hereinafter, 


I 
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Several  intra-  and  inter-molecular  vibrational 

l  bands  are  clues  to  the  structure  of  the  liquid.  For 

instance,  the  infrared  reflectance  spectrum  R(v)  of 

w 

liquid  water,  and  Kramers-Kronig  analysis  of  that  spectrum, 
l-  provide  spectral  values  of  the  phase  shift  $(v)w  of 
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occurs  at  a  higher  frequency  v  than  the  maximum  for 

k(v)  .  Thus  in  the  scientific  literature  there  are 

slight,  tut  implicit,  disagreements  as  to  the  central 

frequency  v  of  sane  vibrational  bands  because  sane 
max 

investigators  work  entirely  with  a(v)w  and  others  -work 
entirely  with  k(v)  .  Noting  this  implicit  difference 
in  obtaining  central  frequencies  vnHX  of  vibrational 
bands,  we  relate  to  k(v)w  the  subsequent  discussions  of 
central  band  frequencies  half  widths  r,  and  band 

strengths  S^. 

In  the  vapor  state  the  H20  monomer  possesses  C_v 
symetry  and  has  three  fundamental  vibrational  modes : 
v^(A^)  symmetrical  0-H  stretch  at  3656  cm~J", 
v2(A^)  OH— 0  bending  motion  at  1591*  cm,  and 
v^CB^)  antisymnetric  OH  stretch  at  375 6  cm-1. 

It  is  very  doubtful  that  the  H20  monomer,  per  se,  exists 
in  the  liquid  state.  Hydrogen-bonded  temporal  clusters 
of  several  H20  molecules,  however,  are  believed  to  exist 
in  liquid  water.  An  infrared  spectrum  of  liquid  water 
therefore  is  interpreted  in  terms  of  bands  associated 
with  the  intra-mo lecular  vibrations  of  H20  molecules, 
and  in  terms  of  inter-molecular  lattice  vibrations  of 
the  clusters.  For  example,  the  infrared  spectrum  of 
liquid  water  at  27°C  possesses  a  band,  designated  as 
\>2,  that  is  associated  with  the  O-K-O  bending  motion. 

The  v~  band  has  a  central  frequency  v  of  16^0  cm-1, 
a  half  width  r  of  110  cm~~,  and  a  band  strength  of 
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PP  P 

(^.5  +0,8)  x  10”  cm  /molecule.  The  symmetric  and 
antisymmetric  G-K  stretching  vibrations,  designated 
and  v,,  respectively,  and  the  first  overtone  2v2 
which  may  be  in  Fermi  resonance  with  one  of  the  C-H 
stretching  vibrations,  provide  broad  overlapping 
bands  in  relatively  close  spectral  proximity  to  one 
another.  The  bands  associated  ’with  v1 ,  v,,  and  2v2 
are  grouped  together  for  purposes  of  analyzing  the 
infrared  spectrum;  together  they  provide  a  very  broad 
band  with  central  frequency  v  of  3395  cm”*,  a 

UlCwC 

half  width  r  of  390  cm”\  and  a  combined  band  strength 

of  (3.65  +  0, 08)  x  10“21  cm2/molecule. 

The  structure  of  the  temporal  clusters  of  H20 

molecules  in  liquid  water  is  not  fully  understood  at 

V 

this  time.  However,  one  recent  model—  is  useful 
when  interpreting  the  infrared  spectrum  of  liquid 
water.  The  model,  regarded  by  its  authors  as  net 
being  definitive  until  further  experimental  verification 
is  obtained,  consists  of  an  r^O  molecule  hydrogen- 
bonded  with  C2v  symmetry  in  a  temporal  cluster  of  four 
first  and  twelve  second  nearest  neighboring  H20 
molecules.  A  normal  coordinate  and  Monte  Carlo  analysis 
of  this  structure  yielded  six  principal  lattice 
vibrational  modes  with  a  total  of  thirty  different 
vibrational  species.  According  to  the  normal  coordinate 
analysis  the  six  principal  lattice  vibrations  occur  at 
ca,  75,  165,  219,  ^50,  550,  and  722  cm-*.  In  the 
infrared  spectrum  of  liquid  ’.rater,  assuming  that  the 

•  T  - — * 
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seventeen  molecule  cluster  with  C2,r  symnetry  has 

significant  validity,  the  lattice  vibrational  modes 

near  450,  550,  and  722  cm  ^  provide  broad  overlapping 

bands  that  collectively  give  a  broader  band  designated 

as  the  librational  band  vT ,  From  infrared  spectra 

obtained  for  liquid  water  at  27°C  the  parameters  for  the 

1  —1 

vT  band  are  v  =  580  cm-x,  r  =  500  cnf~,  and  S.  = 

Li  max  U 

(7.2  +1.5)  X  10~Ci  cm2/  molecule.  Similar  parameters 
for  the  other  three  principal  lattice  vibrations  have 
not  been  measured  because  infrared  reflectance  and 
transmittance  spectra  of  water  are  difficult  to  obtain 
in  the  region  below  300  cm-2. 

Of  particular  relevance  to  the  investigation  of 
ZnCl2  solutions  was  the  fact  that  the  nature  of  water 
molecules  in  the  liquid  state  could  also  be  investigated 
through  effects  that  monatomic  and  polyatomic  ions 
have  on  the  intra-  and  inter-molecular  vibrations  of  the 
temporal  clusters.  The  ions  are  sources  of  electric 
fields  that  interact  with  the  electric  dipole  moments 
of  neighboring  water  molecules  within  the  hydrogen- 
bonded  temporal  clusters.  In  infrared  reflectance 
spectra  R(v)  of  aqueous  solutions,  and  in  associated 
spectra  N(v)s,  the  averaged  effects  of  such  interactions 
are  observed  as  (1)  shifts  in  the  central  frequencies 

of  the  various  intra-  and  inter-molecular  vibrational 
max 

bands,  and  (2)  as  changes  in  the  infrared  band  strength 
associated  with  each  vibrational  band.  Investigations  of 
these  effects  are  described  in  several  previous  articles,— 


3/ 

D.  E.  Irish-  previously  reviewed  the  physical 

and  cherrical  properties  of  fused  salts  and  electrolyte 

solutions  of  ZnClo.  Hainan  spectra  were  the  primary 

basis  for  deducing  the  chemical  properties,  but  several 

references  to  infrared  studies  were  also  included. 

3ased  on  the  review  by  Irish,  and  the  pertinent  publica- 
9-1 3/ 

tions— '—  referred  to  therein,  we  note  that  aqueous 
solutions  of  ZnC^  contain  the  polyatomic  constituients 
ZnCl+,  ZnCl2>  ZnCl,-,  ZnCl^-,  and  Zn(H20)g1+  in  addition 
to  hydrated  species  of  seme  of  these  autoccmplexes. 

The  diatomic  ion  ZnCl+  has  a  vibrational  frequency  of 
ca.  305  cm-1  and  occurs  in  solutions  of  less  than  4M 
concentration.  The  linear  molecule  ZnC^,  which  ccc’urs 
in  solutions  of  greater  than  4M  concentration,  has  three 
fundamental  vibrational  modes:  Cl)  the  v,  Zn-Cl 
symmetric  stretch  at  ca.  305  cm-'1,  (2)  the  v0  Cl-Zn-Cl 
bending  mode  at  ca,  101  cm-1,  and  (3)  the  v,  Zn-Cl  anti¬ 
symmetric  stretch  at  ca.  505  cm-  ,  The  v..  mode  is 
Raman  active  and  the  V2  and  v,  modes  are  infrared  active. 
Mo  evidence  from  Raman  or  infrared  studies  exists  to 

either  prove  or  disprove  the  existence  of  ZnCl^~  in 

2- 

aqueous  solutions  01  ZnC^.  The  anion  ZnClij  ,  which 
exists  in  solutions  of  less  than  10  M  concentration, 
possesses  the  nine  fundamental  vibrational  modes 
characteristic  of  tetrahedral  molecules:  v. (A-,  )  at 
278  cm-',  v-,(Z)  at  ca.  100  cm-1,  at  ca.  280  cm-', 

and  vjj(F;)  at  ca.  122  cm-'.  All  modes  are  Pamar.  active, 
but  only  the  v,(F-,)  and  V]iC?2>  triply  degenerate  modes 
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2+ 

are  infrared  active,  The  octahedral  cation  Zn(H20)g  , 

which  exists  in  solutions  of  less  than  13  M  concentration, 

possesses  three  Raman  active  vibrational  species  A.,  , 

E  ,  and  T0  and  two  triply  degenerate  T-,,,  infrared  active 

vibrational  species.  Irish  observed  a  Raman  band  at 

ca.  390  cn-~  which  was  assigned  to  the  A,  species. 

Infrared  active  modes  attributed  to  Zn CH^O ) g  were 

observed  in  ZnSO^  •  71^0  by  use  of  the  KBr  disk  method  by 

9/ 

Nakagawa  and  Shimanouchi ,—  They  assigned  an  0H2 

wagging  vibration  to  a  band  at  541  cm-^,  a  Zn-0  stretch¬ 
ing  vibration  to  a  band  at  364  cm-1;  and  an  OIL,  rocking 
vibration  that  was  calculated  to  occur  at  620  cm-"  was 
not  observed  possibly  because  of  overlap  with  the  band 
due  to  the  vjj(F2)  mode  of  SO^-.  In  aqueous  solutions 
of  greater  than  10  M  concentration  a  polynuclear 
aggregate  of  ZnCl^  tetrahedra,  comparable  to  that  of 
crystalline  ZnCU,  was  observed  by  Irish  to  possess 
Raman  active  vibrational  modes  at  ca.  90,  230,  335, 
and  280  cm-'''.  Evidence  of  the  water  structure  also 
persists  in  Raman  spectra  for  ZnCl2  solutions  up  to 
about  10  M  concentration, 
c.  R(v)  and  N(v)  Spectra  for  ZnCl2  Solutions 

The  measured  relative  reflectance  spectra  R(v) 
for  the  1.77,  2,8,  and  4, 19  M  ZnCl9  solutions  are  pre¬ 
sented  in  Pigs.  2  and  3,  and  R(v)  for  the  5.75,  3.82, 
and  11.65  M  ZnCl2  solutions  are  presented  in  Figs,  8 


WAVE  NUMBER  (cm*1) 


X  10° 


Figure  Index  of  refraction  n(v)  in  the  spectral  region 
from  270  to  2500  cm-1  for  20.3,  29.4,  and  4o,2 
percent  aqueous  solutions  containing  ZnCU. 
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dashed  line,  and  75?  lonp-dashed  line. 


dashed  line,  and  75$  long-dashed  line. 


65%  short-dashed  line,  and  75%  long-dashed  line. 
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and  9.  The  P.(v)  spectra  were  phase-shift  analyzed  by 
use  of  the  Kramers-Kronig  relation 

co 

a<^(xc)  =  Prin.(x0A)  f  Un[R(x)]/(xAx2)}dx, 

(3) 

where  A(fc(X  1  is  the  difference  in  chase  shift  for 

U 

electromagnetic  waves  of  wavelength  XQ  reflected  from 

the  surfaces  of  the  sample  and  the  water.  The  portions 

of  the  integral  in  Eq.  (3),  beyond  the  region  in  which 

data  were  acquired,  were  evaluated  analytically  by 

assuming  constant  relative  reflectance  R(X)  =  R(X„j^) 

in  the  region  0<X<X  .  and  R(X)  =  R(X  )  in  the  region 
min  rriaX 

x  <x«»1 

Spectral  values  of  N(X)  were  then  obtained  for  each 
solution  by  use  of  the  results  from  the  Kramers-Kronig 
analysis  as  input  to  a  previously  published  algorithm. 
Values  of  N(X)  are  presented  in  Figs,  4-7  for  the  10, 

20,  and  30!?  ZnCl2  solutions,  and  in  Figs,  10-13  for 
the  50,  65,  and  75%  ZnCl2  solutions, 

Discussion  of  Spectra  of  ZnCl?  Solutions 

2+ 

Significant  hydrolysis  of  Zn  dees  not  occur  in 
10  19/ 

aqueous  solutions — 1 -  of  ZnCl2.  Thus  the  infrared 

o  «  ^ 

bands  of  ZnC^COg  at  620,  541,  and  364  cm--1-  provide 
only  weak  contributions  to  the  spectra  shown  in  Figs, 
2-13. 

1  0/ 

Raman  intensity  studies  by  Irish  et  al ,—  indicate 
that  below  4  M  concentration  that  ZhCl  and  ZnCl^  are 
predominate,  and  that  concentrations  of  ZnCl-  are 
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insignificant,  ZnCl+  is  present  in  concentrations  of 

about  1.3,  0.4,  and  0,  M  in  the  1.77,  2.8,  and  4.19  M 

2_ 

ZnCl2  solutions,  respectively,  ZnCl^  is  present  in 
concentrations  of  about  0.3,  0,8,  and  1,6  M  in  the 
1.77,  2.8,  and  4, 19  M  ZnCl,  solutions,  respectively. 

In  the  4,19  M  aqueous  Zr.Clo  solution  the  ZnCl0 
molecule,  per  se,  is  present  at  about  0.1  M  concentra¬ 
tion,  Ln  the  infrared  spectra  obtained  for  the  1.77, 
2.8,  and  4,19  M  solutions  in  the  region  from  270  to 
5000  cm-x  we  expect,  in  addition  to  bands  due  to  the 
water  structure,  broad  infrared  active  bands  at  305  cm”* 
due  to  ZnCl+,  at  280  cm-^  due  to  the  mode  of 

ZnCl^  ,  and  in  the  4.19  M  solution  weak  contributions 
at  505  cra-^  due  to  the  v,  antisymmetric  stretch  of 
ZnCl2.  Additionally,  ions  Zn^  and  Cl-  will  effect  the 
positions  and  strengths  of  infrared  bands  due  to  the 
water  structure. 

The  Raman  intensity  studies  by  Irish  et  al.— ' / 
also  indicate  that,  for  concentrations  above  4  M,  the 

px  p_ . 

species  Zn  ,  ZnCl2,  and  ZnCljj  are  predominate.  ZnCl2 
is  present  in  respective  concentrations  of  about  0.4  and 
1.9  M  in  the  5,75  and  8,82  M  ZnCl9  solutions.  The 
concentration  of  ZnClp  for  the  11, 65  M  solution  is  not 
known  but  probably  exceeds  2M,  ZnCl^  is  present  in 
concentrations  of  about  2.2,  2.7,  and  3,7  M  in  the 
5.75,  3.82,  and  11,65  M  solutions,  respectively.  Thus, 
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in  the  infrared  spectra  for  the  5.75,  3.82,  and  11.65  M 
solutions  we  expect,  in  addition  to  bands  due  to  water, 

broad  infrared  bands  at  280  cm-"  due  to  the  v^CFp) 

2—  — ' 

mode  of  ZnCl^  and  at  505  cm  due  to  the  anti- 

2-1- 

symmetric  stretch  of  ZnClp.  Additionally,  the  Zn 
ions  will  effect  the  positions  and  strengths  of  the 
infrared  bands  due  to  water. 

The  relative  reflectance  spectra  shown  in  Figs, 

2,  3,  8,  and  9  possess  three  general  spectral  features 
of  the  same  type.  For  example,  the  feature  extending 

T  "1 

from  270  to  ca.  1350  cm-  and  centered  at  ca.  800  cm-*. 
This  feature  is  characterized  by  a  decrease  in  the 
magnitude  of  the  relative  reflectance  in  the  region  from 
270  to  about  720  cm-*,  an  increase  in  relative  reflectance 
in  the  region  from  about  720  to  about  900  cm-J_,  and 
then  a  decrease  in  relative  reflectance  in  the  region 
from  about  900  to  about  1350  cm-'1.  The  two  other 
similarly  shaped  spectral  features  are  centered  at  ca. 

1680  and  3500  cm-\  These  features  are  due  primarily 
to  different  water  content  of  the  aqueous  solutions. 

Note  also  that  the  relative  reflectance  of  liquid 
water  would  be  1,0  throughout  the  spectrum, 

I  i 

'  The  three  spectral  features  are  characterized  in 

the  spectra  n(v)  vs,  v  shown  in  Figs.  4,  5,  10,  and  11 

I 

as  regions  of  anomalous  dispersion,  and  as  peaks  in  the 
spectra  k(v)s  vs.  v  shewn  in  Figs,  6,  7,  12,  and  13. 

1 
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In  applying  the  Kramers-Kronig  analysis  to  the 
relative  reflectance  spectra  R(v)  we  assumed  that  in  the 
wave-number  region  from  0  to  270  cm-"  that  R(v)  = 

R(270  cm  x) .  This  is  a  somewhat  poor  assumption  because 
there  are  other  infrared  active  vibrational  modes  of 


polyatomic  constituents  of  Zn  and  Cl  beyond  270  cm-". 
Additionally,  the  full  influence  of  the  ZnCl+  band  at 
305  cm  -  and  the  ZnCl^  band  at  280  cm  extend  to  wave 
numbers  below  270  cm  ",  the  beginning  of  the  relative 
reflectance  spectra  collected  in  these  experiments. 

Also,  the  505  cm  "  band  of  ZnCU  occurs  in  the  spectral 
region  of  the  strong  \u  band  of  water.  Therefore,  it 
is  virtually  impossible  to  extract  any  information  about 
ZnCl  ,  ZnCU,  and  ZnCl^~  from  these  spectra  except  that 
the  rapid  increase  in  relative  reflectance  in  the  region 
from  600  to  270  cm-"  is  undoubtably  due  in  part  to 
infrared  activity  of  these  three  polyatomic  constituents. 


Reflectance  spectra  in  the  region  below  270  cm-", 
obtained  by  Fourier  Transform  techniques  and  then 
joined  with  the  spectra  reported  in  this  paper,  can 

provide  additional  knowledge  of  the  infrared  parameters 

+  2— 

of  ZnCl  and  ZnCl^  In  these  solutions.  The  increase  in 

■I  -7 

relative  reflectance  from  about  700  cm-1  to  2'70  cm-"  is 
also  due  to  the  vT  band  of  liauid  water  shifting  to  lower 


wave  numbers  as  the  ZnG^  concentration  increases. 

The  assumption  that  R(v)  =  R(270  cm-")  throughout 

*1 

the  spectral  region  from  0  to  270  cm-"  will  have  small 
influence  on  values  computed  for  n(v)  and  k(v)  for 

S  o 


"'■‘Ca.-v.' 


TABLE  IV.  Parameters  for  intra-  and  inter-molecular  vibrations  of  water  in  aqueous  solutions 
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greater  than  about  50C  cm-*. 

The  spectra  k(v)  vs,  v  shewn  in  Figs,  6,  7,  12, 
and  13  provide  parameters  kCv)^,  v^,  r,  and  5^ 
listed  in  Table  IV  for  the  infrared  bands  vT ,  vOJ 

u  c. 

and  the  combined  contributions  of  v,,  and  2x>2 

of  liquid  water.  Referring  to  Table  IV  we  note  for  the 

infrared  band  comprising  Intramolecular  vibrations 

v-, ,  v?,  and  2v^  that  as  the  ZnC^  concentration  increases, 

vmax  3^-^'ts  t0  hi?ber  frequencies,  r  increases, 

decreases,  and  3^  decreases  and  then  either  increases 

or  remains  about  the  same.  Similarly,  as  the  ZnCT^ 

concentration  increases,  we  note  for  the  infrared  band 

comprising  intramolecular  vibration  v0  that  v 

2  max 

shifts  to  lower  frequencies,  r  decreases,  ’K(v)-iax 
increases,  and  3^  increases.  We  also  note  from  Fig.  6 
that  the  contribution  of  the  band  to  k(v)  is  super¬ 
imposed  with  a  general  absorption  which  decreases  slowly 
from  1000  to  2 H00  cm~±.  This  general  absorption  is  due 

to  teirporal  water  clusters  and  is  possibly  due  in  part 

2+ 

to  water  molecules  coordinated  with  Zn  ions.  The 
general  ^sorption  also  contributes  to  the  values  tabulated 
for  k(v)may  fbr  the  v 2  band,  For  the  vT  band  com¬ 
prising  inter-molecular  ’/ibraticns,  as  the  ZnCl?  con¬ 
centration  increases,  v  „  shifts  to  lower  frequencies, 
r  remains  about  the  same,  kCv)^^  decreases,  and 
increases.  Two  values  of  S^  are  tabulated  for  the  vT 
band  in  each  solution;  the  lower  value  in  Table  IV  denotes 
deletion  of  tie  general  absorption  which  was  manually 


wir 


I 

I 

I 

I 


extended  into  the  spectral  region  spanned  by  the 


v-  band,  and  the  higher  value  includes  the  contribution 
from  this  general  absorption.  The  integration  indicated 
in  Eq.  (1)  was  made  manually  in  each  case  by  use  of  a 
planimeter.  The  limits  on  the  integrals  were  from  2600 


to  4000  am  x  for  the  v-, ,  v^,  and  2v2  band,  from  1280  to 

1810  cm~~  for  the  v0  band,  and  from  v  to  1140  cm-" 

<;  max 


for  the  band. 


Parameters  r  and  3^  for  the  v.  band 


were  obtained  by  considering  only  the  high-frequency 


side  of  the  band  and  then  multiplying  the  results  by  2 


before  entry  in  Table  IV.  The  increase  in  measured 
for  the  Vy  band  of  liquid  water  shown  in  Table  IV  is 

XJ 

probably  due  to  infrared  activity  of  ZnCH^O)^  at  620 
and  541  cm-"  in  all  three  solutions,  and  also  to 
infrared  activity  of  ZnCl2  at  505  cm-1  in  the  4,19  M 
solution. 


Parameters  r  and  were  not  measured  for  the  vT 
water  band  in  the  5*75,  8,82,  and  11,7  M  solutions, 

because  of  significant  contributions  of  unknown  magnitude 

2- 

due  to  bands  associated  with  ZnCl9  and  ZnCl^  . 


*3 
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I 

2.  Aqueous  Solutions  of  HgPOii 

|  a.  Objectives 

I 

The  principal  objectives  for  this  investigation 
were:  (1)  to  obtain,  for  aqueous  solutions  of  H-PO^, 

!  spectral  values  of  N(v)  =  n(v)  +  ik(v)  which  were  suitable 

L 

for  computing  Mie-scattering  parameters,  and  (2)  to  measure 
central  frequencies,  half-widths,  and  band  strengths 
of  inter-  and  intra-molecular  infrared-active  vibrational 
modes  of  liquid  water,  H3PO4,  and  molecular  products  of 
due  to  dissociation  of  H3PO11  in  aqueous  solution.  The 
methods  used  to  accomplish  these  objectives  were  the 
same  as  those  described  in  Section  V,A,2,a,  Nine 
phosphoric  acid  samples  were  investigated:  5,  10,  20,  40, 
50 ,  65,  75  solutions,  an  85%  solution  at  25°C,  and  an 
85%  solution  at  about  -1Q°C.  In  terms  of  molar  concentra¬ 
tions  the  H3PO4  solutions  were  0,52,  1.07,  2.27,  5.12, 

I  6.81,  9.78,  12,08,  and  14.65M,  respectively, 

b.  Background  Information  for  Discussion  of  EoPO^  Solutions 
Knowledge  of  the  optical  properties  of  aqueous 

► 

solutions  of  K3PO4  is  of  both  practical  and  basic 
importance.  Prom  the  practical  point  of  view,  clouds, 
which  are  ultimately  composed  of  H3PO4,  are  produced 
in  the  earth's  atmosphere  by  combustion  of  white,  yellow, 
or  red  phosphorus.  Computation  of  radiation  transport 
through  such  clouds  relies  primarily  on  Mie-scattering 
theory,  which  requires  prior  knowledge  of  N(v)  of  the 
cloud's  constituents. 


- - - - IV.- 


I 
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The  basic  importance  of  aqueous  H^PO^  solutions 
lies  in  an  Interest  in  the  hydrogen  bonded  structure  of 
liquid  water  and  the  stoichiometry  of  acid  solutions. 
Present;  knowledge  of  the  structure  of  liquid  water,  as 
it  applies  to  this  investigation,  was  reviewed  in 
Section  V.A.2.b.  The  known  structure  and  stoichiometry 
of  aqueous  solutions  of  are  reviewed  here,  with 

emphasis  on  those  properties  providing  greater  under¬ 
standing  of  the  accurate  measurements  of  infrared 

reflectance  spectra  described  hereinafter, 

20/ 

Preston  and  Adams —  investigated  Raman  spectra 
of  aqueous  solutions  of  H^P04  ranging  in  concentration 
from  0.005  M  to  15.6  M,  They  found  that  the  Raman  spectra 
could  be  reasonably  accounted  for  by  liquid  water,  un- 
dissoclated  H3PO4,  and  the  anion  K2P04”,  The  degree  of 
dissociation,  according  to  their  investigations, 
falls  from  about  0.5  at  0,01  M  concentration  to  0,1  in 
the  1-2  M  range  of  concentrations,  and  remains  at  about 
0.1  up  to  concentrations  of  15.7  M,  Based  on  a  molecular 
model  of  H3PO4  of  the  form  OPX^,  where  X  designates  0-H, 
possessing  C^v  symmetry,  Preston  and  Adams  assigned  the 


following 

bands : 
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These  bands  would  also  be  infrared  active. 

The  OPX,  molecular  model,  however,  dees  not  provide 
for  G-H  vibrations.  The  reduced  representation  of 
molecular  H^PO^  with  C~  point  group  symmetry  is 
r  =  5A,  +  A„  +  6E, 

1  C. 

The  At  species  are:  O-H  stretch,  P-0  stretch,  P-0, 
stretch,  PO,  deformation,  P-O-H  inplane  bend,  and  P-O-H 
out-of-plane  bend.  The  A2  species  is  a  P-O-H  out-of¬ 
plane  bend.  And,  the  doubly  degenerate  E  species  are: 

C-H  stretch,  ?-0,  stretch,  PO,  deformation,  P-O-H  inplane 
bend,  and  P-O-H  out-of-plane  bend.  The  A^  and  E  species 
are  infrared  and  Raman  active,  and  the  A2  species  is 
inactive.  In  the  absence  of  previous  normal  coordinate 
analysis  of  H^PO^,  we  are  currently  undertaking  a  normal 
coordinate  analysis  to  better  assign  the  bands  of  H^PO^. 

The  structure  of  the  anion  HjPO^-  is  not  clearly  under¬ 
stood  in  aqueous  solutions,  and  the  possibility  of 
hydrogen  bonding  makes  it  doubtful,  that  the  anion,  per  se, 
exists  in  the  solutions.  However,  as  with  the  HNO, 
solutions  which  are  discussed  in  Section  V.A.4,,  the 
infrared  spectra  of  H,PO^  will  also  be  characterized 
by  bands  due  to  H0FCV~,  The  anion  H0PC^~  carnet  possess 
symmetry  greater  than  C2v,  for  which  the  reduced 
representation  is 

r  =  6At  +  2A2  +  4B-,  +  3B2  • 

The  At,  3. ,  and  B,  species  are  infrared  active.  If 

X  i.  c. 

the  symmetry  is  C2,  rather  than  C-v,  then  the  reduced 
representation  is 


'■’T'+ff  "*-***'  '  — 
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r  =  3A  +  73; 

and  the  A  and  3  species  are  infrared  active, 

21/ 

Chapman  and  Thirlwell—  previously  investigated 
the  infrared  and  Raman  spectra  of  orthophosphates  and 
assigned  the  following  bands  for  H-PO^- : 


2900-3000  cm  ~ 

0-H  stretch 

238C 

-1 

cm 

O-H  stretch 

1830 

cm 

Combination  band 

1230 

cm”J' 

P-O-H  in-plane  deformation 

1076 

cm”1 

PO-,  stretch 

988 

cm~^ 

P03  stretch 

862 

cm"  ^ 

P-O-H  stretch 

520-  590 

cm 

0P0  bending 

455 

cm-^ 

0P0  bending. 

c.  R(v)  and  N(v)  Spectra  of  H^PO^  Solutions, 

The  measured  relative  reflectance  spectra  R(v)  in 
the  spectral  region  from  0-2,500  cm-1  for  5,  10,  20,  and 
40  per  cent  solutions  of  H^PO^  are  presented  in  Fig,  14, 

The  R(v)  spectra  for  these  four  relatively  weak  aqueous 
solutions  varied  little  from  a  value  of  R(v)  =  1,0  through¬ 
out  the  2,500-5,000  cm-"''  region,  and  therefore  are  not 
presented  graphically  in  this  report.  The  R(v)  spectra 
for  50,  65,  75,  and  85  per  cent  aqueous  solutions  of 
H.,?^  are  presented  in  Figs.  19  and  20,  and  the  R(v) 
spectra  for  the  85  per  cent  solution  at  25°C  and  -10cC 
are  presented  in  Figs,  25  and  26. 

The  relative  reflectance  spectrum  of  each  solution 
was  analyzed  by  use  of  the  Kramers -Krcnig  techniques 
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Figure  14.  ^ear-normal  incidence  relative  specular  reflectance 
R(v)  in  the  300-2,500  cm  ^  spectral  range  for 
5,  10,  20,  and  4 o  percent  aqueous  solutions  of 

h3po4. 
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5.  Index  of  refraction  n(v)  in  the  300-2,500  cirT" 
spectral  range  for  5,  10,  20,  and  ko  percent 
aqueous  solutions  of  H^PO^. 
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Figure  16.  Index  of  refraction  n(v)  in  the  2,500-5,000  cm-1 
spectral  range  for  5,  10,  20,  and  40  percent 
aqueous  solutions  of 


0  500  1000  1500  2000  2500 

WAVE  NUMBER  (cm-') 

17.  Extinction  coefficient  k(v)  in  the  300-2,500 
cm”*  spectral  range  for  5,  10,  20,  and 
percent  aqueous  solutions  of  H,PCi,, 
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'igure  20.  Near-nonral  Incidence  relative  specular  reflectance 
R(v)  in  the  2, 500-5 ,000  cm”^  spectral  range  for 
50,  65,  75,  and  85  percent  aqueous  solutions  of 

f^POjj, 
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Figure  22.  Index  of  refraction  n(v)  in  the  2, 500-5, COG  cm' 
spectral  range  for  50,  65,  75,  and  85  percent 
aqueous  solutions  of  H^PO^, 
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gure  23 


.  Extinction  coefficient  k(v)  in  the  300-2,500  cm' 
spectral  range  for  50,  65,  75,  and  85  percent 
aqueous  solutions  of  H^PC^, 
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Figure  2U.  Extinction  coefficient  k(v)  in  the  2,500-5,000  cn~~ 
spectral  range  for  50,  65,  75,  and  85  percent 
aqueous  solutions  of  H^PC^ , 


spectral  range  for  »5  pe 
and  -10°C  (dashed-line). 


(dashed  line). 
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described  in  Section  V.A.2.C,  to  obtain  'J(v)  =  n(v)  +  ik(v) 
for  each  solution.  Spectral  values  of  N(v)  for  the  5, 

10,  20,  and  40  per  cent  solutions  of  H,PO^  are  pre¬ 
sented  graphically  in  Figs.  15-18,  for  the  50,  65,  ”5, 
and  85  per  cent  solutions  of  H^POi^  in  Figs,  21-24,  and 
for  the  85  per  cent  solution  of  H^POjj  at  25°0  and  -1C°C 
in  Figs.  27-30. 

d.  Discussion  of  Spectra  of  H,P0^  Solutions 

The  degree  of  dissociation  is  greater  for  H-^PO^ 
solutions  of  less  than  about  1  M  concentration.  Thus, 
infrared  active  bands  due  to  H^PO^-  should  be  relatively 
more  prominent  in  the  5  and  10$  solutions.  One  of  the 
strongest  bands  for  I^PO^-  occurs  at  10?4  cm”*,  A 
small  shoulder  band  at  this  position  was  observed  in  the 
5,  10,  and  20$  solutions;  for  example,  examine  the  k(v) 
spectra  in  Fig.  17.  The  1074  cm-1  band  is  not  at  all 
apparent  in  the  solutions  of  greater  concentration,  even- 
though  in  the  85$  solution  approximately  8,5$  of  the 
solution  should  be  HgPO^"',  Therefore,  we  interpreted 
the  spectra  of  the  H^PO^  solutions  on  the  basis  of  infrared 
active  bands  of  H^PO^  and  water.  The  central  frequencies 
v(cm_1)  and  tentative  assignments  for  the  infrared  active 
bands  of  H^PO^  solutions  are  presented  in  Table  V. 
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Table  V,  Central  frequencies  v(cm”“)  of  infrared  active  vibrational  bands 
and  tentative  band  assignments  for  H^POh  solutions. 


1  -  - - 

1 

Molecule 

— 

Assignment 

v(cnT-)  for  infrared  bands  of 

■  ■ 

^PQii  Solutions 

10% 
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VL 
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541 

4.  to  ueous  Solutions  of  HNC3 
a.  Objectives 

The  principal  objectives  for  this  investigation  were: 

(1)  to  obtain,  for  aqueous  solutions  of  HNO,,  spectral  values 
of  N(v)  =  n(v)  +  ik(v )  which  were  suitable  for  computing 
Mie-scattering  parameters,  and  (2)  to  measure  central  frequen¬ 
cies,  half-widths,  and  band  strengths  of  inter-  and  intra¬ 
molecular  infrared-active  vibrational  modes  of  liquid  water, 
HNO,,  and  molecular  products  due  to  dissociation  of  HNO-;  in 
aqueous  solutions.  The  methods  whereby  these  objectives 
were  accomplished  where  the  same  as  those  presented  in  Section 
V,A.2.a,  Six  nitric  acid  samples  were  investigated:  0,5,  1, 


*»■«  ■**'<*• 'Xu*,?'  -1t--  , 
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2,  4,  8,  and  15.7  M  solutions  at  27°C. 

Spectral  values  on  N(v)  for  aqueous  solutions  of 
HNO,  are  important  because  this  acid  is  a  common  aerosol 
in  the  Earth’s  atmosphere. 

The  investigations  as  presented  in  this  section,  were 
recently  accepted  for  publication  in  The  Journal  of 
Chemical  Physics. 

b.  Background  Information  for  Discussion  of  HNO^  Solutions. 

Here  we  present  measurements  of  relative  infrared 
reflectance  spectra,  computed  values  of  spectral  complex 
refractive  index  N(v)  =  n(v)  +  ik(v),  central  'wave- 
number  positions  of  infrared  vibrational  bands,  and  measured 
strengths 

^  ■  g  'band  KW  *>  (4) 

of  selected  bands  for  0,5,  1,  2,  4,  8,  and  15.7  M 
aqueous  solutions  of  nitric  acid.  In  Eq,  (4)  S^  is  the 
band  strength,  pQ  is  the  molecular  density  (mol./cm^),  and 
the  integral  is  evaluated  over  the  full  'width  of  the 
band.  There  have  been  many  previous  investigations  of 
the  infrared  and  Raman  spectra  of  nitric  acid.  None  of 
the  previous  investigations,  however ,  were  of  the  reflectance 
spectra,  and  none  presented  spectral  values  of  N(v)  for 
liquid  nitric  acid. 

Investigations,  prior  to  I960,  of  the  infrared  and 

Raman  spectra  of  nitric  acid  were  previously  reviewed  and 

2~>  / 

summarized  by  Stem  et  al.~  Seme  structural  features 
of  the  HNCU  molecule  were  not  well  known  at  that  time 


Figure  31. 


The  structure  of  nitric  acid:  r^  =  0.121  ran, 
V2  =  0.112  nm,  r-,  =  0,1405  nm,  r^  =  0,0961  nm. 
a  =  130°  13',  =  115°  55',  62  =  H3°  52’, 

Y  =  102°  13’  (frcm  Fef.  (23)). 
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gure  32.  Degree  of  dissociation  a  vs,  molar  concentrat 
M  of  aqueous  solutions  of  nitric  acid.  The 
dissociation  of  HNO^  in  water  follows  Eq.  (5) 
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and  the  molecule  was  therefore  treated  by  sene  investiga¬ 
tors  as  an  O'NC^  planar  system  belonging  to  the  point 
group.  Thus  there  were  several  discrepancies  in  assign¬ 
ment  of  the  fundamental  vibrational  bands  of  HNCU. 

2*2/ 

In  1965  McGraw  et  al.-^  reported  a  normal  coordinate 
analysis  of  a  planar  nitric  acid  molecule  belonging  to 
the  C  point  group  and  possessing  the  structure  shown 
in  Fig.  31.  Such  a  molecule  has  nine  fundamental  vibra¬ 
tional  modes.  Seven  modes  are  in-plane  A’  species  and 
two  modes  are  out-of-plane  A"  species.  All  modes  are 
active  in  the  infrared  and  Raman  spectra,  McC-raw  et_  al. 
compared  the  infrared  absorption  spectrum  of  nitric-acid 
vapor  with  results  from  the  normal  coordinate  analysis. 

They  made  assignments  of  the  nine  fundamental  bands  and 
of  several  overtone  and  combination  bands.  Their 
assignments  of  the  fundamental  bands  for  nitric-acid 
vapor  are  valid  yet  today  and  are  presented  in  the  first, 
second,  and  fourth  columns  of  Table  VI. 

The  assignment  of  infrared  and  Raman  vibrational 

bands  for  aqueous  solutions  of  nitric  acid  is  more  difficult 

for  at  least  three  reasons.  One,  the  HNO?  monomer  per 

se  probably  dees  net  exist  in  the  hydrogen-bended  structure 

of  the  liquid.  Two,  the  liquid  structure  is  not  completely 

understood.  And  three,  because  KNG-;  dissociates  in  aqueous 

solutions.  A  graph  of  a  the  degree  of  dissociation  of  HNO, 

2U  p'V 

in  aqueous  solutions — as  a  function  of  molar  con¬ 
centration  M  is  shewn  in  Fig.  32.  In  aqueous  solutions 
the  dissociation  is 


v 11  j  ---  -r  —  ,-^m  ^ - 


Table  VI.  Previous  band  assignments  for  nitric-acid  vapor  and  liquid. 


Mode  (Species) 

Vapor 

Liquid 

Assignment 

v1(A’) 

3550.0 

3410 

HD 

stretch 

v2CA’) 

1708.2 

1675 

no2 

antisym.  stretch 

v,(A’ ) 

J 

1330.7 

1395 

HON 

bend 

v4(A’) 

1324.9 

1303 

no2 

sym,  stretch 

v5(A») 

878.6 

926 

mo2 

deformation 

vgCA') 

646.6 

677 

NO' 

stretch 

V-t(A'  ) 

/ 

579.0 

612 

ONO' 

bend 

c 

CO 

> 

762.2 

771 

nc2 

rock 

Vg(A") 

455.8 

485 

HO 

torsion 

Table  VII.  Molar  concentration  of  nitric  acid  solutions,  volume  density 

p  in  g/cm^,  number  of  water  molecules  per  acid  molecule  N /N„ 

vv  a. 

prior  to  dissociation,  degree  of  dissociation  a,  and  molar  con¬ 
centration  of  HNOo,  H0^+,  NO-3-,  and  K2O  after  dissociation  in 
the  nitric  acid  solutions. 


Solution 

p 

V*a 

0 

HNO3 

HOf1- 

HO,” 

h2o 

0.5  M 

1.015 

109.3 

0.99 

0,005  M 

0.495  M 

0.495  M 

54.16  M 

1.0 

1.032 

53.9 

0.97 

0,03 

0,97 

0.97 

52.85 

2.0 

1.065 

26.1 

0.95 

0.10 

1.90 

1.90 

50.27 

4.0 

1.130 

12.2 

0.35 

0.59 

3.41 

3.41 

45, 31 

8.0 

1.249 

5.2 

0.61 

3.12 

4.88 

4.88 

36.48 

15.7 

1.412 

1.5 

0.165 

13.11 

2.59 

2.59 

20.97 
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HN03  +  H20  t  H30+  +  no3".  (5) 

In  pure  liquid  nitric  acid  (v  24  M)  dissociation 
occ’urs  by  means  of  the  autoionization  process 
(or  ionic  self-dehydration  process) 

2(HN03)  t  H20  +  N 02+  +  N03"  ,  (6) 

thus  accounting  for  the  greater  electrical  conductivity 
of  the  pure  acid  compared  to  that  of  aqueous  solutions 
of  very  high  nitric-acid  concentration.  The  molar 
concentrations  of  HN03,  K30+,  N03~,  and  H20  after 
dissociation  in  0.5,  1,  2,  4,  8,  and  15,7  M  aqueous 
solutions  of  nitric  acid  are  presented  in  Table  VII. 

The  concentrations  of  products  of  the  dissociation  were 
determined  by  use  of  Eq,  (5)  and  the  dissociation  cun/e 
shown  in  Fig.  32. 

Aqueous  solutions  of  intermediate  nitric-acid 

pp  ph  / 

concentration  (8-20  M),  however,  have  Raman  and 

p/f  P'7/ 

infrared  spectra  closely  similar  to  that  of  the 

vapor.  In  the  spectra  of  the  liquids  there  are  also 

other  vibrational  bands  observed  which  correspond  to 

the  dissociation  products  and  to  the  liquid  water. 

2?/ 

Previous  assignments ,  by  McC-raw  et  al . of  fundamental 
vibrational  bands  based  on  Raman  spectra  of  anhydrous 
nitric  acid  (v  24  M)  and  aqueous  solutions  containing  95% 
(22.5  M)  and  60%  (13  M)  nitric  acid  are  shown  in  the  first, 
third,  and  fourth  columns  of  Table  VI. 

t-  c.  R(v)  and  N(v)  Spectra  of  MO,  Solutions. 

Experimental  methods  for  acquisition  of  the  near- 

r 


I 
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normal  incidence  relative  infrared  reflectance  spectra 

of  the  nitric-acid  solutions  were  similar  to  those 

28/ 

described  in  a  previous  paper, —  The  relative 
reflectance  spectra  thus  obtained  are  presented  in 
Pigs.  33  and  3^.  Relative  spectral  reflectance  R(v) 
is  defined  as  R(v)  =  R(v)_/R(v),,,  where  R(v)_  is  the 
spectral  reflectance  of  the  nitric-acid  solution  and 
R(v)w  is  the  spectral  reflectance  of  deionized  and 
filtered  water.  Standard  deviations  based  on  three 
independent  measurements  of  R(v)  were  generally  +  0,005  R(v) 
but  increased  to  about  +  0,01  R(v)  in  spectral  regions 
where  atmospheric  water  vapor  and  carbon  dioxide  are 
strongly  absorbing.  The  temperature  of  the ’water  and  the 
nitric-acid  solutions  was  about  27°C. 

The  relative  reflectance  spectra  were  phase-shift 

pQ  pQ  / 

analyzed  by  use  of  Kramers-Kronig  techniques' =— 1  ■  to 

obtain  spectral  values  of  N(v)  which  are  presented 

graphically  in  Pigs,  35-38,  Uncertainties  for  N(v) 

5  29/ 

obtained  by  these  methods^  are  typically  as  follows: 

+  0.01  n(v),  +  0,02  k(v)  for  k(v)  greater  than  0.05, 
increase  to  about  +0,2  k(v)  for  k(v)  =  0.01,  and 
I  1  become  much  larger  for  k(v)  less  than  0.01, 

d.  Discussion  of  Spectra  of  HNO,  Solutions, 

j  -> 

Water.  Temporal  clusters  composed  of  several 
■  hydrogen-bonded  H2O  molecules  constitute  the  structure 

of  liquid  water,  An  infrared  spectrum  of  liquid  water 
is  therefore  interpreted  in  terms  of  bands  associated 

with  Intramolecular  vibrations  of  Kp0  molecules  and 

*  ^ 
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Figure  33 


.  Near-normal  incidence  relative  reflectance  in 
the  spectral  region  310  to  2,500  cm"*,  for 
0.0,  1,  2,  4,  8,  and  15.7  M  aqueous  solutions 
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Figure  34.  Near-normal  incidence  relative  specular  reflectance 

R(v)  =  R(v)_,/R(v),r  in  the  spectral  region  2,500  to 
s  w 

5,000  cm-1  for  0,5j  1,  2,  4,  8,  and  15.7  M  aa.ueous 


solution  of  nitric  acid. 


0  500  1000  1500  2000  2500 

WAVE  NUMBER  (cnrr1) 


Figure  35.  Index  of  refraction  n(v)  in  the  spectral  region 

310  to  2,500  cm-1  for  0.5,  1,  2,  4,  8,  and  15.7  M 
aqueous  solutions  of  nitric  acid:  n(v)  is  the 
real  part  of  N(v). 
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Figure  36.  Index  of  refraction  n(v)  In  the  spectral  region 

2,500  to  5,000  cm-1  for  0.5,  1,  2,  4,  8,  and  15.7  M 
aqueous  solutions  of  nitric  acid:  n(v)  is  the  real 
part  of  N(v) . 
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'lgure  37.  Extinction  coefficient  k(v)  in  the  spectral 
region  310  to  2,500  cnT1  for  0,5,  1,  2, 


8  and  15.7  M  aqueous  solution  of  nitric  acid 
k(v)  is  the  imaginary  part  of  N(v), 
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Figure  38.  Extinction  coefficient  k(v)  in  the  spectral  region 
2,500  to  5,000  cm-1  for  0.5,  1,  2,  A,  3,  and  15.7  M 
aqueous  solutions  of  nitric  acid:  k(v)  is  the 
imaginary  part  of  N(v). 
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intermolecular  lattice  vibrations  of  the  clusters. 

There  are  three  intramolecular  vibrations :  the  v- 

syimetric  and  v,  antisymmetric  0-H  stretching  motions, 

and  the  v 2  O-H-O  bending  motion.  The  v2  band  has  a 

central  frequency  of  16^0  am-",  a  half-width  of  110  am-", 

and  a  band  strength  of  (4, 5+0, 8)  x  10-^  cnfYmolecule , 

The  \>^  and  stretching  modes  and  the  first  overtone 

2\>2j  which  may  be  in  Fermi  resonance  with  \g  or  v,, 

produce  bread  overlapping  bands.  The  v, ,  v-,,  and  2v2 

bands  are  therefore  grouped  together  in  order  to  analyze 

the  infrared  spectrum  of  liquid  water.  Together  these 

three  bands  form  a  broad  composite  band  of  central 

frequency  ,  half-width  r,  and  band  strength 

as  recorded  in  Table  VIII, 

The  intermolecular  vibrational  modes  are  not  well 

known.  However,  one  structural  model  proposed  by  Bandekar 
80/ 

and  Cumutte^  consists  of  a  central  H20  molecule 

hydrogen  bonded  with  C2,_  point-group  symmetry  in  a 

temporal  cluster  of  k  first  and  12  second  nearest 

neighboring  H20  molecules.  A  normal  coordinate  and 

’0/ 

Monte  Carlo  analysis  of  this  structure,—  neglecting 
the  Intramolecular  motion  and  degeneracies  yielded  six 
principal  lattice  vibrational  modes  with  a  total  of  30 
vibrational  species.  The  principal  lattice  vibrations 
were  75,  165,  219,  ^50,  550,  and  722  cm-",  Spectra 
obtained  during  the  present  investigation  include  only 
the  ^50,  550,  and  720  cm-1  modes.  These  three  principal 
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Table  VIII.  Measurable  parameters  for  Intramolecular  bands  of  selected 

nitric-acid  solutions:  is  the  central  frequency,  r  is 

the  half-width,  k(v)__v  is  the  value  of  k(v)  at  v  .  and  S. 

max  max  D 

is  the  band  strength. 


Solution(M) 

Molecule 

Band 

vCcm"1^^ 

r(cm_1) 

k(v4ax 

( cirr/molecule ) 

15.7 

hno3 

v2(A') 

1672+4 

47+8 

0.292 

1.70+0.l4xl0"21 

15.7 

HNO-3 

v,(A') 

1429+2 

98+6 

0.294 

2.0  +0.1xl0"2i 

J 

j 

15.7 

HNO., 

v4(A') 

1304+2 

50+5 

0.546 

4.7  +0.1xlC“21 

15.7 

HNO, 

v-(A') 

949+2 

49+1 

0.361 

1,6  lO.lxlO-21 

j 

D 

15.7 

HN03 

v6(A') 

691+2 

26+6 

0.256 

0.22+0.03xl0"21 

15.7 

HNCU 

v?(A») 

634+2 

0.283 

— 

15.7 

HNO  n 

Vq(A") 

778+2 

15+2 

0.259 

0,14+0, 04xl0"21 

j 

0 

15.7 

HN03 

2v4 

2633120 

— 

0.082 

: 

15.7 

HN03 

2v3 

2935+20 

— 

0.071 

15.7 

N03" 

^2 ( A2 ) 

820+3 

35+10 

0.198 

15-7 

NO," 

VE) 

730+7 

— 

0.5 

NO  " 

v-(E) 

1350+2 

139+10 

0,065 

l6.7+2.0xl0"21 

J 

J 

1.0 

NO" 

v3(E) 

1350+2 

157+10 

0.099 

I8.1+2.0xl0"21 

9.0 

h3o+ 

v2(A1) 

1210+40 

— 

8.0 

H30+ 

v4(E) 

1742+15 

180+20 

0.148 

4,0+0,5xl0"“~ 

0.0 

H20 

»v3>2v2 

3395+ 

390+ 

0..297 

3. 7+0. lxlO"* 1 

0.0 

h2o 

v2 

1640+ 

110+ 

0.137 

0.^5+0.08xl0"2- 

0.0 

Ho0 

VT 

580+ 

500+ 

0.443 

7.2+1. 5x1 0"2i 
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modes  provide  unresolved  broad  overlapping  bands  that 
are  collectively  referred  to  as  the  libraticnal  band 


v^.  parameters  for  the  v,  band  are  noted  in  Table 


VIII. 


Nitric-Acid  Solutions,  The  molar  concentrations  of 
molecular  constituents  HNO,,  R^0+ ,  NO,”,  and  f^O 
for  the  nitric-acid  solutions  are  listed  in  Table  VII, 

In  the  spectra  for  the  0,5,  1,  2,  and  4  M  solutions  we 
primarily  observed,  in  Figs,  33-38,  spectral  features 
due  to  liquid  water,  In  Figs,  33  and  34  liquid  water 
by  definition  has  a  relative  reflectance  R(v)  =1.0 
throughout  the  spectral  range.  The  bands  are  easily 
seen  in  the  spectra  shown  in  Figs.  37  and  38,  There 
were  also  secondary  spectral  features  due  to  the 
strongest  Intramolecular  vibrational  modes  of  NO  ”  and 

h3o+. 

The  NO^-  ion  is  planar  with  point-group 
symmetry,  and  thus  possesses  three  infrared  active 
vibrational  mode  v,(A,),  v^CE),  and  v^(E)  at  about 
820,  1350,  and  730  cm-*,  respectively.  The  strongest 
vibrational  modes  is  v^CE)  which  is  doubly  degenerate. 
The  v,(E)  mode  is  quite  evident  in  the  spectra  for  the 
0,5,  1,  2,  and  4  m  solutions.  Although  the  v^CE)  mode 
is  doubly  degenerate  it  appears  in  the  present  spectra 
as  a  doublet.  Coordination  of  NO,-  with  H-0  molecules 
is  probably  responsible  for  removing  the  degeneracy, 
Parameters  r,  and  measured  for  the  v,(E)  mode 
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of  N0^~  in  the  0,5  and  1  M  solutions  are  listed  in 
Table  VIII,  The  v2(A2)  and  v^(E)  modes  of  NC^~  are 
quite  weak.  They  were  not  observable  in  previous  infrared 

h/ 

reflectance  spectra^-  of  0.5  M  solutions  of  NaMO,, 

The  and  \>^(E)  modes,  however,  were  observed 

during  the  present  investigation  in  the  spectra  for 
the  8  and  15.7  M  solutions,  e,g.,  see  Fig.  38  for 
very  weak  bands  at  820  and  730  cm-1. 

The  H^0+  ion  is  pyramidal  with  CL  point-group 
symmetry,  and  thus  possesses  four  infrared  active 
intramolecular  vibrational  modes  v^CA-,),  v2  (A-, ) , 
v^(E),  and  v^(E).  The  v^(A^)  and  v^(E)  bands  are  very 
broad  and  extend  from  about  2650  to  338O  cm-1.  They 
were  not  observed  specifically  in  the  present  spectra 
for  the  nitric  acid  solutions  because  the  v., ,  v^,  and 
2\>2  modes  for  Yip  occur  in  the  same  spectral  region. 

The  v^CA-, )  band  of  occurs  at  a  central  frequency 
vmax  =  ^10+^0  cm~’L  ^  spectrum  for  the  8.  M  solution, 
and  is  best  observed  as  a  low- frequency  shoulder  on 
the  vjjCA')  band  of  HNO^  which  is  discussed  in  the 
following  paragraph.  The  presence  of  the  v2(A^)  band 
is  also  notable  in  the  spectra  for  all  the  other  solutions, 
e.g.  see  Fig.  38.  The  doubly  degenerate  v^(E)  band  of 
H,0+  occurs  as  a  shoulder  on  the  high  frequency  side 
of  the  v2  band  of  Yip.  Parameters  ,  r,  and  for 
the  v2(A^)  and  v1|(E)  bands  of  H^0+  are  listed  in  Table  VIII, 
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The  structure  of  the  HNO,  molecule  was  discussed 

J 

in  Section  V.A.4.b,  Intramolecular  vibrational  bands 
associated  with  HNO,  were  particularly  evident  in  the 
spectra  for  the  8  and  15.7  M  solutions.  Of  the  nine 
intramolecular  bands  only  the  v-,(A’)  and  vQ(A")  bands 
were  not  clearly  identified.  Both  of  these  bands  occur 
in  spectral  regions  where  liquid  water  also  has  bands, 

—1 

The  Vy(A')  band  at  central  frequency  =  634  +  2  cm” 
was  the  weakest  band  of  HNO^  that  was  observed.  The 
Vy(A’ )  band,  however,  was  net  resolved  to  the  extent 
that  its  band  strength  could  be  measured.  Parameters 
N^max*  r,  and  S^  for  the  six  other  bands  were  measured 
and  appear  in  Table  VIII. 

The  band  strengths  were  determined  from  graphs  such 
as  those  in  Fig.  38  by  manual  decomposition  of  overlapping 
bands  when  necessary,  by  subtracting  the  continuum 
absorption,  and  by  use  of  a  manual  planimeter.  The 
strength  of  the  v^(E)  band  of  H^C+  was  determined  by 
planimetering  only  the  high-frequency  half  of  the  band 
and  multiplying  the  result  by  a  factor  of  2.  The  strength 
of  the  v?(A')  band  of  HNO.,  was  determined  by  planimetering 
the  combined  areal  contribution  of  the  v^CA')  band  of 
HNO,,  the  \>2  band  of  f^O,  and  the  v^CE)  band  of  H,0+, 

The  band  strengths  due  to  of  ^0  and  v^(E)  of  K,0+ 
at  their  appropriate  molecular  density  were  then  subtracted 
from  the  measurement  of  the  combined  band  strength.  The 
resultant  band  strength  for  v2(A')  of  HN0,  was 
Sb  =  (1.70+0.14)  x  10“  ^  cm2 /molecule . 
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5. 


Organic  Liquids  and  Solids 
a.  Objectives 

Clcuds  ccmposed  cf  water  droplets,  which  are 
stabalized  against  evaporation  by  surface  coatings 
of  long-chain  alcohols,  clouds  composed  of  oil  droplets, 
and  clouds  composed  of  small  acrylic  particles  are  of 
practical  interest  as  obscurants  in  the  infrared 
spectral  region.  Radiation  transport  through  such 
clouds  can  be  computed  by  use  of  classical  or  modified 
Mie- scattering  theory.  The  spectral  refractive  index 
N(v)  =  n(v)  +  ik(v),  and  a  particle-size  distribution, 
of  the  cloud's  constituients  are  required  to  begin 
such  computations.  Values  of  N(v)  in  the  infrared 
are  available  for  water,  but  values  of  N(v)  for  oils, 
alcohols,  and  an  acrylic  were  net  previously  available. 
The  objective  of  the  investigation  reported  in  this 
section  was  to  obtain  the  near-normal  incidence  relative 
reflectance  spectra  R(v)  and  values  cf  N(v)  for  three 
oils,  three  solid  alcohols,  and  polymethy 1-methacrylate. 

Oils  selected  for  the  investigation  were  commercial 
diesel  fuel,  No.  100  pale  oil  (light  fog  oil),  and 
tetramethyl-tetraphenyl-trisiloxane  (TMTPTS),  The  TMTFTS 
is  commercially  available  as  DON  Coming  No,  70^  Oil, 

In  addition  to  its  ’use  as  an  obscurant,  diesel  fuel 
was  chosen  because  it  is  a  common  aerosol  which  is  pre¬ 
sent  in  the  atmosphere  due  to  incomplete  combustion  in 
ir.t emal-ccmbus t ion  engines.  Mo,  100  pale  oil  was  chosen 
because  it  is  used  by  the  Military  as  a  screer.ing-fcg 
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material.  TMTPTS  was  chosen  because  we  wished  to  extend, 

Into  the  infrared  spectral  region,  seme  previous  investiga- 
81  82/ 

tlons-  ->■  ■  of  its  optical  properties  in  the  ultraviolet 
spectral  region. 

Long-chain  alcohols  selected  for  investigation 
were  1-hexadecanol  (palmityl),  1-octadecar.ol  (stearyl), 
and  1-eicosar.ol  (arachidyl),  The  melting  temperatures  of 
these  alcohols  are  49,  58,5,  and  65.5°C,  respectively.  Their 
densities  are  0.8176,  0,8124,  and  0,8405,  respectively,  These 
particular  alcohols  were  of  interest  because  of  their 
use  in  water- fog  generating  equipment,  wherein  they  are 
melted  in  water  at  about  95°C.  Compressed  air  bubbling 
vigorously  through  the  water-alcohol  mixture  in  the  generator 
then  produces  steam-like  clouds  composed  of  alcohol-coated 
water  droplets.  As  the  droplets  cool  in  the  atmosphere 
the  alcohol  solidifies,  thus  forming  irregularly  shaped 
particles  which,  depending  on  the  ambient  temperature, 
possesses  a  liquid  water  or  ice  core. 

Polymethyl-methacrylate  (PMMA)  was  chosen  because  it 
is  a  readily  available  commercial  acrylic. 

Also,  of  basic  interest  was  an  interpretation  of 
the  R(v)  and  N(v)  spectra  of  these  materials  in  terms  of 
intramolecular  vibrational  modes, 

b.  Background  Information  for  Discussion  of  Organics 
Diesel  feel  is  composed  of  a  host  of  different 
hydrocarbons.  Its  primary  const ituient s ,  however,  are 
paraffins  £^^+2*  wkere  n  =  12-16.  No,  100  pale  oil  has 
the  consistency  of  a  light  lubricating  oil  and  is 
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probably  composed  primarily  of  paraffins  where  n  =  15-20, 

The  long-chain  alcohols  are  closely  related  in  their 

chemical  composition  and  structure  to  the  paraffins. 

Their  chemical  formulas  are:  1-hexadecanol  CH,(CHo)  CH-OH, 

o  d  Hi  d 

1-octadecanol  CH-(CH0)  CPL,0H,  and  1-eicosanol 
3  2  16  2 

CH,(CH0)  CHo0H, 

3  2  18  2 

In  1963  Snyder  and  Schacht Schneider  published  two 
3V 

papers—  on  vibrational  analysis  of  the  n-paraffins. 

A  synopsis  of  their  results  appears  in  Table  IX, 

Paraffins  CnH,^+2  with  n  odd  possess  C2v  point  group 

symmetry  and  those  with  n  even  possess  C2h  point  group 

syrmetry.  Infrared  active  species  of  C2v  are  A-,,  B-, ,  and 

B2  and  of  C2^  are  Au  and  B  .  The  assymimetric  and 

symetric  stretching  modes  of  CH2,  and  CH^  to  a  lesser 

extent,  in  the  2,929-2,849  cm”'*'  region,  and  the  HCH 

angular-bending  modes  of  CH2  and  CH^  in  the  1,376-1,473  cm”** 

region  were  of  primary  interest  in  this  investigation 

of  diesel  fuel,  fog  oil,  and  the  alcohols. 

PMMA  is  the  first  (k  =  1)  in  the  chemical  series  of 

polymethacrylates  [  CH2C(CH,)  (C02  P^k+l^  •  ITie  point 

n 

group  symnetry  of  PMMA  varies  with  structure  which  can 
be  either  isotactic,  syndictactic,  or  atactic.  Additionally, 
the  structure  can  vary  from  one  form  to  >_  -.her  in  a 

particular  PMMA  sample.  The  infrared  spectra  and  structure 

34/ 

of  PMMA  remains  as  an  open  research  topic,—  and  assign¬ 
ments  of  vibrational  species  for  infrared  bands  of  PMMA 
could  not  be  made  during  this  investigation.  However, 
most  of  the  infrared  bands  of  PMMA  are  characteristic  of 


various  chemical  groups  shown  in  the  preceding  chemical 
formula. 


Table  IX.  Infrared  active  vibrational  modes,  species,  and  frequencies  of 

n-paraffir.s  compiled  from  the  normal  coordinate  analysis  reported 
in  Ref.  (33). 


Vibrational  mode 

Species* 

n-odd  n-even 

Frequency  of  Frequency 
Range  (cm-1) 

CH^  asym.  C-K  stretch 

A1*B1 

2967 

CH0  asym.  C-H  stretch 

B2 

Au 

2912-2929 

CK^  sym.  C-H  stretch 

Ai>Bi 

B 

u 

2884 

CH2  sym.  C-H  stretch 

A1,B1 

Bu 

2849-2861 

CHo  out-of-plane 
*HCH  bend 

B2 

A 

u 

1465 

CH-,  in-plane  <HCH  bend 

A1’B1 

Bu 

1446-1473 

CH2  <HCH  bend 

A1’31 

Bu 

1446-1473 

CH^  sym.  <HCK  bend 

A1,B1 

Bu 

1376 

CH2  wagging 

Al’^l 

3u 

1174-1411 

C-C  stretch 

A1’B1 

3u 

885-1132 

CH,  rocking 

A1*B1 

3u 

895 

<C— C— C  bend 

A1’B1 

Bu 

0-535 

CH2  twist-rock 

Bo 

Au 

1175-1310 

CH2  rock-twist 

B2 

\l 

719-1060 

CH,-CH~  torsion 

B2 

Au 

0-153 

Crl-CH-  torsion 

d  2 

B2 

25-140 

^Paraffins  CrH2n+2  with  n  odd  belong  to  the 

C2v  point  group  and  these 

with  n  even  belong  to  t 

:he 

point  group. 

Infrared  active  species  of 

C^.  are  Ay,,  B, ,  32  and 

of  Cov, 
c.n 

are  A.  and  B, 

,.  All  species  of  C2y  are 

Raman  active,  but  only 

species 

i  Ag  and  Bg  0 

f  Coy,  are  Raman  active. 

! 
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The  chemical  formula  for  TPr^TS  is  (CH-0  (CgHe)  Si,0„, 
and  the  infrared  spectrum  of  this  liquid  is  characterized 
by  bands  due  to  the  methyl  (CH,)  groups,  the  phenyl 

J 

(CgH^)  groups,  Si-(CH^)  linkages,  3i  =  (CgH^)  linkages, 

Si-0,  and  many  structural  modes.  At  least  35  infrared 
active  bands  were  observed  for  this  material.  A  helpful 
review  of  infrared  bands  due  to  these  chemical  groups 
is  provided  by  References  (35)  and  (36). 

c,  R(v)  and  N(v)  Spectra  for  Diesel  Fuel,  Fog  Oil,  and 
TMTPTS. 

The  measured  near-normal  incidence  relative  reflectance 
spectra  for  diesel  fuel,  fog  oil,  and  TM’RTS (DC-704)  in 
the  spectral  regions  300-2,500  cm-"''  and  2,500-5,000  cm-* 
are  presented  in  Figs,  39  and  40,  respectively.  The 
relative  reflectance  spectrum  of  each  material  was 
analyzed  by  use  of  the  Kramers-Kronig  techniques  described 
in  Section  V.A.2.C.  to  obtain  N(v)  =  n(v)  +  ik(v)  for 
each  material.  Spectral  values  of  N(v)  for  diesel  fuel, 
fog  oil,  and  TMIPTS  are  presented  graphically  in  Rigs.  41-44. 

d.  Discussion  of  Spectra  for  Diesel  Fuel,  Fog  Oil,  and  TMTPTS. 

In  the  relative  reflectance  spectra,  the  spectral 

if  _1 

,  features  near  3,600  cm  and  1,640  cm  are  due  primarily 

to  liquid  water  which  was  the  reflectance  standard.  Liquid 
water  also  has  the  infrared  active  intermolecular 
vibrational  band  of  half-width  500  cm  centered  at  580  cm  . 
The  ether  features  in  the  R(v)  spectra  correspond  to 
infrared  active  bands  of  the  diesel  fuel,  fog  oil,  or 
TYT7TS ,  The  spectral  features  due  to  water  do  net  appear 

j  1 


WAVE  NUMBER  (cm-') 


Figure  39.  Near-normal  incidence  relative  specular  reflectance 
R(v)  in  the  300-2,500  cm”^  spectral  region  for  DC-71 
(IWTPIS),  fcg  oil,  and  diesel  fuel. 
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Figure  1*0,  Near-normal  Incidence  relative  specular  reflectance 
R(v)  in  the  2,500-5,000  cn-^  spectral  region  for 
DC-704  (TMTPTS),  fcg  oil,  and  diesel  fuel. 
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Figure.  41.  Index  of  refraction  r.(v)  in  the  300-2, 5CC  cm” 
spectral  region  for  DC-704  (TT'TTFrS),  fog  oil, 
and  diesel  fuel. 
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Figure  42.  Index  of  refraction  n(v)  in  the  2,500-5,000  cm-' 
spectral  region  for  DC-704  (TTITPTS),  fog  oil,  and 
diesel  fuel. 
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Figure  43.  Extinction  coefficient  k(v)  in  the  3CC-2.5C0  cm_1 
spectral  region  for  CC-TO^  (TJ-TTPrS),  fog  oil,  and 
diesel  fuel. 
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Figure  44.  Extinction  coefficient  k(v)  in  the  2,500-5,000  cm' 
spectral  region  for  DC-704  (TWTPTS) ,  fcg  oil,  and 
diesel  fuel. 
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in  the  M(v)  spectra,  However,  the  N(v)  spectra 
should  be  uniformly  smoothed  through  the  3,300-3,800 
cm-1  spectral  region  where  the  strong  absorption  due 
to  atmospheric  water  vapor  rendered  these  data  less 
reliable  than  in  adjoining  spectral  regions. 

The  spectra  for  diesel  fuel  and  fog  oil  were  quite 
similar  as  one  would  expect  from  prior  knowledge  of  their 
chemical  composition.  The  infrared  active  intramolecular 
bands  observed  in  the  spectra  for  diesel  fuel  and  fog 
oil  are  listed  in  Table  X.  The  molecular  vibrations 
corresponding  to  these  bands  are  listed  in  column  3  of 
Table  X,  and  were  obtained  by  consulting  Table  IX.  The 
bands  are  all  relatively  weak. 

Table  X.  infrared  active  intramolecular  bands  observed  in  the  spectra 
for  Diesel  Fuel  and  Fog  Oil. 
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There  are  at  least  35  infrared  active  bands  in  the 
spectra  for  TT4TFTS,  A  positive  identification  of  all 
the  bands  was  not  possible  without  a  normal  coordinate 
analysis  of  TMTPS.  However,  we  did  assign  the  more 
obvious  bands  by  use  of  References  (35)  and  (36). 

Our  assignments  were  tabulated  and  appear  in  Table  XI. 

We  have  not  determined  the  half-widths  and  strengths 
of  the  infrared  bands  of  these  materials,  but  we  will 
during  the  follow-on-contract  period. 


Table  XI.  Infrared  active  intramolecular  vibrational  bands  observed  in 
the  spectrum  of  tetramethyl-tetraphenyl-trisiloxane. 


v(cm-^) 

v(cm  "*■) 

3,056(W) 

Phenyl  C-H  stretch 

1,082  (M) 

Si-O-Si  str.  or  phenyl  in  plane  deform. 

3,012(W) 

Phenyl  C-H  stretch 

1,041 (Vs) 

Si-O-Si 

2,956(W) 

CH,  Asym.  stretch 

1,024 (Vs) 

Si-O-Si  or  phenyl  in  plane  deform. 

2,899(VW) 

CH,  sym,  stretch 

1,000(W) 

Phenyl  ring  breathing  (Star  of  David) 

2,853(VW) 

915(M) 

CH^  rocking 

1,961 (VW) 

Comb,  or  overtone 

840(S) 

(CH3)-Si-(CH3) 

1,877(W) 

Comb,  or  overtone 

791 (Vs) 

(CH3)-Si-(CH3) 

1, 58500 

Phenyl  skeletal  stretch 

774  (Vs) 

(CH3)-Si-(CH3) 

1,492(W) 

Phenyl  skeletal  stretch 

72 3 (Vs) 

Phenyl  out  of  plane  bend 

1,427 (M) 

Phenyl  skeletal  stretch 

693 (Vs) 

Phenyl  out  of  plane  skeletal  band 

1,409 (w) 

Asym.  CH^  deformation 

663(M) 

Si-(CgH5) 

1,375(W) 

Sym.  Ch^  deformation 

595  (M) 

Si-O-Si 

1,323(VW) 

570(M) 

Si-O-Si 

l,301(vw) 

468 (Vs) 

Si-(CgH5) 

1,256(S) 

Si-Ch,  &  Si-(CHg)  str. 

5  ->  2 

434  (Vs) 

(CK3)-Si-(C6H5)^ 

1,187  Or) 

Phenyl  in  plane  deform. 

374  (M) 

Si-IO-CCHg)]  2 

1,156(VW) 

Phenyl  in  plane  deform. 

351(M) 

(CH3)-Si-(Cg^5) 

1,115 (VS) 

Si-CgH^  stretch 

2 

Vs  =  very  strong,  S  =  strong,  M  =  medium,  W  =  weak,  VW  =  very  weak 
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Figure  ^5.  Near-normal  incidence  relative  specular  reflectance 
R(v)  in  the  300-2,500  cm-'*'  spectral  region  for 
1-hexadecanol  (polmityl) ,  1-octadecanol  (stearyl), 
and  1-eiconsanol  (arachidyl). 
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Figure  ^6.  Near-normal  incidence  relative  specular  reflectance 
R(v)  in  the  2,500-5,000  cm  *  spectral  region  for 
1-hexadecar.ol  (palmityl),  1-octadecanol  (stearyl), 
and  1-eicosanol  (arachiayl). 
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Figure  47.  Index  of  refraction  n(v)  in  the  300-2,500  cr.-1 

spectral  region  for  1-hexadecar.ol  (palrnityl), 
1-cctadecanol  (stearyl),  and  1-eicosanol  (arachidyl). 
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Figure  48.  Index  of  refraction  n(v)  in  the  2,500-5,000  cm-1 

spectral  region  for  1-hexadecanol  (palMtyl),  1-octadecanol 
(stearyl),  and  1-eicosanol  (arachidyl). 
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Figure  49.  Extinction  coefficient  k(v)  in  the  300-2,500  cm-1 
spectral  region  for  1-hexadecanol  (palmityl), 
1-cctadecar.ol  (stearyl),  and  1-eicosar.ol  (arachidyl). 


I 

J 


— - - 16- - - - 


_  i 


I 

I 


WAVE  NUMBER  (cm-1) 

Figure  50,  Extinction  coefficient  k(v)  in  the  2,500-5,000  cm-1 
spectral  region  for  1-hexadecanol  (palmityl) , 
1-octadecancl  (stearyl),  and  1-eicosanol  (arachidyl). 
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e. 


R(v)  and  M(v)  Spectra  of  for  Solid  Alcohols, 


The  measured  near-normal  incidence  relative 
reflectance  spectra  R(v)  for  1-hexadexcanol  (palmity), 
1-octadecancl  (stearyl),  and  1-eicosanol  (arachidyl) 
in  the  300-3,500  and  2,500-5,000  cm-"  regions  are 
presented  in  Figs.  45  and  46,  The  R(v)  spectrum  of 
each  alcohol  was  analyzed  by  use  of  Xraners-Krcnig 
techniques  to  obtain  N(v)  =  n(v)  +  ik(v)'  of  each  material. 
Spectral  values  of  N(v)  for  these  alcohols  are  presented 
graphically  in  Figs,  47-50, 

Discussion  of  Spectra  for  Solid  Alcohols. 

Ln  their  structure  and  chemical  composition  these 
alcohols  are  very  closely  related  to  the  n-paraffins,  in 
that  one  hydrogen  atom  of  the  n-paraffin  is  replaced  by 
an  0-H  member.  The  spectra  of  the  alcohols  thus  include 
many  features  that  are  similar  to  the  spectra  of  the 
n-paraffins.  The  strongest  bands  observed  in  these 
spectra  are  listed  in  Table  XEI,  The  number  of  vibrational 
modes  for  these  alcohol  molecules  is  9  n,  where  n  =  16, 

18,  and  20,  Thus  some  of  the  bands  listed  in  Table  XII 
are  actually  an  envelope  of  several  closely  spaced 
bands. 

Although  the  point  group  symmetry  of  the  alcohols 
is  net  the  same  as  that  of  the  corresponding  n-even  (C^) 
n-paraffin,  the  bands  in  the  alcohols  can  be  assigned 
by  using  References  (33)  as  a  guide. 
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Table  TIL  Wave-number  positions  and  assignments  for  the  strongest  infrared 
bands  observed  in  the  spectra  of  the  alcohols. 


1-hexadecanol 

(Palmityl) 

1-octadecanol 

(stearyl) 

1-eicosancl 

(Arachidyl) 

Vib rational 

Assignment 

3,223  cnfx 

3,223  cm"1 

3,228  cm-1 

C-K  stretch 

2,920 

2,920 

2,920 

CH2  asym,  stretch 

2,353 

2,853 

2,853 

CH2  sym,  stretch 

1,464 

1,464 

1,464 

-CH2-  deformation 

1,061 

1,063 

1,060 

C-C  stretch 

1,039 

1,038 

1,034 

C-C  stretch 

801 

801 

800 

CK2  rocking-twisting 

718 

719 

718 

CK2  rocking-twisting 

572 

578 

576 

C-C-C  structural  mode 

I  ♦ 


R(v)  and  N(v)  spectra  for  Pc lymethy 1-methacrylate 
(PMMft). 

The  measured  near-normal  incidence  relative  reflectance 
spectrum  for  PMiA  in  the  spectral  regions  300-2,500  cm-1 
and  2,500-5,000  cm-1  are  presented  in  Figs,  51(a)  and 
51(b),  respectively.  The  relative  reflectance  spectrum 
was  analyzed  by  use  of  the  Kramers-Kronig  techniques 
described  in  Section  V.A.2.C.  to  obtain  N(v)  =  n(v)  +  ik(v) 
for  this  material.  Spectral  values  of  N(v)  are  presented 
graphically  in  Figs.  52(a),  52(b),  53(a),  and  53(b). 


spectral  region  for  PMMA, 


for  PMMA 
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Discussion  of  Spectra  for  PMMA, 

The  chemical  formula  for  PMMA  is 


The  most  prominent  Infrared  active  bands  observed  in 
the  spectrum  of  PMMA,  and  their  assignments,  are  listed 
in  Table  XIII.  There  are  weaker  overtone  and  combination 
bands  which  are  not  listed  in  Table  XIII,  The  bands  at 
3,591  and  3,552  cm-1,  are  believed  to  be  due  to  a 
surface  layer  of  water  taken  on  by  the  PMMA.  No  attempt 
was  made  to  dry  the  PMMA  in  order  to  confirm  this 
assumption. 
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Table  XIIL  Wave-number  positions  and  assignments  for  infrared  bands 
observed  in  FMMA. 


v  ( cm”1 ) 

Assignment 

3.59KW) 

H20  asym.  stretch 

3,552(W) 

HjO  sym.  stretch 

2,99400 

CH2  Sc  CH^  asym.  stretch 

2,94900 

CH2  &  CH,  sym.  stretch 

l,730(s) 

C=0  stretch 

1,47900 

0-CH- 

1,44300 

C-CH^  asym.  deformation 

1,38500 

C-CH^  sym.  deformation 

1,26700 

C-O-C 

1,23900 

C-O-C 

l,198(s) 

C-O-C 

l,l47(s) 

C-O—  c 

1,062(W) 

C-C-C  stretch 

98500 

CH,  rock  or  twist 

966  (M) 

CHn  rock  or  twist 

917  (W) 

9 

840  (M) 

9 

833(W) 

9 

810  (W) 

9 

75000 

CH^  rock  or  twist 

63900 

O 

586  00 

C-CH,  wagging 

47600 

C-CH.,  deform.  +  C=0  deform.  +  ( 

35900 

C-O-CH^  deform.  +  C-CH^  deform 
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6,  Natural  Rutile 

a.  Objectives 

Clouds  composed  of  dust  particles  constitute 
naturally  occuring  obscurants.  Radiation  transport 
through  such  clouds  is  of  current  interest.  Thus, 
during  this  Grant  period  we  began  an  investigation 
of  the  optical  properties  of  natural  minerals.  The 
objective  of  this  investigation,  which  continues  during 
a  follow-on  contract,  is  to  measure  the  near-normal 
incidence  relative  reflectance  spectra  R(v)  of  selected 
natural  minerals,  and  to  compute  values  of  N(v)  by  use 
of  Kramers-Krcnig  techniques.  Natural  Rutile  (TiCU)  was 
one  of  the  minerals  selected  for  investigation. 

b.  Background  Information  for  Discussion  of  Rutile. 

The  infrared  reflectance  spectra  of  Rutile,  which  is 
an  optically  uniaxial  crystal,  has  been  the  subject  of 
several  prior  investigations.  The  more  significant 
of  these  investigations  can  be  found  in  References  (37) 
and  (38),  During  these  previous  investigations  Rutile 
crystals  were  cut  and  polished  with  faces  both  parallel 

37/ 

and  perpendicular  to  the  C-axis,  In  one  instance^ 
the  polished  crystals  were  etched  in  concentrated 
H,P0^  to  remove  an  anomalous  surface  layer.  In  the 
other  case  the  crystals  were  annealed  several  hours  at 
1700  K  to  regenerate  the  surface  after  polishing. 

During  our  investigation  a  Rutile  crystal  was  cut 
perpendicular  to  the  C-axis  and  polished.  The  crystal 
was  not  etched  or  annealed,  because  the  unetched  or 
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non-annealed  surface  nxcst  likely  corresponds  to  Rutile 

in  its  natural  environment.  Previous  investigations  of 
RQ/ 

Corundum^  demonstrated  a  decrease  in  reflectance  in 
the  spectral  region  of  phonon  modes  for  unetched  or 
nonannealed  crystals  in  comparison  to  the  reflectance 
of  crystals  which  were  etched  or  annealed.  The  reflectance 
of  the  unetched  crystals  can  diminish  to  50$  of  that 
for  the  etched  crystals  in  the  spectral  vicinity  of 
the  lattice  phonons,  This  is  precisely  the  effect  we 
observed  in  our  reflectance  spectrum  for  Rutile  when 
compared  to  the  reflectance  spectra  in  References  (37) 
and  (38), 

The  allowable  infrared  active  phonon  modes  for 
Rutile  are  A2  +  3Eu<  Additionally,  there  is  a  Raman 
active  mode.  The  Eu(T0)  modes  occur  at  ca.  189, 

382,  and  508  cm”*.  The  A2u(T0)  mode  occurs  at  ca.  172 
cm-1  and  the  A~  (LO)  mode  at  ca.  796  cm-1.  The  A, 

-*-0 

mode  occurs  in  the  500-600  cm-''-  region.  The  Eu  modes 
are  observed  when  the  incident  electromagnetic  radiation 
is  polarized  perpendicular  to  the  c-axis,  and  the  A0 
mode  when  the  radiation  is  polarized  parallel  to  the 
c-axis. 

c.  R(v)  and  N(v)  Spectra  for  Rutile. 

The  measured  near-normal  incidence  relative  reflectar.c 
spectrum  R(v)  for  natural  Rutile  in  the  300-5,000  cm”* 
spectral  region  is  presented  in  Fig.  5^.  The  relative 
reflectance  spectrum  was  analyzed  by  use  of  the 
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Kramers-Krcnig  techniques  described  in  Section  V.A.2.C. 
to  obtain  N(v)  for  this  material.  Spectral  values  of 
N(v)  are  presented  graphically  in  Figs.  55  and  5c. 

These  spectra  are  for  the  incident  electromagnetic 
radiation  polarized  perpendicular  to  the  C-axis  of 
the  Rutile  sample.  The  sample  was  net  etched  or 
annealed  after  it  was  polished. 

d.  Discussion  of  Spectra  of  Rutile. 

These  spectra  snow  to  of  the  strong  doubly  degenerate 
modes  at  503  and  400  cm-".  The  features  are 
best  observed  in  Fig,  56.  Additionally ,  there  is  a 
secondary  band  at  594  cm  which  has  previously  been 
attributed  to  the  A,  mode.  Although  the  A.  mode  is 

-Lfe  -*-0 

not  infrared  active  in  a  pure  Rutile  crystal,  it  is 
obser/ed  in  the  infrared  spectrum  of  natural  crystals 
due  to  the  presence  of  impurities,  She  A-  mode  could 

-*-o 

also  be  due,  in  part,  to  the  lack  of  etching  or  annealing 
of  the  polished  surface,  There  is  also  another  weak 
unidentified  mode  at  670  cm-". 

When  the  Kramers-Kronig  analysis  was  applied  to 
R(v)  for  Rutile,  small  negative  values  of  k(v)  were 
obtained  in  the  spectral  regions  (v  %  800  cm  where 
this  material  is  a  -weak  absorber.  In  plotting  the 

graph  of  k(v)  in  Fig,  56  the  negative  values  were 

\ 

supressed  thus  causing  the  unusual  appearance  of  the 
k(v)  spectrum  in  the  region  v  £  300  cm-  ,. 


S-rv*!* 
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At  the  time  these  data  were  acquired  we  did  not 
possess  an  adequate  polarizer  and  therefore  could  not 
obtain  the  spectrum  of  the  k^,  mode.  We  new  have  a 
polarizer  and  will  continue  this  investigation  during 
the  follow-on  contract. 

Limestone 
a.  Objectives 

The  advent  of  remote  sensing  by  use  of  optical 
instruments  mounted  in  ground  stations,  aircraft,  and 
satellites,  and  an  interest  in  the  thermal  budget  of 
the  planets,  has  produced  an  increasing  interest  in 
spectral  values  of  the  complex  refractive  index 
N(A)  =  n(A)  +  ik(A)  of  materials  commonly  found  on 
terrestial  and  other  planetary  surfaces.  The  present 
investigation  was  undertaken  to  obtain  spectral  values 
of  for  limestone  which  consists  primarily  of 

randomly  oriented  microcrystals  of  cal cite.  Other 
constituents  of  limestone  are  randomly  oriented  micro¬ 
crystals  of  dolomite  and  organic  fossil  material. 

Limestone  is  relatively  abundant  in  the  terrestial 
environment,  but  would  be  on  other  planetary  surfaces 
only  if  biological  conditions  as  we  know  them  on  earth 
cr.ee  existed  on  these  planets. 

In  the  present  investigation  spectral  values 
N(A)^  are  obtained  in  the  following  manner,  First  the 
near-normal  incidence  relative  reflectance  R(.\)  =  R(\)./R(\) 

i  fi 

of  limestone  l,  with  water  w  as  the  reflectance  standard, 
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was  measured  throughout  the  0.2-32.3-um  wavelength 
region.  Second,  a  Kramers-Kronig  phase-shift  dispersion 

1  liQ/ 

analysis— —  of  R(A)  provided  a  spectrum  Aij>(A)  which 

is  the  difference  in  phase  change  for  electromagnetic 

waves  of  wavelength  A  reflected  from  the  limestone 

and  from  the  water  standard.  An  algorithm-* Uy//  that 

uses  R(A),  A<j>(A),  the  measured  angle  of  incidence  8, 

and  known  spectral  values  for  the  complex  refractive 

index  N(A)  of  water  then  provided  scectral  values 
w 

of  N(A)  .  Additionally,  a  classical  Lorentz  dis¬ 
persion  analysis  of  the  absolute  reflectance  spectrum 
R(  A ) ^  of  limestone  yielded  spectral  values  of  the 
complex  dielectric  function  eCA)^  from  which  N(a) 
were  readily  obtained. 

The  optical  properties  of  natural  minerals  and 
rocks  have  been  investigated  by  several  other  scientists, 

ill/ 

Aronson  and  Strong—  obtained  Lorentz  parameters  in 

the  infrared  for  muscovite  mica,  anorthosite,  diopsidic 

pyroxenite,  almandite-pyrcpe ,  garnet,  and  soda  lime 

il?/ 

glass.  Pollack  et_  al. — •  obtained  Lorentz  parameters 
for  obsidian,  basaltic  glass,  basalt,  and  andesite, 

ill/ 

Spitzer  and  Kleinman— -  obtained  Lorents  parameters 

for  the  infrared  lattice  bands  of  quartz.  Holland 
4  4/ 

et  al. — •  measured  infrared  reflectance  spectra  in  the 
2-50-ym  wavelength  region  and  computed  N(a)  by  use  of 
Kramers-Kronig  techniques  for  quartz,  calcite,  dolomite, 
fluorite,  galena,  sphalerite,  brucite,  magnetite, 
goethite,  and  tite.  Analysis  of  the  optical 


properties  of  minerals  and  rocks  in  the  visible 
■and  infrared  are  also  provided  by  Refs.  (45)-(48), 

b.  Preparation  of  Samples. 

Three  samples  of  the  Bethany  Falls  limestone, 

iiq/ 

a  prominent  member  of  the  Swope  formation,—  were 
collected  from  the  approximate  center  of  Section 
7-R32W-T49N  in  northeast  Kansas  City,  Mo.  All  samples 
were  prepared  by  standard  polishing  techniques.  Follow¬ 
ing  a  final  buffing  with  0,1  urn  Linde  B,  the  samples 
were  examined  for  scratches  and  pits  by  use  of  a 
Vickers  metalore  microscope.  The  three  polished  lime¬ 
stone  samples  were  slightly  marred  by  pits  and  cracks 
due  to  the  presence  of  natural  vugs  and  joints.  Modal 
analyses  were  performed  by  point  counting  in  order  to 
assess  the  impact  of  the  pits  and  cracks  on  subsequent 
reflectance  spectra.  Modal  analyses  revealed  2.96  to 
7,89  areal  percent  cavities  with  a  mean  areal  percent 
of  4.60.  Consequently,  the  reflectance  spectra 
described  hereinafter  were  increased  by  a  factor  of 
1.04  tc  compensate  for  reduced  specular  reflectance 
due  to  the  presence  of  pits  and  cracks,  although 
this  may  be  a  function  of  wavelength  and  not  a  constant. 

The  Bethany-Falls  samples  consisted  of  dense 
fine-grained  buff  to  gray  mottled  limestone.  Calcite 
constituted  greater  than  99.7%  of  the  samples,  the 
remainder  was  sparse  pyrite.  Because  the  Bethany 
Falls  limestone  commonly  contained  dolomite  in  abundance, 
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this  mineral  was  specifically  searched  for.  However, 
results  from  both  stained  thin  sections  and  x-ray 
diffraction  disclosed  no  dolomite, 

Texturally  the  samples  were  all  very  similar, 
typically  containing  5-20  pm  irregular  interlocking 
grains  and  no  readily  apparent  preferred  orientation. 
Seme  grains  were  as  small  as  two  microns  and  others 
were  rare  and  widely  dispersed  clots  of  exceptionally 
large  ^0-200  -pm  grains  which  appeared  to  be  fossil 
fragments.  A  minor,  but  ubiquitous,  dusting  of  two 
urn  diameter  pyrite  grains  was  also  present.  Analysis 
of  polished  sections  disclosed  from  0.14  to  G.2 
percent  pyrite. 

c.  R(A)  and  N(A)  Spectra  for  limestone, 

In  the  present  investigation  we  measured  relative 
reflectance  R(A )  =  R(A)^/R(A)w,  where  l  and  w  denote 
limestone  and  water,  respectively.  Data  were  acquired 
in  the  following  order:  Z-W-L-V7-1-V/-L-V/-Z ,  where 
Z  indicates  a  spectral  scan  made  when  the  beam  of 
radiant  flux  was  blocked  just  in  front  of  the  infrared 
source,  W  indicates  a  spectral  scan  made  with  the  water 
standard  w  in  the  sample  position,  and  each  1  indicates 
a  sequence  of  three  scans,  one  for  each  of  the  three 
limestone  samples,  After  data  we re  acquired  throughout 
the  0,2-32.8  pm  spectral  region  the  relative  reflectance 
of  the  limestone  samples  was  calculated  as  follows. 
First,  the  average  of  the  two  spectral  scans  Z  was 
substracted  point  by  point  from  the  other  thirteen 


spectral  scans  denoted  by  W  and  L.  Second,  the 
three  spectral  scans  denoted  by  each  1  were  averaged 
and  ratioea  to  the  average  of  the  spectral  scans  W 
made  just  before  and  after  spectral  scans  L,  This 
procedure  provided  three  relative  reflectance  spectra, 
one  for  each  of  the  three  groups  of  three  spectral 
scans  L,  which  we re  then  averaged  to  provide  a  com¬ 
posite  spectrum  and  standard  deviations.  The  composite 
spectrum  was  then  multiplied  by  1,04  as  indicated  in 
the  preceding  Section,  The  absolute  reflectance 
R(A)^  of  the  limestone  was  provided  by  R(A)2  =  R(A)  R(A)  , 

where  R(A)  was  computed  by  use  of  the  Fresnel  equations 
2/ 

and  known  values—  for  N(A)  ,  and  R(A)  designates  the 

w 

composite  relative  reflectance  spectrum.  The  absolute 
reflectance  spectrum  R(A) ^  of  limestone  is  shown  in 
Figs.  57  and  58.  The  standard  deviations  were 
+  0.06  R(a)^  in  the  spectral  region  beyond  25  vm,  and 
reduced  to  about  +  0,01  RCA)^  in  the  spectral  region 
from  0.2  to  20  um. 

Evidence  of  birefringence  was  not  observed  as 
the  limestone  samples  were  rotated  while  placed  in 
the  ref lactometer;  thus  we  treat  limestone  as  a 
homogeneous  isotropic  material, 

Spectral  values  of  the  complex  refractive  index 
N(A),  were  obtained  by  ’use  of  the  Kraners-Krcnig 
algorithm  described  in  Refs,  (1)  and  (40).  Values  for 
N(A)a  are  presented  in  graphical  form  in  Figs,  59  and  60, 
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Figure  52.  Near-normal  Incidence  absolute  reflectance  R(v) 
in  the  17-32.8  pm  spectral  region  for  Bethany 
Falls  Limestone. 
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Figure  60.  Complex  refractive  index  N(A)  +  ikO)  in  the 
17-32.8  urn  spectral  region  for  Bethany  Falls 
Limestone. 
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Evaluation  of  the  Kramers-Kronig  dispersion  relation 

| 

i  j  requires  an  integration  of  R(A)  over  the  entire  electro¬ 

magnetic  spectrum.  Our  data,  however,  extended  only 
from  0.2  to  32.3  pm,  Thus  we  made  the  approximations 

I 

>-  that  R(A )  =  R(0.2  pm)  throughout  the  0  <_  A  <_  0.2-pm 

region,  and  R(A)  =  R(32.S  um)  throughout  the  32.8  <_  A  <_  °°-pm 
region.  The  latter  assumption  is  open  to  serious 
inquiry  because  calcite,  the  primary  constituent  of 
limestone,  possess  other  infrared  active  lattice  bands 
at  wavelengths  greater  than  32.8  pm. 

Spectral  values  of  N(A)£  were  also  obtained  by 
use  of  classical  Lorentz  dispersion  theory.  According 
to  the  Lorentz  theory  the  complex  dielectric  function 
e(v)  at  wavenumber  v  =  A~“  is  given  by 

4wp  .  v . 2  [ ( v . )+iy .  v.,  v ] 
e(v)  =  Ea+  Z  - -J  11  ^  ,  (1) 

j  (Vj2-V2)  +  (YjVjV) 

where  is  the  high  frequency  dielectric  constant, 
p,  is  the  strength  of  the  j  n  infrared  active  band, 

v.,  =  A.”  is  the  central  wave-number  position  of  the 

J  J 

j  '  Infrared  active  band,  and  y^  is  a  damping  coefficient. 
Lorentz  parameters  determined  for  limestone  appear  in 
Table  XU.  The  relation  between  N(  A)  ^  and  e(A)t  is 

e(A)^  =  H(A)£2  ,  (2) 

from  which  the  real  -and  imaginary  parts  of  N ( A ) }  =  nCA^+ikCA) 

i 

can  be  determined, 

I" 

I  I. 


125 


Table  XT/.  Classical  dispersion  theory  parameters  for  carbonate 
minerals:  Limestone,  Calcite,  and  Dolomite. 


(um)  v.(cm-1)  p. 


1 

6.55 

1517.5 

0,005 

0.020 

2 

6.83 

1464.1 

0.119 

0.031 

3 

7.045 

1419.4 

0.245 

0,025 

4 

11.47 

871.8 

0.052 

0.011 

5 

14.11 

708.7 

0.010 

0.0075 

6 

27.43 

364.6 

0.013 

0.050 

7 

33.80 

295.9 

1.220 

0.050 

8 

36.30 

275.5 

1.020 

0.092 

=  2.38, 

Limestone 

1 

7.11 

1406.5 

0.634 

0.0055 

2 

14.11 

708.7 

0,0183 

0.0046 

3 

28.95 

3^5.4 

0.034 

0.069 

4 

3^.50 

289.9 

2.5 

0.042 

=  3.00, 

Calcite  (EiC) 

1 

6.80 

1470.6 

0.025 

0.010 

2 

6.9^ 

1440.9 

0.650 

0.004 

3 

13.82 

723.6 

0.030 

0.012 

4 

14.22 

703.2 

0.0275 

0.007 

5 

25.00 

400.0 

0.009 

0.01C 

6 

27.00 

370.4 

0.075 

0.025 

<7 

28.10 

355.9 

2.220 

0.041! 

=  3-50,  Dolomite  (Ei_C) 
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A  comparison  cf  values  for  N(X)  obtained  by 
Kramers-Kronig  techniques  and  by  use  cf  the  Lorentz 
theory  shewed  reasonable  agreement  in  the  spectral 
region  from  0,2  to  about  25  um  and  increasingly  poor 
agreement  as  wavelength  increased  beyond  25  um.  Both 
the  real  n(X)£  and  imaginary  k(X)^  parts  cf  NCX)^  were 
higher  in  the  25  -  32.8  pm  region  for  the  Lorentz  theory 
than  for  the  Kramers-Kronig  techniques.  This  difference 
in  N(X)^  beyond  25  um  is  attributed  to  the  assumption 
that  R(x)  =  R(32,8  pm)  for  the  Kramers-Kronig  technique 
in  the  spectral  region  32, 3  <_  X  <_  °°  um.  The  values  cf 
N(x) „  in  the  region  beyond  25  um  as  computed  from  the 
Lorentz  theory  are  therefore  preferred  over  these 
obtained  from  the  Kramers-Kronig  technique, 

The  regions  of  anamalous  dispersion  shown  in  the 
spectrum  for  n(X)^,  and  strong  absorption  shown  in 
the  spectrum  of  k(X)^,  centered  at  wavelengths  7.00, 

11.4,  and  14,05  in  Pigs.  59  and  60  are  characteristic 
of  the  ordinary  and  extra-ordinary  rays  in  birefringent 
calcite  (CaCO^). 

8.  Natural  Waters 

I  ,  a.  Objectives 

The  complex  refractive  index  is  denoted  by  N  =  n  +  ik, 
where  italics  denote  a  complex  number,  n  is  the  index  of 
refraction,  k  is  the  extinction  coefficient,  and  i 
represents  the  square  root  of  minus  one,  A  quantitative 

1 

*  '  '  .  v.  . 

&  -  — —  -x.— - - - - 
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knowledge  of  spectral  values  of  N  for  natural  waters 

is  presently  of  considerable  interest  to  many  scientists. 

Such  knowledge  has  application  to:  (1)  computer 

simulations  of  the  transport  of  electromagnetic  radiation 

through  hydrcsols,  aerosols,  and  maritime  fogs, 

(2)  theoretical  investigations  of  the  energy  and  hydro- 

logical  balance  of  the  earth's  lakes,  oceans,  and 

atmosphere,  (3)  remotely  sensing  the  chemical  quality  of 

environmental  waters,— ^  and  (4)  solar-energy  development.— ^ 

i  4  52/ 

We  previously  reported-*-- -  measurements  of  the  relative 
reflectance  spectra  in  the  2-20-um  wavelength  region  of 
the  infrared  for  individual  aqueous  solutions  prepared 
with  distilled  water  in  the  laboratory.  The  solutions 
contained  NaCl,  rLSO^,  MH^P^PO^,  NaNO,,  or  NaKCO,. 

Kramers-Kronig  phase-shift  dispersion  analysis  of  the 
reflectance  spectra  for  those  solutions  provided  spectral 
values  for  N.  We  observed  that  monatomic  ions  such  as 
Wa+,  Cl”,  and  K+  perturbed  the  intra-  and  inter-molecular 
bonds  of  the  liquid-water  substances  sc  that  corresponding 
wavelength  positions  of  infrared  active  fundamental 
vibrational  and  libraticnal  (lattice  vibrational)  modes 
were  moderately  shifted  to  shorter  and  longer  wavelengths, 
respectively.  We  also  explicitly  identified  twenty-one 
infrared  active  vibrational  modes  associated  with  SC7~, 

WHa+,  HSCU”,  K,0+,  HjPOjj”,  MO,”,  HCO,~,  or  HjO,  The 
half  width  r,  central  spectral  position  of  each  mode  X, 
mode  (band)  strength  £>.,  and  maximum  value  of  k  were 


Table  XV  (cont 


tabulated  for  each  of  the  twenty-one  modes.  Ihese 


parameters  are  presented  collectively  for  the  first 
time  in  Table  XV  Band  strengths  were  evaluated  by 


use  of 


^i  =  *1  J"band  ^i  dv  *  (1) 


where  a^  (mole/liter)  was  the  molecular  density  of  the 

til 

i  constituent  in  the  solution,  Ak^  ’was  the  part  of  k 

th 

due  to  a  fundamental  vibrational  mode  of  the  iu  con¬ 
stituent,  and  ~  =  X"1  with  ~  as  the  wave  number  (crrf~) 
and  X  the  wavelength  (cm).  The  integrals  were  evaluated 
over  the  spectral  region  of  an  infrared  active  band  by 
planimetering  appropriate  areas  on  large  graphs  of  Ak 
vs.  ~,  In  Table  XV  relative  band  strengths  RS^  were 
compiled  based  on  the  arbitrary  choice  RS^  =  1,00  for 
the  combined  v-^CA-^v^B^)  modes  of  H20;  the  absolute 
band  strength  of  these  combined  modes  was  =  2.196 
(cm  moie/liter)-~.  v^(A^)  denotes  the  symmetric  OH 
stretch  of  the  H_0  monomer  and  v,(B.,)  denotes  the  ar.ti- 

d  Ox 

symmetric  CH  stretching  mode. 

The  investigations  reported  here  had  three  objectives : 
(1)  to  measure  the  relative  reflectance  spectra  in  the 
infrared  for  saline  natural  waters  containing  both 
monatomic  and  polyatomic  ions,  (2)  to  compute  and  tabulate 
spectral  values  of  N  for  these  samples,  and  (3)  to 
compare  the  relative  reflectance  spectra  for  the  natural 
waters  with  reflectance  spectra  for  aa.ueous  solutions 
containing  a  single  Inorganic  salt.  Natural  saline  waters 
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selected  for  investigation  were  obtained  from  the  San 

Francisco  3ay,  Pacific  Ocean,  Atlantic  Ocean,  Great 

Salt  Lake  (Utah),  Dead  Sea  (Israel),  and  an  effluent 

pit  adjacent  to  a  phosphate  nine  in  the  Alafia  river 

region  of  Central  Florida.  All  samples,  except  the  one 

from  the  Atlantic  Ocean,  were  supplied  by  members  or 

associates  of  the  U.S.  Geological  Survey,  EROS  Program 

Office.  Chemical  and  spectrographic  analyses  of  the 

samples  performed  by  chemists  at  USGS  laboratories,  and 
5*3/ 

by  Neev  and  Energy—  for  the  Dead-Sea  sample,  are 
presented  in  Table  XVI, 

b.  R(X)  Spectra  for  Natural  Waters. 

The  procedure  for  measuring  the  relative  reflectance 
spectra  of  the  natural  water  samples  was  identical  to 

5  2/ 

that  used  during  our  investigations  of  Nacl  in  water.— 
Relative  reflectance  R+  is  defined  as  R+  =  R+s/R+w ,  where 
the  +  subscript  denotes  radiant  flux  linearly  polarized 
with  the  e]ec trie- field  vector  perpendicular  to  the 
plane  of  incidence,  and  R+s  and  R+w  denote  absolute 
specular  reflectance  of  the  natural  water  sample  and  the 
distilled  water,  respectively.  Relative  reflectance 
spectra  throughout  the  2-20-um  region  for  the  six  water 
samples  investigated  are  shown  in  Figs,  61-63%  The  angle 
of  incidence  was  9  =  70,03C+G,23°.  Each  spectrum  shown 
in  Figs.  6 1-6 3  represents  smoothed  data  obtained  In  the 
following  manner.  Averages  <R+>  of  three  measurements  of 
E  were  determined  at  consecutive  0,05— uR  intervals 
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Figure  6l.  Measured  relative  reflectance  spectra  in  the  2-2C-um 

(5000-500-cirfx)  spectral  region  for  surface  water  sanples 
from  San  Francisco  Bay,  Atlantic  Ocean,  and  Pacific  Ocean. 
The  angle  of  incidence  was  70. 03°+0.23°.  Distilled  water 
was  the  standard  reflector;  its  relative  reflectance  is 
1.0  at  all  wavelengths.  The  radiant  flux  was  linearly 
polarized  with  the  electric-field  vector  perpendicular  to 
the  plane  of  incidence.  Standard  deviations  are  stated 
in  the  text. 
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Figure  62.  Measured  relative  reflectar.ce  spectra  in  the  2-20-ym 
(5000-500-cn-*)  spectral  region  for  surface  water 
samples  from  the  Great  Salt  lake  (Utah)  and  the  Dead 
Sea  (Israel).  The  angle  of  incidence  was  70,03°+0.S3°. 
Distilled  water  was  the  standard  reflector;  its  relative 
reflectance  is  1.0  at  all  wavelengths.  The  radiant 
flux  was  linearly  polarized  with  the  electric-field 
vector  perpendicular  to  the  plane  of  incidence. 

Standard  deviations  are  stated  in  the  text. 
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Figure  63.  Measured  relative  reflectance  in  the  2-20-ym  (5000-500-cnf 1) 
spectral  region  for  a  surface  water  sarcple  from  an 
effluent  pit  of  a  phosphate  mine  in  the  Alafia  River 
Basin  of  Central  Florida.  The  angle  of  incidence  was 
70.03°+0.23°.  Distilled  water  was  the  standard  reflector; 
its  relative  reflectance  is  1.0  at  all  wavelengths.  The 
radiant  flux  was  linearly  polarized  with  the  electric-field 
vector  perpendicular  to  the  plane  of  Incidence.  Standard 
deviations  are  stated  in  the  text, 
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throughout  the  2-20-wn  region.  The  standard  deviation 
for  <R+>  was  about  +0.008  <R+>,  In  narrow  spectral  regions 
centered  about  the  central  positions  of  the  fundamental 
Infrared  bands  of  atmospheric  water  vapor  at  2.75  and 
6.09  pm,  and  atmospheric  C02  at  4.3  and  14,7  ym,  the 
standard  deviations  increased  to  about  +0,025  <R,  > . 

Our  confidence  in  the  data  was  then  increased  significantly 
by  a  critical  examination  and  manual  smoothing  of  a 
point -by-point  plot  of  <R+>  vs,  wavelength  A, 

The  spectral  resolution  A/AA  was  from  125  to  200 
throughout  the  2-20-pm  region;  AA  represents  the  spectral 
slitwidth  of  the  monochromator,  Care  was  taken  to  keep 
the  air-sample  interface  clean  and  free  of  standing  waves. 
The  temperature  of  the  samples  was  about  300K°, 

c.  N( A)  Spectra  for  Natural  Waters 

Spectral  values  of  Ng  for  the  natural  water  samples, 
designated  here  by  the  subscript  s,  were  calculated  by 
•use  of  a  Kraners-Kronig  (KK)  algorithm  applied  to  the 
relative  reflectance  spectra  (Hale  et,  al,  1973,  Querry 
et.  al.  1974).  KK  analysis  of  a  relative  reflectance 
spectrum  provides  A<(>+  =  4>+u  -  <(>  the  difference  in 
phase  shifts  4>+  for  a  monochromatic  electromagnetic  wave 
reflected  from  s  the  air-sample  and  u  the  air-water 
interfaces ; 

A  00  m[<R(x)  ,>] 

A*(A0)  »Prin.  -m  /  —5 - 5 - 

0  v  - x 

(2) 


INDEX  OF 
REFRACTION 


WAVE  NUMBER  (cm-')  xIO5 


Figure  64.  Index  of  refraction  n,  in  the  2-20-um  (5000-500-cirf -1) 
region,  for  distilled  water  and  for  the  six  surface 
water  samples,  The  algorithm  consisting  of  Eqs.  (2)- 
(12)  was  used  to  compute  n  from  the  measured  relative 
reflectance  spectrum  for  each  surface  water  sample. 
The  uncertainty  is  about  +0.01n,  Data  from  Kale  and 
Querry  (1973)  was  used  to  plot  n  for  water. 
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Figure  65.  Extinction  coefficient  k,  in  the  2-20-wn  (5000-50C-cm_1) 
region,  for  distilled  water  and  for  the  six  surface 
water  samples.  The  algorithm  consisting  of  Eqs,  (2)- 
(12)  was  used  to  compute  k  from  the  measured  relative 
reflectance  spectrum  for  each  water  sample.  The  uncertainties 
are  about  +0,02k  for  k  greater  than  0,045,  increase  to 
about  +0.15k  for  k  =  0.01,  and  become  much  larger  for  k 
less  than  0.01.  Data  from  Hale  and  Querry  (1973)  was 
used  to  plot  k  for  water. 


I 

I 


•  *  — T-r - 

- - - v - - - 


Prin.  denotes  the  Cauchy  principal  value  of  the  integral. 

The  smoothed  reflectance  spectra  shown  in  Figs.  61-63 

were  used  to  compute  values  for  •  Use  of 

Eq,  (2)  requires  a  quantitative  knowledge  of  the  relative 

reflectance  spectrum  <R(x)+>  for  all  X,  Measurements, 

however,  were  made  only  in  the  2-2 C- pm  region.  Because 

<R(X)+>  was  the- relative  reflectance  of  two  aqueous 

substances,  rather  than  the  absolute  reflectance,  only  a 

small  negligible  error  in  A<Kxo)  was  introduced  by 

assuming  for  each  spectrum  that  <R(X)+>  was  equal  to 

<R(2  pm)+>  and  <R(20  pn)+>  throughout  the  X  <_  2  pm  and 

X  <_  20  pm  regions,  respectively. 

The  refractive  index  N  =  n  +  ik_  of  the  saline 

water  samples  were  then  determined  by  use  of  the  algorithm 

from  Querry  et  al.,-^  and  by  use  of  N  =  n  +  ik  for 

2/ 

distilled  water  as  tabulated  by  Kale  and  Querry.—  The 
angle  of  incidence  was  9  =  70,03°.  The  uncertainties 
were  about  +0,01n  ,  about  +0.C2k  for  k  greater  than 
0.045,  increasing  to  about  0,15ks  for  ks  =  0,01,  and 
became  much  larger  for  k_  less  than  0.01.  Graphical 
presentations  of  n  and  k_  vs.  X  appear  in  Figs,  64  and 

S  S 

j  65 j  respectively, 

d.  Discussion  of  R(x)  Spectra  for  Natural  Waters 

The  relative  reflectance  spectra  for  the  six  natural 
water  samples  are  compared  here  with  absolute  and  relative 
reflectance  spectra  previously  obtained  in  cur  laboratory 
for  individual  1,  3,  or  5M  aqueous  solutions  of  NaCl 
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Figure  66.  Measured  relative  reflectance  (lower  panel)  and  calculated 
absolute  reflectance  (upper  panel)  spectra  in  the  2-2C-yin 
region  for  aqueous  solutions  of  1,  3,  and  5M  NaCl  content. 
The  angle  of  incidence  was  70.03°+0,23°.  Distilled  water 
was  the  standard  reflector;  its  relative  reflectance  is 
1.0  at  all  wavelengths.  The  radiant  flux  was  linearly 
polarized  with  the  electric-field  vector  perpendicular 
to  the  plane  of  incidence.  The  uncertainty  in  the  relative 
reflectance  is  about  +0,008  <R+>,  Reproduced  by  permission 
of  M.  R.  Querry  et.  al. ,  J.  Dpt.  Soc.  An.  62,  949-855, 

1972. 
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Figure  67.  Measured  relative  reflectance  (lower  panel)  and  calculated 
absolute  reflectance  spectra  (upper  panel)  in  the  2-20- pm 
region  for  individual  0.5M  aqueous  solutions  containing 
NaNO^,  or  I^SO^,  The  angle  of  incidence  was 

70. 03°+0,23°.  Distilled  water  was  the  standard  reflector; 
its  relative  reflectance  is  1,0  at  all  wavelengths.  The 
radiant  flux  was  linearly  polarised  with  the  electric- 
field  vector  perpendicular  to  the  plane  of  incidence.  The 
uncertainty  in  the  relative  reflectance  was  about  +0,008  <R+> 
before  the  data  were  smoothed, 
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and  for  individual  0.5M  aqueous  solutions  of  NaNO, , 
KpSO^,  or  NH^HpPOjj,  Reflectance  spectra  for  these 
aqueous  solutions  appear  in  Figs.  66  and  67.  First, 
we  very  briefly  review  sillent  characteristics  of  those 
reflectance  spectra  of  the  solutions  that  contain  an 
individual  chemical  compound, 

(1)  Individual  Aqueous  Solutions 
(a)  Liquid  Water 

The  KpO  monomer  has  three  fundamental 
infrared  active  vibrational  modes  which  are 


evident  at  different  spectral  positions : 
v-|  (A^) »  .  .  symmetric  OH  stretch.  ...  3.10  pm, 

VpCA-jJ.  .  .  breathing  motion.  .  .  6.10  um,  and 
v3(A1).  .  .  antisymmetric  OH  stretch.  .  .  2.90  pm. 

Liquid  water  also  possesses  infrared  active 
lattice  vibrations  which  are  known  as  libraticns. 

One  libraticnal  mode  is  due  to  hindered  oscilla¬ 
tions  of  an  entire  HpO  monomer  that  is  instan¬ 
taneously  hydrogen  bended  in  a  tetrahedral  cage 

64  SR/ 

of  its  four  nearest  neighbors. — This 
librational  mode  is  designated  vT  and  Is 

j—I 

spectrally  centered  at  17.2  ym.  Other  libraticnal 
modes  are  centered  at  wavelengths  greater  than  i 

20  ym. 

Parameters  for  the  oscillatory  modes  of 
water  comprise  the  first  three  lines  of  Table 


I 


XV,  In  Figs.  66  and  67  the 
electromagnetic  radiation  with 


&■'  * 


interaction  of 
the  vibrational 


j 


and  librational  modes  of  water  cause  the 


respective  decrease  ar.d  increase  in  absolute 
reflectance  on  the  short  and  long  wavelength 
sides  of  the  central  spectral  position  for  each 
mode.  We  previously  tabulated  the  spectral 
complex  refractive  index  of  water^  at  about 
298°K, 

(b)  Sodium  Chloride  in  Water 

Sodium  chloride  goes  into  solution  as 
monatomic  ions  Wa+  and  Cl“.  There  are  no 
molecular  vibrational  modes  for  the  monatomic 
ions.  As  NaCl  content  af  an  aqueous  solution 
is  increased,  however,  the  central  positions  of 
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(c)  Potassium  Sulfate  in  Water 

Potassium  sulfate  goes  into  aqueous  solu¬ 
tion  as  monatomic  K+  and  polyatomic  SOjTh  The 
sulfate  ion  has  infrared  active  v^CPj)  and 
vjjCFp)  modes  centerd  at  9.C5  and  6.90  um, 
respectively.  We  note  from  Table  XV  that  the 
infrared  strength  of  the  v-,(F2)  mode  of  SC^ 
is  about  11.76  times  greater  than  that  of  the 
v,.(F„)  mode,  In  the  relative  reflectance  soectrum 

a  ^ 

shorn  in  Fig,  67  for  0.5M  K2S0^  in  water,  the 
predominant  feature  is  the  effect  of  the  strong 
v^(F2)  band  in  the  8-12  ym  region.  The  much 
weaker  v^(F2)  band  is  barely  evident  in  the 
14,9-13  um  region,  In  the  absolute  reflectance 
spectrum  the  v^(F2)  mode  of  SO^.  causes  the 
respective  increase  and  decrease  in  reflectance 
on  the  short  and  long  wavelength  sides  of  this 
erode ' s  central  position.  One  mole  of  K2S0^ 
contains  78. 2g  of  K  and  96. Ig  of  SC^, 

(d)  .Ammonium  Phosphate  (Monobasic)  in  Water 

NHjjHgPOjj  goes  into  aqueous  solution  as 
polyatomic  ions  NKii+  and  PL,PC^,  The  anmonium 
ion  has  infrared  active  v,(F2)  and  v^(F-,)  modes 
centered  at  3.15  and  6.90  um,  respectively. 

Again  the  v,(F2)  mode  has  greater  infrared  band 
strength  than  the  v^(F2)  mode  (see  column  seven, 

-  Table  XV).  The  effect  of  the  v^(F2)  mode  of  HN^ ' 

*  is  difficult  to  see  in  spectra  of  aqueous 

I 

A '  1 

J  ' 
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solutions  because  it  is  an  NH-oscillaticn  occurr¬ 
ing  at  nearly  the  same  spectral  position  as  the 
combined  v,\>,  modes  of  H^O,  The  effect  of  the 
v4(F2)  mode  of  NHj,  appears  at  about  6.9C  urn  as 
a  small  dip  and  then  increases  in  both  the  re¬ 
lative  and  absolute  reflectance  spectra  shown 
in  Fig,  67  for  a  0.5M  solution  of  NH^PO^, 

KgPCTJ  possesses  four  infrared  active 
vibrational  modes  of  moderate  band  strength. 

They  are  listed  in  Table  XV  as  v-jCA^),  v^CA-,), 
v,(B),  v^(B),  The  effects  of  these  vibrational 
modes  are  easily  identified  when  their  spectral 
position  listed  in  Table  XV  are  compared  with  the 
absolute  and  relative  reflectance  spectra  in 
Fig.  67.  One  mole  of  NH^H^PC^  contains  l8g  of 
NH4  and  97g  of  H2POu, 

(e)  Sodium  titrate  in  Water 

Sodium  nitrate  dissolves  in  water  as  Na+ 
and  N0~,  The  nitrate  ion  possesses  a  v^(E) 
mode  centered  at  7.4  ym,  The  effect  of  the 
v^CE)  mode  appears  as  a  respective  decrease  and 
then  increase  in  absolute  and  relative  reflectance 
on  the  short  and  long  wavelength  sides  of  the 
central  position  at  7.4  ym.  This  effect  is 
easily  seen  in  the  spectra  shown  in  Fig.  67  fc’-’ 
Q.5M  NaNO,.  One  mole  of  NaNC,  contains  22,99?  of 

7  Na  and  62g  of  NO,. 

1  J 

1 

J  I 
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(2)  Natural  Saline  Waters 

The  foregoing  discussion  of  the  reflectance 
spectra  of  aqueous  solutions  containing  a  single 
chemical  compound,  together  with  the  molecular 
parameters  in  Table  XV  ,  and  the  chemical  analyses 
shewn  in  Table  XVI  for  the  natural  waters,  make  the 
analysis  of  the  reflectance  spectra  for  the  natural 
waters  a  relatively  simple  task. 

(a)  San  Francisco  Bay,  Pacific  Ocean,  Atlantic  Ocean 
The  water  samples  from  the  San  Francisco 
Bay  and  the  Pacific  Ocean  contained  monatomic 

ions  Ca+2,  Mg+2,  Na+,  K+,  and  Cl"  and  polyatomic 

-  -2 

ions  HCO^  and  SO^  ,  The  concentration  of 
monatomic  ions  was  equivalent  to  about  0.45M 
for  the  San  Francisco  Bay  sample,  and  about 
0.53M  for  the  Pacific-Ocean  sample.  The  con¬ 
centration  of  polyatomic  ions  was  about  0.002M 
for  HCO^~  and  about  0.025M  for  SO^2.  We  believe, 
although  it  was  not  measured,  that  the  Atlantic- 
Ocean  sample  had  a  similar  chemical  composition. 

The  relative  reflectance  spectra  are  shown 
in  Fig.  6l  for  these  three  water  samples.  The 
spectra  exhibit  ail  the  spectral  characteristics 
or  relative  reflectance  spectra  for  aqueous 
solutions  containing  monatomic  ions,  e.g.  the 
spectra  for  1,  3,  and  5M  NaCl  solutions  shown 
in  Fig.  66,  At  2.75  urn  is  a  prominent  peak 
due  to  shifting  of  the  combined  v,v,  modes  of 
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HpO  when  the  monatomic  ions  are  present  in  the 
sample.  In  the  8-11-ym  intervals  the  relative 
reflectance  of  these  three  samples  increases, 
then  decreases  in  the  ll-12ym  interval,  and  then 
begins  to  generally  increase  again  beyond  12  ym; 
this  is  caused  in  part  by  the  monatomic  ions 
shifting  the  vT  modes  of  liquid  water  to 

greater  wavelengths.  The  effect  of  the  v,(F0) 

_2 

mode  of  SO^  at  9.05  ym  was  also  partially 
responsible  for  the  increased  reflectance  in  the 
8-11- interval  of  the  relative  reflectance 
spectra  for  these  three  samples.  The  effect  of 
the  HCO,-  ion  is  not  discemable  in  these 
spectra. 

(b)  Great  Salt  lake  (Utah) 

The  water  sample  from  the  Great  Salt  Lake 

4*2  4*2  41 

contained  monatomic  ions  Ca  ,  Mg  ,  Na  ,  K  , 

_  _  _? 

and  Cl  and  polyatomic  ions  HCO^  and  SO^  . 

The  total  concentration  of  monatomic  positive 
and  negative  ions  was  about  2.65M,  The  concent ra- 
tion  of  KC07  was  0.008M;  and  for  SO^  was 
0.125M. 

The  relative  reflectance  spectrum  for  this 
sample  is  shewn  in  Fig.  62,  The  spectrum  con¬ 
tains  the  prominent  peak  at  2,75  pm  due  a  shifting 
to  shorter  wavelengths  of  the  combined  vn v,  modes 
of  H~0,  In  the  8-11-ym  interval  the  relative 

<L 

reflectance  increases,  then  decreases  in  the 
11-12. 5- ym  interval,  and  then  generally  increases 
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again  beyond  12, 6yn;  this  is  caused  in  part  by 
the  Vr  mode  of  liquid  water  being  shifted  to 
greater  wavelengths  in  this  sample.  The  effect 
of  the  v,(F2)  mode  of  £0^  centered  at  9.05  um 
is  particularly  evident  in  the  8-11-um  interval 
when  the  relative  reflectance  spectrum  for  this 
sample  is  compared  to  that  for  shown  in 
Fig.  67.  The  effect  of  the  v2  mode  of  HCO^- 
is  just  beginning  to  appear  at  6,10  ym  in  the 
spectrum  for  this  sample.  Also,  the  combined 
effect  of  v^vjj  at  7.3  um  for  H  CO^-  is  barely 
discerable  in  the  spectrum, 

(c)  Dead  Sea  (Isreal) 

Chemical  composition  for  the  sample  from  the 

Dead  Sea  was  obtained  from  the  publication  by 
5V 

Neev  and  Emery, because  spectrographic  techniques, 
exclusively,  were  utilised  at  the  USGS 
Laboratory  to  analyse  this  sample;  see  Table  XTV, 

The  sample,  therefore,  probably  contained 

+2  +2  +  +  -  - 

monatomic  ions  Ca  ,  Mg  ,  Na  ,  K  ,  Cl  ,  Br  , 

plus  relatively  small  amounts  of  metals.  The 
concentration  of  monatomic  ions  in  such  a  water 
sample  is  equivalent  to  about  5.60M  each  for 
total  negative  ions  and  for  total  positive  ions, 

HCO~  content  was  probably  about  0.003M,  and 

_p 

SOh  content  about  0.006M, 
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The  relative  reflectance  spectrum  for  this 
sample  is  shown  in  Fig.  62.  The  spectrum 
possesses  only  these  features  characteristic  of 
a  solution  containing  monatomic  ions.  There  is 
a  great  similarity  in  the  spectrum  for  this  sample 
and  the  relative  reflectance  spectrum  shewn  in 

Fig.  66  for  a  5M  NaCl  solution.  The  effects  of 

-  -2 

vibrational  modes  for  HOC,  and  SO^  are  not 
discemable  in  the  spectrum  for  the  sample  from 
the  Dead  Sea. 

(d)  Phosphate  Effluent  Pit  (Central  Florida) 

This  water  sample  contained  monatomic  ions 

Ca+“ ,  Fg+2 ,  Na+,  K+,  Cl-,  and  F-,  The  relative 

reflectance  spectrum  shown  in  Fig.  3  Per  this 

sample  has  the  sharp  peak  at  2,75  um  which  is 

due  to  shifting  of  the  combined  modes  of 

1^0.  The  presence  of  the  smaller  F“  ion  causes 

the  combined  modes  of  f^O  to  shift  to 
56/ 

longer  wavelengths , —  ’whereas  Cl  causes  a 
shift  to  shorter  wavelengths. 

This  water  sample  also  contained  polyatomic 

_2 

ions  SO,,  in  relatively  small  concentrations 
and  contained  PC^J  and  I^PO^  in  much  larger 
concentrations,  The  PCiV  ion  has  infrared  active 
modes  at  about  10  urn  and  ^(F?)  at  19.^  un. 

The  v  (F2 )  mode  is  stronger  in  infrared  activity 
than  the  vli(F2)  mode.  H^PO-  has  the  four  infrared 
active  modes  listed  in  Table  XV  ,  ana  has  additional 


infrared  active  modes  v^(A^) ,  \^(B.)»  and 

v^(Bp)  which  were  not  observed  in  the  spectrum 

of  NH^H^FCjj  shown  in  Big,  67.  The  prominent  peak 

in  the  8-12  urn  interval  of  the  reflectance 

spectrum  shown  in  Big.  63  is  caused  by  the 

v,(B2)  mode  of  and  the  v-^CA^),  v,(A-,), 

v^(B^),  v.,(B2)  modes  of  H2PO^, 

In  the  relative  reflectance  spectrum,  as 

shown  in  Fig.  63  for  this  sample,  there  is  a 

small  peak  in  the  12-1^-um  interval,  The  small 

peak  is  tentatively  attribtued  to  a  vibrational 

mode  of  Si02,  which  also  appears  as  a  strong  band 

in  the  absorption  spectra  of  natural  silicates 

67/ 

such  as  ci-quartz,  Chert,  and  Opal.— 
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B.  Investigations  of  Infrared  Transmittance  of  Powdered  Materials, 

1.  Objectives 

As  previously  stated  in  Section  V.A.b, ,  the  complex 
refractive  index  N(v)  of  natural  minerals  is  a  current 
interest  to  some  scientists.  To  obtain  adequate  specular 
reflectance  measurements,  however,  requires  a  crystal  sample 
of  approximate  minimum  dimensions  of  1  cm.  Many  natural 
minerals,  particularly  the  clays,  do  not  occur  as  large 
homogeneous  crystals  that  can  be  cut  and  polished.  Samples  of 
such  materials  must,  therefore,  be  investigate  in  a  powdered 
form. 


The  objectives  of  this  investigation  are:  (1)  to 
measure  either  a  reflectance  or  transmittance  spectrum  of 
the  powdered  material,  and  (2)  to  obtain  spectral  values  of 
N(v)  from  the  measured  spectrum. 
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2.  Background  Information  for  Discussion  of  Powdered  Materials 
The  acquisition  and  analysis  of  optical  reflectance  or 
transmittance  spectra  of  polished  single  crystals  is  relatively 
simple  when  compared  to  acquisition  and  analysis  of  optical 
spectra  for  powders  which  are  composed  of  small  amorphous 
or  crystalline  particles.  Although  group-theory  analysis 
remains  a  valid  means  for  predicting  the  optical  activity 
of  molecular  vibrations  or  phonons  in  small  amorphous  or 
crystalline  particles,  the  acquisition  of  appropriate  optical 
spectra  and  its  subsequent  analysis  is  complicated  by  diffraction, 
interference  and  total  internal  reflectance  effects  which  are 
characteristic  of  the  particle-size  distribution,  particle 
shape,  random  or  preferred  orientation  of  ’uniaxial  or  biaxial 
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crystals,  and  the  relative  spatial  distribution  of  the 
particles  within  the  powdered  material.  The  Kramers-Kronig 
analysis  is  a  valid  means  for  analyzing  specular  reflectance 
spectra  and  fractional  transmittance  spectra,  but  is  not  validly 
applied  to  diffuse  reflectance  spectra  which  is  typical 
of  powders.  Dispersion-theory  analysis  is  a  "/alia  means  for 
determining  N(v)  of  the  powders  provided  there  exists  an 
adea.uate  theory  for  computing  the  measured  optical  properties 
of  the  powders.  A  wholly  statisfactcry  method  of  acquisition 
and  analysis  of  optical  spectra  must  account  for  all  of  these 
things.  Currently  no  such  wholly  satisfactory  method  exists, 
but  several  methods  are  now  practiced  with  varying  degrees 
of  success  in  research  laboratories. 

cQ/ 

One  method,  developed  by  Ebersol—  and  his  colleagues, 
requires  simultaneous  measurements  of  the  particle-size  dis¬ 
tribution  and  of  scattering  and  attenuation  of  laser  radiation 
by  aribcme  particulates  which  continuously  flow  through  a 
test  chamber.  The  particle-size  distribution  is  measured  by 
use  of  a  Khcller.berg  counter.  Mie-scattering  computations, 
assuming  spherical-particles  are  present,  are  then  empirically 
adjusted  until  a  best  fit  to  the  observed  optical  data  is 
obtained;  thereby  providing  a  value  for  N(v)  at  the  single 
frequency  v  of  the  laser  radiation,  Presently  the  principal 
objections  to  this  method  are  (1)  its  reliance  on  the 
Kncllenberg  counter,  (2)  the  spherical-particle  assumption, 
and  (3)  its  limited  spectral  extent  in  the  infrared. 


5°/ 

A  second  method-*-  is  application  of  Fourier  transform 

spectroscopy  in  the  300-1^00  cm~'L  spectral  range  to  measure 

diffuse  emittance  or  reflectance  spectra  of  loose  powders 

whose  chemical  composition,  crystalline  structure,  and 

spectral  values  of  N(v)  are  known.  The  spectral  values  of 

N(v)  are  known  from  previous  investigations  of  the  optical 

properties  of  large  single  crystals.  The  measured  spectra 

go  / 

are  then  compared  with  spectra  computed  from  a  theory — • 
for  the  spectral  reflectance  and  emittance  of  particulate 
materials.  The  theory,  which  is  initiated  with  the  known 
values  of  M(v)  and  known  particle  sizes  and  shapes,  yields 
a  favorable  comparison  with  measured  spectra.  In  principal 
this  theory  could  be  mathematically  inverted  so  that  its 
application  to  the  measured  spectra  would  yield  N(v),  In 
practice,  however,  the  inversion  would  be  extremely 
difficult  and  one  would  probably  resort  to  empirical  adjust¬ 
ments  of  N(v)  until  a  best  fit  to  the  measured  spectra  is 
obtained  from  the  theory. 

A  third  method  is  to  press  the  particulate  material  in 
a  transparent  matrix  and  then  measure  the  fractional  trans¬ 
mittance  of  the  matrix  plus  the  particulate  material  relative 
to  the  transmittance  of  the  matrix  alone.  In  the  Infrared  KBr 
and  Csl,  and  in  the  far  Infrared  polyethylene,  are  the  pre¬ 
ferred  substances  for  the  transparent  matrix.  This  method 
was  used  by  Volz— 2— ^  to  obtain  H(v)  in  the  250-^  ,C00-cnf~ 
spectral  range  for  Sahara  dust,  volcanic  pumice,  ammonium 
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sulfate,  flyash,  arid  other  aerosols  such  as  soot  and  sea 
salt.  The  principal  objections  to  this  method  are  (1)  the 
particulate  material  is  embedded  in  a  matrix  whose  refractive 
index  thereby  alters  the  relative  refractive  index  of  the 
particulate  material,  (2)  the  attenuation  coefficient  obtained 
from  the  fractional  transmittance  is  the  sum  of  the  absorption 
and  scattering  coefficients,  and  (3)  the  scattering  coefficient 
is  altered  due  to  the  presence  of  the  matrix.  However,  this 

method  does  yield  quantitative  data  over  a  wide  spectral 

♦ 

range.  This  was  the  method  whereby  powdered  materials 
were  investigated  during  the  last  months  of  this  grant. 

A  fourth  method  is  a  variation  of  the  third  method  which 
was  outlined  in  the  proceeding  paragraph.  Father  than 
pressing  the  powder  within  a  matrix  as  in  method  three,  one 
spreads  the  powder  loosely  (about  5  per  cent  packing  fraction) 
on  a  substrate.  Electrostatic  attraction  binds  the  powder 
to  the  surface  of  the  substrate,  which  is  typically  KRS-5, 
crystalline  quartz,  or  soft  polyethylene.  The  transmittance 
spectrum  of  the  substrate  with  the  particulate  material  is 
measured  relative  to  the  transmittance  spectrum  of  the  sub¬ 
strate  alone.  This  avoids  objections  (1)  and  (3)  as  mentioned 
for  the  third  method,  but  objection  (2)  still  holds  for 
short  wavelength  electromagnetic  waves.  At  longer  wave¬ 
lengths,  i.e.  wavelengths  much  greater  than  the  diameters 

go/ 

of  the  individual  particles,  Maxwell-Gamett  theory—'  allows 
one  to  compute  the  measured  transmittance  spectrum  of  the 
powder  if  N(v)  is  known  for  the  bulk  material.  In  those 
cases  where  N(v)  is  net  known  for  the  bulk,  i.e.  in  the 


proposed  investigations,  one  empirically  adjusts  N(v)  'until 

a  best  fit  to  the  measured  spectrum  is  obtained. 

A  fifth  method  is  to  compact  the  particulate  material 

under  pressure  thereby  producing  a  solid  sample  that  can  be 

polished,  and  then  measure  the  specular  reflectance  spectrum, 

66/ 

Schatz —  utilized  this  method  to  measure  the  total 
reflectance  of  unpolished  compacted  powders  composed  of 
transparent  or  opaque  particles  as  a  function  of  the  compacting 

o  Q  o 

pressure.  Using  pressures  from  6.9  x  10  M/irr  to  4.8  x  10  N/m 
3  4 

(10  to  7  x  10  psi)  he  found  that  powders  composed  of 
transparent  particles  underwent  a  large  decrease  in  reflectance 
as  the  pressure  was  increased,  and  powders  composed  of 
opaque  particles  increased  in  reflectance  as  the  pressure  was 
increased.  As  a  powder  is  compacted  its  optical  properties . 
approach  those  of  a  large  bulk  sample  of  the  amorphous  or 
poly crystalline  material.  This  method  has  some  merit  in 
investigations  of  optical  properties  provided  (1)  some 
practical  compacting  pressure  is  reached  above  which  the 
reflectance  does  not  vary,  (2)  the  compacted  powder  can  be 
pressed  with  a  polished  die  or  can  be  polished  so  that 
specular  reflectance  can  be  measured,  (3)  the  crystalline 
structure  remains  unchanged  for  the  individual  small  crystals, 
and  (4)  the  small  crystals  are  known  to  be  either  randomly 
or  preferrentially  oriented.  The  Xramers-Krcnig  analysis 
-  would  be  valid  when  applied  to  the  specular  reflectance 

*  spectrum  of  the  polished  compacted  powder.  This  is  the 

|  method  whereby  we  are  investigating  powdered  materials  during 

the  follow-on  contract. 
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3.  T(v)  Spectra  for  Powdered  'Materials  Pressed  in  KBr  Pellets 

Materials  investigated  were  these  listed  in  Table  '7/ II. 
Relative  transmittance  spectra  T(v)  =  T(v)  /T(v)m  were 
measured,  where  T(v)  denotes  the  fractional  transmittance 
for  a  KEr  pellet  containing  the  powdered  material  and  T(v)m 
denotes  the  fractional  transmittance  for  a  KEr  pellet  without 
powdered  material.  The  T(v)  spectra  are  presented  in  graphical 
fora  for  the  various  materials  in  the  Figures  listed  in 
column  4  of  Table  XVII, 

Table  XVII.  Powdered  materials  pressed  in  KBr  pellets  for  which  relative 
transmittance  T(\>)  was  measured. 


Material 

Concentration* 

Spectral  Range 

T(v) 

Mcntmorillonite 

io~2,  io"3,  io_il 

300-5,000  cm"1 

Fig.  63 

Colemanite 

10-2,  10~3,  io-4 

300-5,000 

Fig.  69 

Kemite 

IT2,  10-3,  10“4 

300-5,000 

Fig.  70 

Kaolin 

10-2,  10“3,  10-4 

300-5,000 

Fig.  71 

Lampblack 

2  x  10-5 

300-5,000 

Fig.  72 

Active  Charcoal 

2  x  10" 5 

300-5,000 

Fig.  73 

Wood  Charcoal 

2  x  10"14 

300-5,000 

Fig.  74 

Aluminum 

2  x  10~2 

1,665-5,000 

Fig.  75 

Kaolin 

2  x  10-2 

1,665-5,000 

Fig.  76 

Mcntmorillonite 

2  x  10"2 

1,665-5,000 

Fig.  77 

Colemanite 

2  x  ID-2 

1,665-5,000 

Fig.  73 

Illite 

2  x  10"2 

1,665-5,000 

Fig.  79 

Kemite 

2  x  10~2 

1,665-5,000 

Fig.  80 

*See  following  section  for  a  discussion  of  the  concentrations. 


4.  Discussion  of  T(v)  Spectra  for  Powdered  Materials 

The  KEr  pellets  containing  the  powdered  materials  were 
prepared  as  follows.  For  the  first  four  materials  listed  in 
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region  for  Lampblack  at  a  concentration  of  2  x  10  ;  in  KBr  pellets. 


region  for  activated  characoal  at  a  concentration  of  2 
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‘or  Kemite  at  a  concentration  of  2  x  10  in  KBr  pe 


Table  XVII,  a  stock  mixture  consisting  of  1  part  powder  per 
100  parts  XBr  (6g  cf  XBr)  was  prepared.  This  provided  the 
10~c  concentration.  Serial  dilution  of  portions  cf  the  stock 

_2i 

mixture  then  provided  the  10  '  and  10  concentrations.  A 
similar  serial  dilution  process  provided  the  concentrations 
for  lampblack,  activated  charcoal,  and  wood  charcoal.  The 
remaining  six  materials  in  Table  XVII  were  prepared  by  mixing 
0.12g  of  powdered  material  with  6g  of  XBr. 

Each  pellet  was  prepared  by  use  of  300  mg  of  KBr  or 
powder  +  KBr.  The  pellets  were  pressed  in  a  polished  die  at 
a  pressure  of  87,000  psi.  The  pellets  were  1.30+0.03  mm  in 
diameter  and  0,8t2+0.015  mm  thick.  Just  prior  to  pressing 
the  pellet,  about  5CCmg  cf  the  powder — KBr  mixture  was 
re ground  by  use  of  a  vibrcmill.  The  300mg  required  for  the 
pellet  was  then  obtained  from  this  reground  material.  All 
powdered  materials,  except  the  KBr,  were  obtained  from 
Wards  Scientific. 

Montmorillonlte  is  an  optically  biaxial  material  whose 
chemical  composition  is  Al2Si|+0||(0H)2,hH20.  The  T(v)  spectra 
for  this  material  are  presented  in  Figs.  68  and  77, 

Colemanite  is  an  optically  biaxial  material  whose  chemical 
composition  is  Ca^^C The  T(v)  spectra  for  this 
material  are  presented  in  Figs.  6°  and  78.  The  colemanite 
and  montmorillonlte  spectra  both  show  the  O-H  stretching 
bands  in  the  3,300  cm~  region.  We  did  net  make  a  group 
analysis  of  these  minerals,  and  therefore  did  not  assign  the 
bands  in  the  lower  wave— number  region. 


Kemite  is  an  optically  biaxial  material  whose  chemical 
composition  is  4^0,  The  T(v)  spectra  for  this 

material  are  shewn  in  Figs,  70  and  80,  The  0-H  stretch 
bands  are  also  apparent  in  the  spectra  of  this  material. 

Kaolin  is  an  optically  biaxial  material  whose  chemical 
composition  is  A^S^ClCOH)^.  The  T(v)  spectra  for  this 
material  are  shown  in  Figs.  71  arid  76. 

lampblack  consists  of  opaque  carbon  soot  particles.  The 
T(v)  spectrum  for  this  material  is  shewn  in  Fig.  72. 

Note  the  absence  of  infrared  active  bands  for  this  material. 

Activated  charcoal  also  consists  of  opaque  particles. 

Its  T(v)  spectrum  is  shown  in  Fig.  73.  Again  there  is  an 
absence  of  infrared  active  bands. 

Wood  charcoal  was  more  transparent  than  lampblack  or 
activated  charcoal,  and  there  is  one  weak  infrared  active 
band  near  1,430  cm-'''.  The  T(v)  spectrum  for  this  material 
is  shown  in  Fig.  74. 

Aluminum  powder  consists  of  optically  opaque  particles. 
The  T(v)  spectrum  for  this  material  is  shorn  in  Fig.  75. 

Again  there  is  an  absence  of  Infrared  active  bands, 

Illite  is  an  optically  biaxial  material  whose  variable 
chemical  composition  is  K,_xAl2Si,+xAl1_xO-l0(0H)o.  The  T(v) 
spectrum  for  this  material  is  shown  in  Fig.  75. 

Attempts  to  determine  attenuation  cross  sections  from 
these  T(v )  spectra  yielded  inconsistent  results  for 
reasons  discussed  in  Section  V.B.2. 

The  powdered  materials  whose  T(v )  spectra  have  teen 
presented  here  are  currently  being  investigated  by  use  of 


reflectance  methods  during  the  follow-on  contract. 


The  T(v)  spectra  shown  in  Figures  68-80  constitute 
raw  data  which  have  not  been  corrected  for  random  errors 
picked  up  during  data  transmission  (anomalous  sharp 
spikes),  and  the  use  of  Fractional  transmitting  sector 
wheels  during  data  acquisition  (abrupt  discontinuities 
in  the  spectra). 
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D.  Other  Investigations 


1.  Split-pulse  laser  method  for  measuring  attenuation  coefficients 

of  transparent  liquids:  Application  to  deionized  filter  water  In  the 
visible  region 

Marvin  R.  Querry,  Philip  G.  Cary,  and  Richard  C.  Waring 
Department  of  Physics,  University  of  Missouri-Kansas  City, 

Kansas  City,  Mo,  64110 

a.  Abstract 

Attenuation  coefficients  a  (A)  for  collimated  quasimono chroma tic 
radiation  passing  through  deionized  filtered  water  were  measured  throughout 
the  418.6  to  640.3  nm  wavelength  region  by  use  of  a  split-pulse  laser 
method  \^ich  employs  reference  and  sample  cells  arranged  in  a  geometry 
similar  to  that  of  a  Michelscn  interferometer.  The  radiant  source  was  a 
pulsed  wavelength- tunable  dye  laser  possessing  a  relatively  short 
coherence  time.  This  paper  includes  descriptions  of  theoretical  and 
experimental  techniques  applicable  to  the  split-pulse  laser  method,  and 
includes  a  tabulation  of  a(A)  measured  for  deionized  filtered  water  at 
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b.  Introduction 

From  a  survey  of  the  scientific  literature^  on  the  complex 

refractive  index  N(X)  =  n(^)  +  ik(x)  of  liquid  water,  and  from  a  recent 

2/ 

article  by  Kopelevich,—  we  became  increasingly  aware  of  considerable 

3-12/ 

variations  in  the  values  measured  by  various  investigators -  for 

the  attenuation  coefficient  a(x)  =  4irk(A) /A  of  liquid  water  in  the  visible 

3  A 

region  where  ct(X)  is  the  order  of  from  10~  to  10"  .  Although  such 
variations  in  a  (A)  are  in  part  directly  attributable  to  varying  quantities 
of  particulates  and  dissolved  chemicals  in  water  samples,  different 
experimental  methods  and  procedures  can  also  introduce  variations  in 
measured  values  of  ct(A).  Some  problems  that  arise  in  conventional 
absorption  spectroscopy  are  the  following:  (1)  incorrect  determination 
of  reflectance  of  cell  windows  by  measuring  reflectance  losses  for  an 
empty  cell  and  applying  those  measurements  to  a  filled  cell,  (2)  insufficient 
collimation  of  radiant  flux  from  an  extended  source  thus  precluding  the 
use  of  adequate  length  cells  to  provide  attenuation  of  from  20  to  80%  of 
the  incident  radiant  flux,  (3)  insufficient  collimation  thus  allowing  radiant 
flux  to  be  reflected  from  the  side  walls  of  long  cells,  (4)  temporal 
variations  in  the  brightness  of  the  radiant  source  during  single-beam 
long-path  investigations  where  transmittance  must  be  measured  at  different 
times  through  two  or  more  cells  of  different  length,  and  (5)  additional 
signal  from  a  photomultiplier  detector  due  to  stray  light  reaching  the 
detector.  Although  these  problems,  if  anticipated,  can  be  overcome  in 
conventional  spectroscopy  by  proper  experimental  techniques  they  do 
constitute  points  of  serious  concern  in  the  measurement  of  the  attenuation 
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coefficients  of  transparent  liquids. 

The  principal  objective  of  this  investigation  was  to  develop  an 
experimental  method  that  overcame  the  five  problems  listed  above.  The 
problems  were  overcome  in  the  following  manner.  A  wavelength  tunable  dye 
laser  was  chosen  as  the  radiant  source  thus  assuring  proper  collimation  of 
radiant  flux  and  elimination  of  stray  light.  Reflectance  at  the  end  windows 
of  the  absorption  cells  was  accounted  for  by  choosing  a  split-beam 
experimental  apparatus  that  was  similar  in  geometry  to  a  Michelson  inter¬ 
ferometer.  Temporal  variations  in  the  brightness  of  the  radiant  source  were 
considered  by  choosing  a  pulsed  dye  laser,  rather  than  a  CTv  dye  laser,  as 
the  radiant  source  and  applying  repetitive  single-pulse  ratioing  techniques 
by  use  of  a  dual  channel  boxcar  integrator.  Also,  a  pulsed  dye  laser  with 
a  relatively  small  coherence  length  was  chosen  thereby  avoiding  interference 
fringes  in  optical  components  of  the  experimental  apparatus.  The  split-pulse 
laser  method  that  we  developed,  based  on  the  preceeding  choices,  for 
measuring  a(X)  of  highly  transparent  liquids  is  presented  in  the  following 
sections.  The  method  is  presented  by  describing  its  application  to  measuring 
a(X)  for  deionized  filtered  water. 
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Figure  8l 

A  block  diagram  of  the  instrumentation  used  to 
measure  the  attenuation  coefficient  a(A)  of  deionized 
filter  water  in  the  4i8,6-  to  640.3-nm  region.  Further 
details  are  provided  in  Section  II. 
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c .  Experimental  Methods 

The  radiant  source  for  this  experiment  was  a  wavelength  tunable  dye 
13/ 

laser —  which  was  optically  pumped  by  an  AVCO  C950  pulsed  nitrogen  laser. 
Pulses  from  the  dye  laser  were  of  about  10  nsec  duration,  were  from  1  rn 
to  2  rn  in  spectral  width,  and  had  with  a  repetition  rate  of  50  pps. 

Radiant  flux  from  the  dye  laser  was  passed  through  a  limiting  aperature 
and  was  then  further  collimated  by  a  converging  lens  in  order  to  retain 
the  radiant  flux  well  within  the  sample  and  reference  cells  described 
hereinafter.  After  passage  through  the  convergent  lens  the  beam  of  radiant 
flux  encountered  the  beam  splitter  depicted  in  Fig.  8l  where  it  was  split 
into  two  beams.  One  beam,  that  reflected  by  the  splitter,  then  passed 
through  a  reference  cell  of  length  Lr  =  24,3  cm,  and  the  other  beam,  that 
transmitted  by  the  splitter,  passed  through  a  sample  cell  of  length 
Ls  =  499.1  cm.  The  reference  and  samples  cells  were  both  filled  with 
purified  deionized  water  which  was  previously  passed  through  a  0.22  ym 
Millipore  filter  to  remove  macroscopic  particulates.  The  beam  passing 
through  the  reference  cell,  and  the  beam  passing  through  the  sample  cell, 
then  encounter  plane  mirrors  and  M2,  respectively,  where  they  were 
reflected  back  through  the  respective  cells  to  again  encounter  the  beam 
splitter.  After  the  second  encounter  with  the  beam  splitter  that 
portion  of  the  reference-cell  beam  transmitted  by  the  splitter,  and  that 
portion  of  the  sample-cell  beam  reflected  by  the  splitter,  were  incident 
on  a  single  ITT  FV114A  photodiode  detector.  The  cells  were  slightly 
misaligned  relative  to  the  direction  of  propagation  of  radiant  flux  so  that 
radiant  flux  reflected  from  both  surfaces  of  the  cell  windows  was  not 


incident  on  the  detector. 


184 


Due  to  the  difference  in  length  of  the  two  cells,  the  pulse  of  radiant 
flux  passing  through  the  sample  cell  lagged  the  pulse  passing  through  the 
reference  cell  by  about  40  nsec  at  the  detector.  The  signal  from  the 
detector  was  applied  to  a  Princeton  Applied  Research  (PAR)  model  162  dual 
channel  boxcar  integrator  equipped  with  two  PAR-164  integrator  modules. 

Gate  A  of  the  integrator  was  adjusted  to  sense  the  signal  associated  with 
the  pulse  passing  through  the  sample  cell,  and  gate  B  was  adjusted  to 
sense  the  signal  associated  with  the  pulse  passing  through  the  reference 
cell.  The  gate  widths  were  about  10  nsec.  The  integrator  provided  an 
output  signal  A(X)/B(X)  to  the  input  of  a  PAR-260  analog/digital  converter, 
where  A(X)  and  B(X)  are  proportional  to  the  radiant  intensity.  The  output 
from  the  PAR-260  was  then  applied  to  a  Hewlett-Packard  9820  programmable 
calculator  where  data  were  printed  on  paper  tape  and  were  also  stored  on 
cassette  magnetic  tapes.  The  calculator  was  programmed  to  compute  the 
attenuation  coefficient 


a  a) 


2~(l' -l  r  • 


(1) 


As  explained  in  the  following  Section,  Eq.  (1)  is  deemed  to  be  valid  for 
analysis  of  data  from  this  experiment. 

Each  cell  was  constructed  from  cylindrical  pyrex  tubing  of  4.1  cm 
nominal  inside  diameter.  End  windows  for  all  cells  were  cut  from  the  same 
plate  of  pyrex  glass  so  they  would  have  the  same  transmission  and  reflection 
properties.  The  end  windows  were  fused  to  the  pyrex  tubing  by  the  glassblowe 
thus  avoiding  the  use  of  0-rings  or  any  material  other  than  glass.  An 
inlet/outlet  tube,  also  of  pyrex  glass,  was  fused  to  the  side  wall  near 
each  end  of  each  ceil.  The  cells  were  initially  cleaned  with  chromic  acid. 
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then  with  Alconox,  and  were  then  flushed  many  times  with  clean  deionized 
filtered  water.  The  entire  cleaning  procedure  was  applied  twice  to  each 
cell.  Care  was  taken  to  keep  the  outer  surfaces  of  the  windows  clean  and 
dust  free. 

Clean  deionized  filtered  water  was  sealed  inside  the  cells  for  about 
180  days  prior  to  serious  attenuation  measurements,  i:hus  allowing 
microscopic  air  and  vapor  bubbles  to  dissipate.  The  bubbles  were  initially 
observed  in  the  path  of  the  laser  beam  as  profuse  scattering  centers. 

Sullivan,—^  who  used  a  60-cm  long  cell  made  from  acrylic  plastic  and 
equipped  with  pyrex  windows  sealed  with  o-rings,  noted  that  in  24  hours 
the  fractional  transmittance  of  freshly  distilled  water  decreased  from 
88.5%+lZ  to  80*  or  less  at  a  wavelength  of  470  nm.  Because  Sullivan 
attributed  the  decrease  in  fractional  transmittance  to  the  deterioration 
of  the  optical  purity  of  distilled  water,  we  were  concerned  about  possible 
leaching  of  the  pyrex  cells  by  deionized  water.  However,  no  systematic 
increase  or  decrease  in  the  fractional  transmittance  of  the  deionized  water 
was  observed  throughout  the  418.6-  to  640, 3-nm  region  during  the  four 
months  when  we  were  making  several  measurements  of  the  spectral  attenuation. 
If  leaching  did  occur,  then  it  reached  equilibrium  during  the  initial 
180  days  that  the  deionized  water  was  sealed  in  the  pyrex  cells,  or  the 
leached  borosilicates  did  not  affect  the  absorption  measurements. 

Three  types  of  beam  splitters  were  tried  before  a  satisfactory  one  was 
obtained.  The  first  was  a  thin  pellicle  with  a  partially  reflecting  coating 
on  one  surface.  The  pellicle  was  unsatisfactory  because  its  thickness 
was  less  than  the  coherence  length  of  the  laser  pulse,  and  interference 
fringes  were  observed  as  the  laser  was  scanned  through  the  region  of 
wavelength  tunability.  The  second  was  a  mirror  type  (Edmund  Sci.  if 578)  beam 
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splitter  with  a  partially  reflecting  coating  on  one  side  of  a  5-nan 
thick  glass  flat.  The  second  beam  splitter  was  unsatisfactory  because 
the  sample-cell  and  reference-cexl  laser  beams  were  not  subject  to  exactly 
the  same  net  optical  treatment  after  interacting  twice  each  with  the 
beam  splitter.  A  6 -mm  thick  glass  optical  flat  proved  to  be  a  very 
satisfactory  beam  splitter.  Additionally,  mirrors  and  M,,  were  inter¬ 
changed  several  times  and  the  measured  attenuation  remained  the  same  each 
time. 

The  temperature  of  the  water  was  maintained  at  26.4+l,7°C. 
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d.  Theoretical  Background 

The  macroscopic  form  of  Maxwell's  equations  and  the  corresponding 
constitutive  relations,  in  cgs  gaussian  units,  applied  to  an  infinite, 
homogeneous,  linear,  isotropic,  magnetic,  conducting  medium  which  is  free 
of  net  electric  charge  provides  the  wave  equations 


(v2  Ampc  3_ 
K  ~  2  3t 


ye  3 _ ,,, 

~~2  ~2^ 

c  at 


(2) 


where  is  any  one  of  the  components  of  the  electric  field  intensity  E  or 
che  magnetic  field  intensity  H,  c  is  the  speed  of  electromagnetic  waves  in 
vacuum,  y  is  the  relative  magnetic  permeability,  e  is  the  dielectric 
function  which  originates  from  the  physical  behavior  of  bound  charge, 
o  is  the  electrical  conductivity  which  is  attributed  to  free  charge  carriers 
responding  to  the  incident  electromagnetic  wave,  and  u  =  2mc/X  is  the 
angular  frequency.  Equation  (2)  has  a  plane  monochromatic  wave  solution 
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where  subscripts  r  and  i  denote  real  and  imaginary  parts,  respectively.  The 
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complex  refractive  index  N  =  n  +  ik  of  the  material  medium  is  defined  as 
the  ratio  of  c,  the  speed  of  electromagnetic  waves  in  vacuum  where  c  =  1, 
y  =  1,  and  o=0  for  all  u,  to  the  phase  velocity  =  w/k  of  the  electro¬ 
magnetic  wave  in  the  material  medium: 

n  +  ik  =  c(<r  +  iK^)/u  ,  (9) 

For  simplicity  we  have  dropped  the  functional  notation.  The  absorption  of 

electromagnetic  flux  moving  in  the  +z-direction  from  zq  to  z  is 

-a  (z-z  ) 

I  =  I0e  a  °  >  (10) 


where  =  2k^,  and  I  and  IQ  are  radiant  intensities  at  z  and  z0,  respectively. 
The  parameter  a  in  Eq.  (10)  represents  a  true  absorption  coefficient 

a 

corresponding  to  the  electromagnetic  properties  y,  e,  and  o  of  the  material 

responding  to  the  incident  electromagnetic  wave.  In  the  visible  region  of 

the  electromagnetic  spectrum  we  assume  for  liquid  water  that  y  =  1,  a  =  0, 

2 

<<  er,  and  =  n  so  that;  by  use  of  Eqs.  (5),  (7),  and  (8); 


k2  =  (u2  +  i2bw)/a, 

(ID 

where 

=  c2/er,  and 

(12) 

b  =  (mei)/(2er)  . 

(13) 

In  practical  situations,  especially  in  investigations  of  liquid  water, 
the  incident  electromagnetic  radiation  is  absorbed  and  is  also  scattered. 
The  measured  attenuation  of  electromagnetic  radiation  by  water,  therefore, 
is  described  by 


I  =  I0e 


- (a  +  a  ) (z-z  ) 
a  s  o 


% 


where  is  the  absorption  coefficient  of  Eq.  (10),  and  as  is  the  attenuation 
coefficient  due  to  scattering.  The  scattering  coefficient  ag  is  comprised 
of  the  sum  of  two  contributions: 
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a  =  a  +  a 
s  ps  ms 


(15) 


where  ct  ^  is  an  attenuation  coefficient  due  to  scattering  by  microscopic 
particles,  and  a  is  an  attenuation  coefficient  due  to  molecular  scattering, 
i.e.,  Rayleigh  scattering  due  to  local  thermal  fluctuations  in  the  molecular 
density. 

The  attenuation  coefficient  a(w)  for  electromagnetic  radiation  of 

ms 

angular  frequency  u  may  be  calculated  by  use  of  the  Einstein-Smoluchovski 
equation 


<x(w) 


ms 


k"Tg  .u. 

6tt  V 


[e(u)-l]  [e (oi)+2 ] 


(16) 


where  k"  is  the  Boltzmann  constant,  T  is  the  absolute  temperature,  and 
-11  2 

3  =  4.508  x  10  cm  /dyne  is  the  isothermal  compressibility  of  water  at 

2 

26.4°C.  In  the  visible  region  where  e(u)  =  n(w)  for  water  we  have 

2 


a(a°ms  =  ^  (^)4  {[n(m)2-l][n(m)2+2]} 


(17) 


Values  of  n(“)  for  water  are  readily  interpolated  from  those  given  in  Ref. 

(1). 

In  the  present  investigation  we  measure,  by  use  of  Eq.  (14)  the 
total  attenuation  coefficient  a  =  a  +  a  of  liquid  water  for  electromagnetic 
radiation  from  a  pulsed  dye  laser.  Pulses  from  the  dye  laser  are  about  10  nsec 
duration  and  of  spectral  width  Aw  =  8  x  1022  sec  2  at  w  =  2.73  x  TO2"*  sec  2 
(640nm) ,  and  Am  =  1.7  x  1022  sec  ^  at  u  =  4.5  x  1022  sec  2  (418nm) .  The 
coherence  time  for  the  laser  radiation  is  tQ  =  2ir/Aw.  As  such  pulses  of 
electromagnetic  radiation  propagate  through  the  water  they  are  attenuated 
and  the  shape  of  the  pulses  may  also  be  distorted.  We  now  examine  the  possible 
distortion.  The  pulses  are  regarded  as  quasimonochromatic  radiation  in  that 
mathematically  they  are  a  Fourier  superposition  of  a  finite  but  continuous 
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14  / 

spectral  distxi  bution  of  plane  electromagnetic  waves.  Stratton —  has 
provided  the  mathematical  rigor  for  examining  the  behavior  of  such  pulses; 
to  wit 


Y(z,t)  =  y  ek°/af (t+z/a)  +  i  e  bZ/^af  (t-z/a) 

t+z/a  _ _ 

+  f  e"bt  /  f(£)ebeJo  [(b/a)/z2-a2(t-B)2  ]dg 
t-z/a 

t+z  /  a 

+  fe“bt  /  F(0)ebe  Jo[(b/a)/z2-a2(t-6)2'  ]dB,(18) 
t-z/a 

where  iKz.t)  is  a  cartesian  component  of  E  or  H,  z  is  the  distance  measured 
from  z=0  along  the  axis  of  propagation,  J0 [ (b/a)  *4^-a^(t-S)^]  is  the  zero- 
order  Bessel  function  and  its  argument,  Jq"1  is  the  partial  derivative  of  J0 
with  respect  to  z,  f(t)  =  K0,t),  F(t)  =  (3^/3z)z=0>  and  b  and  a  are  defined 
by  Eqs.  (13)  and  (12),  respectively.  We  assume  a  simple  model  for  quasi- 
monochromatic  radiation  of  the  form 


£  =  £  gi(<*  r-ut+4>(t) ) 
o 


(19) 


where  $(t)  is  a  random  step  function  of  period  tq  the  coherence  time,  and 
of  range  0  <_  <j>(t)  £  2m.  Equation  (19)  represents  a  sequence  of  finite  waves, 
each  with  wave  vector  tc ,  coherence  time  to  =  2m/Aoj,  but  with  randomly 
differing  phases  <f>(t).  Next,  we  consider  one  of  the  finite  waves  in  the 
sequence  indicated  by  Eq.  (19).  We  represent  a  component  of  E  or  H  by  a 


wave  moving  in  the  +z-direction  as 
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A/  iK2  -iicz.  -iut 
iKz.t)  ■  A(e  +  e  )e 

for  j  1 1  _<  t  / 2 ,  and 

Kz,t)  =  G 

for  jtj  >tq/2.  We  thus  find  that  Eq.  (18)  becomes 

Kz.t)  =  A(e  £bz/a-iu(t+z/ a)  ]+e  t“bz/  a-im°(t-z/  a)  ]  > 
t+z/a 

+  aAe~bt  /  e(b-ia))8J^^(b/a)V/z2_a2(t_e)2'  jde>(22) 

t-z/a 

where  the  pulse  moving  in  the  +z-direction  is  the  one  of  interest  for  our 
application.  The  integral  in  Eq.  (22)  represents  the  distortion  of  the 
pulse  in  time  t  at  position  z.  Ey  use  of  data  from  Ref.  (1)  we  note  that 
2  x  10^  cm/sec  <  a  <  3  x  10^  cm/sec,  (23) 

2b/a  =  aa  <  4  x  10  ^  cm  ^  .  (24) 

In  the  present  investigation  the  pluse  propagates  a  maximum  distance  of  10m 

in  the  water,  but  by  conservatively  using  z  =  20  m,  we  find  the  value  of  the 

integral  in  Eq.  (22)  to  be  on  the  order  of  A  x  10  However,  the  amplitude 

A  of  the  pulse  propagating  20  m  in  water  is  attenuated  to  only  something 

_2 

on  the  order  of  A  x  10  .  Thus,  for  the  present  investigation,  any  distortion 
of  the  pulse  by  the  water  is  at  least  four  orders  of  magnitude  less  than  the 
attenuation,  and  we  can  quite  accurately  represent  the  pulse  moving  in  the 
+z-direction  by 


(20) 

(21) 


A  e~bz/a  eiu (t-z/a) 

♦  C*,t)+  =  I 

0 


1=1  i  V2 

M  >  V2  • 


(25) 
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The  attenuation  of  radiation  described  by  Eq.  (25)  is 


1=1  e 
o 


-2bz/a 


The  absorption  coefficient  for  water,  therefore,  is  a  =  2b/a,  and 
2b/a  =  2k^  provided  p  =  1,  a  =  0,  and  e(w)^  <<  e(ai)r  as  was  previously 
assumed  for  water  in  the  visible  region. 

As  a  final  point,  we  note  that  when  quasimonochromatic  radiation  is 
incident  on  two  boundary  surfaces  between  three  electromagnetically 
dissimilar  media  (e.g.  air/glass/water)  multiple  reflections  at  the 
boundaries  of  the  intermediate  material  can  provide  constructive  and 
destructive  interference  of  the  electromagnetic  radiation.  No  interference 
is  observed,  however,  when  the  coherence  time  tq  is  less  than  the  time 
required  for  the  radiation  to  traverse  two  thicknesses  of  the  intermediate 
material.—^  The  smallest  separation  of  two  such  boundaries  in  the  present 
investigation  was  3.17  mm,  the  thickness  of  the  pyrex  windows  at  the  end 
of  the  cells  containing  the  water.  The  minimum  time  for  radiation  to 
propagate  through  a  single  thickness  of  a  glass  window  was  about  10to. 
Interference  effects,  therefore,  were  not  observed  during  this  experiment. 


L9 4 


TABLE  XVIII.  Wavelength  A,  attenuation  coefficients  a(A),  Rayleigh 
scattering  coefficients  a(A)ms,  the  sum  ot  the  absorption 
coefficient  a (A)  and  the  particle  scattering  coefficient  a (A)  , 

and  wavenumber  A-l  for  deionized  filtered  water.  F 


A 

a(A) 

a(A)s 

a(A)a+ct(A) 

A'1 

nm 

cm-^ 

cm~l 

r 

cm“l 

cm-^ 

418.6 

5.161xl0-4 

5 , 015xl0-5 

4.659xl0-4 

23889.2 

421.9 

4.960 

4.855 

4.475 

23702.3 

425.2 

4.890 

4.701 

4.420 

23518.3 

428.4 

4.944 

4.557 

4,488 

23342.7 

431.7 

5.088 

4,415 

4.647 

23164.2 

434.9 

4.975 

4.232 

4.547 

22993.8 

438.2 

4.937 

4.150 

4.522 

22820.6 

441.5 

5.040 

4.022 

4.638 

22650.1 

444.7 

5.046 

3,904 

4.656 

22487.1 

448.0 

4.811 

3.786 

4.432 

22321.4 

450.4 

4.108 

3.703 

3.738 

22202.5 

453.5 

4.353 

3.599 

3.993 

22050.7 

455.9 

4,473 

3.521 

4.121 

21934.6 

458.3 

4.045 

3.445 

3.700 

21819.8 

460.6 

3.840 

3.374 

3.503 

21710.8 

463.0 

3.753 

3.302 

3.423 

21598.3 

465.4 

3.692 

3.232 

3.369 

21486.9 

468.3 

4.076 

3,150 

3.761 

21353.8 

470.7 

4.061 

3.084 

3.753 

21245.0 

473.0 

4.220 

3.022 

3.913 

21141.6 

475.0 

3.889 

2.969 

3.592 

21052.6 

477.3 

3.997 

2.910 

3.706 

20951.2 

479.6 

4.008 

2.853 

3.723 

20850. 7 

481.8 

3.981 

2.799 

3.  701 

20755.5 

484.2 

3.933 

2.742 

3.659 

20652.6 
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486.6 

3.935 

2.686 

3.666 

20550.8 

489.0 

3.919 

2.632 

3.656 

20449.9 

491.4 

3.933 

2.578 

3.675 

20350.0 

493.8 

3.985 

2.527 

3.732 

20251.1 

495.7 

4.374 

2.487 

4.125 

20173.5 

497.9 

4.618 

2.441 

4,374 

20084.4 

500.3 

4.754 

2.393 

4.515 

19988.0 

502.5 

5.031 

2.349 

4.796 

19900.5 

504.8 

5.046 

2.305 

4.815 

19809.8 

507.2 

5.412 

2.260 

5.186 

19716.1 

509.5 

5.662 

2.218 

5.440 

19627.1 

511.8 

6.084 

2.177 

5.866 

19538.9 

514.2 

6.198 

2.135 

5.985 

19447.7 

516.6 

6.699 

2.093 

6.490 

19357.3 

518.9 

6.631 

2.055 

6.425 

19271.5 

521.2 

6.662 

2.017 

6.460 

19186.5 

523.6 

6.722 

1.979 

6.524 

19098.5 

525.9 

6.408 

1.943 

6.214 

19015.0 

528.2 

6.507 

1.908 

6.316 

18932.2 

530.4 

6.553 

1.875 

6.365 

18853.7 

532.8 

6.650 

1.840 

6.466 

18768.8 

535.0 

6.943 

1.809 

6.762 

18691.6 

536.6 

6.851 

1.787 

6.672 

18635.9 

533.2 

6.942 

1.764 

6.766 

18580,5 

539.7 

7.413 

1.744 

7.239 

18528.3 

541.2 

7.134 

1.724 

6.962 

18477.5 

542.3 

7.278 

1.703 

7.108 

18423.0 

544.3 

7.255 

1.683 

7.087 

18372.2 
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545.8 

7.179 

1.664 

7.013 

18321,7 

547.8 

7.539 

1.639 

7.375 

18254.8 

550.1 

7.783 

1.610 

7.622 

18178.5 

552.4 

7.994 

1.584 

7.836 

18102.8 

554.5 

8.135 

1.560 

7.979 

18034.3 

556.8 

8.226 

1,534 

8.073 

17959.8 

559.0 

8.314 

1.510 

8.163 

17889.1 

561.2 

8.425 

1.487 

8.276 

17819.0 

563.4 

8.572 

1.464 

8.426 

17749.4 

565.6 

8.748 

1.441 

8.604 

17680.3 

567.7 

8.930 

1.420 

8.788 

17614.9 

570.0 

9.092 

1.397 

8.925 

17543.9 

580.3 

1. 085xl0-3 

1.298 

1. 072xl0**3 

17232.5 

582.5 

1.139 

1.278 

1.126 

17167.4 

584.6 

1.222 

1.259 

1.209 

17105.7 

586.8 

1.323 

1.239 

1.311 

17041.6 

588.9 

1.434 

1.^20 

1.422 

16980.8 

591.0 

1.562 

1.202 

1.550 

16920.5 

593.2 

1.698 

1.184 

1,686 

16857.7 

595.4 

1.864 

1.166 

1.852 

16795.4 

597.5 

2.085 

1.148 

2.074 

16736.4 

599.7 

2.372 

1.131 

2.361 

16675.0 

601.9 

2.639 

1.114 

2.628 

16614.1 

604.0 

2.798 

1.099 

•  2, 787 

16556.3 

606.2 

2.865 

1.083 

2.854 

16496.2 

609.1 

2.913 

1,063 

2.902 

16417.7 

611.2 

2.939 

1.048 

2.929 

16361.3 

613.3 

2.959 

1.034 

2.949 

16305.2 

615.5 

2.985 

1.019 

2.975 

16247.0 
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617.6 

3.010 

1.005 

3,000 

16191. 7 

619.7 

3.028 

9.  917xl0-6 

3.018 

16136.8 

621.8 

3.051 

9.784 

3.041 

16082.3 

623.9 

3.077 

9.652 

3.067 

16028.2 

626.1 

3.111 

9.514 

3.101 

15971.9 

628.2 

3.139 

9.381 

3.130 

15918.5 

630.3 

3.163 

9.250 

3.154 

15865.5 

632.4 

3.200 

9.121 

3.191 

15812.8 

634.5 

3.224 

8.995 

3.215 

15760.4 

636.7 

3.263 

8.865 

3.254 

15706.0 

638.9 

3.283 

8.737 

3.274 

15651.9 

640.3 

3,311 

8.657 

3.302 

15617.7 

T 


e.  Experimental  Results 


The  HP9820  calculator  was  programmed  to  compute  attenuation  coefficients 

a(A)  =  a(A) ^  +  a(X)^s  +  a(A)^s  ^  use  Eq,  0-)  •  Five  independent 

measurements  of  A(A)/B(A)  were  made  at  each  of  several  wavelength  positions 

from  which  five  values  for  a (A)  were  obtained.  Values  of  a (A)  presented 

in  the  second  column  of  Table  XVIII  are  averages  of  the  five  independent 

measurements.  The  experimental  error,  expressed  as  percentages  of  a(A), 

were  +  13%  in  the  wavelength  region  from  418,6  to  509.5  nn,  and  +  8% 

in  the  wavelength  region  from  511.8  to  640.3  nm.  Scattering  coefficients 

a(A)Qs  due  to  thermal  density  fluctuations  were  computed  by  use  of  Eq.  (17). 

Values  thus  obtained  for  a (A)  appear  in  the  third  column  of  Table  XVIII. 

ms 

Values  for  the  coefficient  a(A)  +  a(A)  appear  in  the  fourth  column  of 

cl  P  S 

Table  yvill.  Wavelength  A  and  wave  number  A  ^  appear  in  the  first  and  fifth 
columns  of  Table  XVIII,  respectively. 

The  attenuation  coefficients  a(A)  measured  during  the  present  investiga¬ 
tion  are  designated  by  solid  circles  in  Fig.  §>,  where  they  are  graphically 
compared  with  values  of  a(A)  and  a (A)  obtained  by  previous  investigators. 

The  open  circles  in  Fig.  2  represent  a  (A)  of  distilled  water  measured  by 

Sullivan.—^  The  open  triangles  represent  a (A)  for  sea  water  which  were 

9/ 

obtained  from  absorption  spectra  measured  by  Clark  and  James.—  The 
solid  squares  designate  a (A)  computed  for  clear  lake  and  clear  ocean 
water  by  Tyler  et  al.— ^  The  solid  triangles  designate  et(A)  =  4irk(A)/A 
computed  from  k(A)  chosen  by  Hale  and  Querry— ^  for  distilled  water  at  25°C. 


I 

I 

I 


192 


And,  the  open  squares  designate  data  obtained  for  deionized  filtered  water 
by  Hass  and  Davisson^-^  who  used  adiabatic  laser  calorimetry  to  measure 
only  a  (A) 

u 

In  the  spectral  region  from  418.6  to  580  nm,a(A)  from  the  present 
investigation  are  intermediate  to  those  previously  computed  for  clear 
natural  water— ^  and  distilled  water. This  is  due  to  scattering  by 
microscopic  particles  w’hich  are  generally  more  abundant  in  natural  water 
than  in  deionized  filtered  water,  and  are  more  abundant  in  deionized 
filtered  water  than  in  distilled  water.  In  the  spectral  region  from  532 
to  610  nm  ci(A)  from  the  present  investigation  are  slightly  less  than  ct (A) 
obtained  by  Sullivan  for  distilled  water,  and  are  less  than  a  (A) 
computed  by  Tyler  et  al.  for  clear  natural  water.  In  the  spectral  region 
from  610  to  640.3  nm  a (A)  from  the  present  investigation  agree  very  well 
with  ct(A)  obtained  by  Sullivan  for  distilled  water,  are  about  11  percent  greater 
than  a  (A)  from  Hale  and  Querrv,  and  are  less  than  a (A)  computed  for 
clear  natural  water. 

Therefore,  the  split-pulse  laser  method  described  in  this  paper 
provides  a  convenient  means  for  measuring  the  attenuation  coefficients  of 
highly  transparent  liquids. 


... 
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Optical  Properties  of  Crystalline  Graphite  (E]C). 
a.  Introduction 

A  search  of  the  scientific  literature  on  the  optical  pro¬ 
perties  of  carbon  in  any  spectral  region  provided  experimental  data 
from  which  the  imaginary  part  of  N(£)  =  n(£)  +  ik(£)  could  be 
determined  in  the  10- ^ev  <  c  <  lO^ev  spectral  region.  By  use  of 
high  and  low  energy  approximations  the  k(£)  spectrum  thus  obtained 
was  defined  piecewise  analytically  over  the  entire  electromagnetic 
spectrum.  Compton  scattering  and  pair  production  were  ignored  in  the 
high  energy  approximations.  Kramers-Kronig  analysis  of  the  k(£) 
spectrum  provided  values  for  n(£),  which  along  with  knowledge  of 
k(£)  were  used  to  compute  the  reflectance  spectrum,  the  relative 
dielectric  function,  and  the  energy  loss  function.  The  evaluation 
of  five  different  sum  rules  showed  the  data  to  be  of  high  quality. 

Graphite  is  an  optically  uniaxial  material.  In  the  infrared 
spectral  region,  therefore,  two  k(£)  spectra  are  required  to 
completely  specify  the  optical  constants.  The  search  of  the 
scientific  literature,  however,  provided  a  k(£)  spectrum  only  for 
the  electric  vector  of  the  electromagnetic  wave  linearly  polarised 
perpendicular  to  the  c-axis  of  the  graphite  crystal.  The  results 
here  are  only  for  the  case  EjC.  In  the  energy  region  above  £  =  1  ev 
this  constraint  of  ElC  is  not  pertinent. 

These  investigations  of  graphite  are  based  on  a  M.S,  Thesis 
(UMKC,  1979)  by  Mr.  Roy  Holt. 


LOC  ENERGY  <EV) 

A  graphical  presentation  of  log  [k( C ) ]  vs.  log(s) 
In  the  10--3  -  10  ev  spectral  region.  Values 
for  k(c)  were  compiled  from  the  scientific 
literature. 
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b.  Spectra  for  Graphite 

The  k(£)  spectrum  for  graphite  obtained  from  about  30 
papers  in  the  scientific  literature  is  presented  in  Fig,  33. > 
where  the  base  10  log  rithn  of  k(£)  vs,  £  is  plotted.  For 
purposes  of  comparison  with  other  data  the  convesion  x(ym)  = 

(1.24  m  ev/£(ev))  may  be  helpful  to  the  reader.  There  is  seme 

-1  -2 

structure  in  the  far  infrared  region  (10  -10  ev),  valence 
and  conduction  electronic  bands  near  5ev  and  15ev,  and  the 
x-ray  k-absorption  edge  in  the  10  -10  J  region.  The  k(£) 
spectrjm  in  the  10  -40ev  range  is  presented  in  Fig.  84, 
where  the  electronic  bands  are  very  vivid  at  ca.  5  ev  and  10  ev. 
One  sum  rule  evaluated  to  test  the  quality  of  the  k(£) 
spectrum  was 


Z(£q)  - 


m 

h^ne2 


/  £  k(c)  d£  , 

0,001  ev 


where  Z(£Q)  is  the  effective  number  of  electrcns/atom  contributing 
to  the  k(£)  spectrum  in  the  region  0,001  ev  to  £  ,  m  is  the 
electron  rest  mass,  h  is  Plancks  constant  divided  by  2r ,  n  is 
the  number  of  atoms  per  unit  volume,  and  e  is  the  charge  of  the 
electron.  In  the  limit  as  £0  -►  «*,  for  graphite  Z(£Q-+«)  =  6 
the  number  of  electrons/carbon  atom.  A  spectrum  depicting 
the  evaluation  of  the  Z (£ „ )  sum  rule  for  the  k(£ )  spectrum  of 
graphite  is  shown  in  Fig.  85,  where  Z(£Q)  vs.  log(£0)  is  plotted. 
Mote  that  in  the  region  between  10  ev  and  lGJev  the  Z(£  )  reaches 
a  value  of  6  and  then  remains  constant  for  about  2,5  orders  of 
magnitude  in  £  .  This  sum  rule,  together  with  the  evaluation  of 
four  other  different  sum  rules,  gave  us  considerable  confidence 
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in  the  k(£)  spectrum, 

A  Kramers-Kronig  analysis  of  the  k(c)  spectrum  provided 
nCc)  for  graphite.  A  plot  of  logfnCc)]  vs.  log(c)  in  the 
energy  range  from  10-:1  ev  to  10^  ev  is  presented  in  Fig,  86. 

And,  n(£)  vs.  £  in  the  energy  range  from  10~^  ev  to  ^0  ev  is 
presented  in  Fig,  87,  The  normal  incidence  reflectance  spectrum 
R(5)  was  computed  by  use  of  the  Fresnel  equations  and  the 
knowledge  of  k(£)  and  n(£).  The  spectrum  RC C )  vs.  £  in  the 
energy  range  from  10  ev  to  40  ev  is  presented  in  Fig,  88. 

The  reader  is  referred  to  Mr.  Holt's  M.S,  Thesis  for  a 
more  detailed  discussion  of  the  physics  associated  with  this 
invest igaticn. 


ar.d  is  the  contribution  from  the  K-edged  absorption  in  carbon 


ENERGY  (EV) 


Figure  87.  A  graphical  presentation  of  n(c)  vs.  £  in  the 
10  -^0  ev  spectral  region.  n(£)  is  the  real 
part  of  the  complex  refractive  index 
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3.  Thin  Wedge-Shaped  Cell  for  Highly  Absorbent  Liquids* 

I.  L.  Tyler,  G.  Taylor,  and  Marvin  R.  Querry 
University  of  Missouri-Kansas  City, 

Department  of  Physics,  Kansas  City,  Mo.  64110 

a.  Abstract 

We  describe  the  construction,  physical  function,  and  use 
of  a  wedge-shaped  cell  for  measuring  the  Lambert  absorption 
coefficient  of  highly  absorbent  liquids.  Liquids  are  held  in 
the  cell  by  surface  tension  thus  avoiding  the  use  of  seals. 
The  cell  is  simple  in  its  design  and  is  about  as  convenient  to 
use  as  an  ordinary  cuvette. 
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b.  Introduction 

The  fractional  spectral  transmittance  T(v,z)  for  electro¬ 
magnetic  radiation  of  frequency  v  passing  through  an  absorption 
cell  containing  a  volume  of  liquid  of  thickness  z  is 

T(v,z)  =  t(v)^  t(v)^  exp[-2da(v)w-za(\j)zJ  ,  (1) 

where  t(v)^  is  the  spectral  transmittance  at  the  air-window 
interfaces,  t(v)^  is  the  spectral  transmittance  at  the  window- 
liquid  interfaces,  a('J)w  and  aCv)^  are  the  Lambert  absorption 
coefficients  of  the  window  material  w  and  the  liquid  1,  and  d 
is  the  thickness  of  each  of  two  similar  cell  windows.  Two 
measurements  of  T(v,z)  made  with  electromagnetic  radiation  of 
the  same  frequency  v  and  with  the  same  cell  windows,  but  with 
different  thicknesses  z^  and  z^  °£  the  absorbing  liquid,  provide 

a(v)4  =  ln[TCv,z)/TCv,z2)]/C22-^1)  (2) 


In  spectral  regions  where  liquids  are  highly  absorbent  measure¬ 
ments  of  a(v)^  require  differences  in  thicknesses  in  the  range 
0< (z2-z1)f 200ym. 

Robertson  and  Williams— ^  previously  described  the  use  of 
a  unique  absorption  cell  wherein  liquid  water  was  contained  in 
a  small  wedge  shaped  volume  between  two  optically  flat  infrared 
transmitting  windows  of  CaF2  or  KRS-5.  Neoprene  strips  were 
used  to  seal  the  edges  of  the  cell.  The  wedge  varied  uniformly 
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from  zero  thickness  at  its  vertex  to  20-ym  thickness  at  its 
base,  which  was  determined  by  use  of  interferometric  techniques. 

The  cell  was  then  oriented  with  the  window  surfaces  perpendicular 

to  the  optic  axis  of  a  spectrophotometer  system.  As  the  cell 

was  translated  in  a  direction  perpendicular  to  the  optic  axis  the  wedge- 

shaped  volume  presented  uniformly  increasing  thicknesses  z  of 

water  to  a  convergent  beam  of  infrared  radiant  flux.  In  this 

manner  Robertson  and  Williams  used  Eq.  (2)  to  determine  ct(y)^  for 

water  throughout  a  broad  region  of  the  infrared  spectrum. 

The  cell  described  by  Robertson  and  Williams  was  ,  at 

inception,  novel  in  its  application  and  it  served  their  purposes 

well.  There  are  at  least  three  things,  however,  that  limit  its 

versatility  for  general  use  with  liquids  other  than  water.  First, 

the  cell  was  assembled  and  then  sealed  as  a  fixed  unit  with  a 

relatively  small  vertex  angle,  thus  requiring  a  Beckman 

2/ 

variable-pathlength  cell—  for  use  with  liquids  possessing  a 
Lambert  absorption  coefficient  less  than  about  50  cm-'*'.  Second, 
the  presence  of  the  neoprene  seals  would  make  it  difficult  to 
flush  the  cell  when  changing  samples,  and  the  seals  would  also 
be  a  potential  source  of  contaminants.  The  seals,  however,  would 
prevent  evaporation  of  the  liquid  sample.  And  third,  during 
use  the  cell  was  oriented  so  that  possible  interference  fringes 
would  be  parallel  to  the  entrance  slit  of  the  spectrophotometer. 
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This  orientation  required  a  careful  evaluation  of  the  effect 
of  interference  fringes  during  measurements  on  liquids  for  which 
the  Lambert  absorption  coefficient  was  relatively  low  and  for 
which  the  index  of  refraction  was  significantly  different  from 
that  of  the  cell  windows.  The  cell  we  have  designed,  constructed, 
and  used  avoids  the  foregoing  shortcomings  by  use  of  means 
described  in  the  following  sections. 


C.  Description  of  the  Cell 

An  exploded  view  of  our  wedge  shaped  absorption  cell  and 

cell  holder  is  shown  in  Fig.  The  dimensions  provided  in  the 

following  description  are  dictated  by  the  dimensions  of  the  cell 

windows;  thus  minor  alterations  may  be  required  if  one  chooses 

windows  with  different  dimensions.  The  cell  comprises: 

a  square  stainless  steel  base  1  measuring  7.6  x  7.6  x  0.95  cm 

with  a  rectangular  hole  measuring  5.4  x  2.7  x  0.95  cm  centered 

therethrough  as  shown  in  Fig.  89; 

two  2.2  x  2.2  cm  riser  mounts  2'  and  2*'  attached  to  base  1  by 

use  of  machine  screws  as  shown  in  Fig.  39; 

3  / 

toggle  clamps—  3’  and  3M  mounted  on  2'  and  2'',  respectively: 

|  4/ 

spring  plungers—  4*  and  4''  mounted  in  the  heads  of  3'  and  3,r, 

respectively; 

two  rectangular  stainless  steel  blocks  5'  and  5''  measuring 

2.5  x  1.3  x  0.95  cm  with  two  threaded  holes  through  their  sides 
for  receiving  set  screws  6’  and  7’,  and  6"  and  7",  respectively, 
and  with  5'  and  5 ' *  mounted  on  base  1  as  shown  in  Fig.  89; 

1 

\ 

___  __  a*.  r  ■  — 

J 

43*  ^  ^  y,  *'  *■**►*->■  -  "*•»<»' ' — ■ -«.■ —  •  — 


'The  windows  determine  the  cell  dimensions 
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two  rectangular  stainless  steel  blocks  8'  and  8"  measuring 

2.5  x  1.3  x  0.95  cm,  with  two  threaded  holes  in  each  for 

3/ 

receiving  threaded  ball  plungers—  9'  and  10',  and  9M  and 
10'r,  respectively,  with  a  slot  0.8  cm  deep  and  0.16  cm 
wide  in  each  block,  and  with  8’  and  8”  mounted  on  base  1 
as  shown  in  Fig.-  89; 

two  stainless  steel  posts  11’  and  11"  each  0.32  cm  in  diameter 
and  1.6  cm  long  mounted  2.5  cm  and  2.2  cm  from  edges  of 
base  1  as  shown  in  Fig. 

two  posts  12'  and  12"  each  0.32  cm  in  diameter  and  0.95  cm 
long,  and  mounted  0.6  cm  and  1.9  cm  from  edges  of  base  1 
near  blocks  8'  and  8",  respectively; 
two  square  optically  flat  windows^  13'  and  13"  measuring 
5.09  x  5.08  x  0.64  cm; 

a  segment,  of  machine  shop  feeler  stock— ^  14  measuring  about 

9.5  x  1.3  cm  and  of  a  thickness  suitable  to  the  needs  of  the 
investigator;  and 

a  handle  15  centered  and  mounted  on  a  side  surface  of  base  1 
which  provides  a  means  for  holding  the  cell  for  translation. 
The  cell  is  assembled  empty  as  follows.  The  bottom  surface 
of  window  13'  rests  flat  against  the  top  surface  of  base  1  so 
that  13'  is  squarely  centered  on  and  covering  the  rectangular 
hole  in  base  1.  The  rounded  ends  of  set  screws  7’  and  7" 
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contact  one  edge  of  window  13',  The  ball  ends  of  ball  plungers 
10'  and  10 ’ '  contact  the  opposite  edge  of  13'  thus  acting  as  a 
resilient  means  for  holding  13'  in  place.  Window  13'  easily 
snaps  in  place  when  a  slight  force  is  applied  near  the  edge  that 
contacts  10*  and  10''.  The  slots  in  8'  and  8"  receive  the 
segment  of  feeler  stock  14  as  indicated  in  Fig  89,  A  portion 
of  the  lower  surface  of  14  rests  flat  against  the  top  surface 
of  window  13'  as  indicated  in  Fig.  89.  Posts  12'  and  12"  retain 
14  in  slots  8'  and  8".  The  second  window  13"  is  then  gently 
pressed  in  position  so' that  one  bottom  edge  of  13 "  contacts  the 
top  surface  of  window  13'  thus  forming  the  vertex  of  a  wedge- 
shaped  volume  between  13*  and  13".  The  opposite  edge  of  window 
13"  rests  flat  against  the  upper  surface  of  feeler  stock  14 
thus  forming  the  base  of  the  wedge-shaped  volume  between  13 '  and 
13".  Pvounded  ends  of  set  screws  6'  and  6"  contact  one  edge  of 
window  13".  The  ball  ends  of  ball  plungers  9'  and  9"  contact 
the  opposite  edge  of  13"  thus  providing  a  resilient  means  for  hold¬ 
ing  13"  in  place.  Set  screws  6'  and  6"  are  advanced  slightly 
more  than  screws  7'  and  7"  so  that  the  edge  of  window  13"  rests 
on  the  flat  upper  surface  of  window  13'.  One  side  of  windows  13' 
and  13"  are  each  in  contact  with  posts  11*  and  11".  Toggle 
clamps  3'  and  3"  are  then  closed  so  that  spring  plungers  4* 
and  4"  gently  but  resiliently  contact  points  on  the  top  surface 
of  window  13";  this  secures  contact  between  1  and  13',  13'  and 
14,  13"  and  14,  and  13'  anc  13". 


I 

I 
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After  assembly,  the  vertex  angle  of  the  wedge-shaped  volume 
between  windows  13’  and  13  *  *  was  measured  by  use  of  a  He-Ne  laser 
interferometer.  We  determined  by  repeatedly  assembling  and 
disassembling  the  cell  that  the  vertex  angle  of  the  wedge  shaped 
volume  was  reproducible  to  better  thair-.+  0.3  per  cent  if 
windows  13'  and  13 ' ’  were  not  rotated  or  interchanged. 

A  liquid  is  placed  in  the  cell  during  assembly  in  the  following 
manner.  A  few  droplets  of  liquid  are  deposited  on  the  upper 
surface  of  window  13’  near  the  edge  that  forms  the  vertex  of 
the  wedge-shaped  volume  between  13'  and  13".  Window  13" 
is  then  gently  pressed  in  position  thus  deforming  the  liquid 
so  it  fills  the  wedge-shaped  volume  between  13’  and  13".  Toggle 
clamps  3*  and  3"  are  then  closed.  Most  feeler  stock  is  made 
of  ordinary  steel  so  care  must  be  exercised  to  avoid  overfilling 
the  cell  when  corrosive  liquids  are  used , 

During  us^  we  place  the  cell  just  in  front  of  the  entrance 
slit  of  a  Perkin  -Elmer  El  Ebert  monochromator  so  that  the 
vertex  is  uppermost  and  so  that  window  13’  is  perpendicular  to 
a  beam  of  convergent  radiant  flux  that  passes  through  the  rectangular 
hole  in  base  1.  The  limiting  aperture  of  the  optical  system  is 
formed  horizontally  by  the  entrance  slit  of  the  monochromator, 
which  is  continuously  adjustable  from  10  to  10^  pm,  and  is  formed 
vertically  by  a  0.64  cm  wide  slit  milled  in  a  plate  placed  between 
the  cell  and  the  entrance  slit.  Thus,  possible  interference  fringes 


I 

I 


from  the  wedge-shaped  volume  between  windows  13'  and  13’ ' 


are 
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perpendicular  to  the  entrance  slit.  The  cell  is  translated 
vertically  by  use  of  an  xy-translating  stage.  Interference 
effects  are  thereby  avoided  when  the  entrance  slit  is  narrowed 
for  investigations  requiring  greater  spectral  resolution. 

d.  Calculation  of  Maximum  Cell  Thickness 

Fig. 90  is  two  cross-sectional  drawings  wherein  the  windows, 

feeler  stock  of  thickness  D,  and  liquid  sample  are  shown  in  the 

vertical  orientation  in  which  the  cell  is  used  for  measurements 

of  a(v)^.  The  liquid  is  retained  in  the  cell  by  surface  tension 

Ywi  of  the  liquid  in  contact  with  the  cell  windows,  y^  of  the 

liquid  in  contact  with  its  vapor  and  y  of  the  vapor  in  contact 

wv 

with  the  cell  windows.  Fig.  90a  shows  the  liquid  completely 
filling  the  vertex  region  as  it  does  when  retained  in  the  cell. 

Fig 90b  shows  a  virtual  downward  displacement  x  of  the  top  surface 
of  the  liquid  from  the  vertex.  The  corresponding  virtual  change 
in  surface  potential  energy  Ac ,  per  unit  length  of  liquid  in  a 
direction  perpendicular  to  the  plane  of  Fig.90,  and  to  first  order 
in  x ,  is 

“  ■  1Dt1v/L  +  .  <3> 

where  x  and  L  are  shown  in  Fig.  2.  Ac  is  positive  for  liquids 
that  wet  the  windows  because  (y^-y...)  £  y ^  for  such  liquids. 

r 


l 


(a)  (b) 

Figure  90.  Cross-sectional  view  of  the  cell  windows,  liquid 
.-rsampie,  and  feeler-stcck  spacer  shown  in  the 
vertical  orientation  in  which  the  cell  is  used. 

In  (a)  the  sample  fills  the  vertex  region  between 
the. windows.  In  (b)  the  sample  is  virtually  dis¬ 
placed  a  distance  x  from  the  vertex. 


i 
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Furthermore,  in  the  right-hand  side  of  Eq.  (3)  the  first  term 
is  orders  of  magnitude  less  than  the  second  term.  The  upward 
force  on  the  liquid  due  to  surface  tension  is  provided  by  the 
negative  derivative  of  Ac  with  respect  to  x.  The  weight  W  of 
the  liquid,  per  unit  length  is 

W  =  DS2pg/2L  ,  (4) 

where  S  is  shown  in  Fig.  2A,  p  is  the  volume  density,  and  g  is 
the  acceleration  due  to  gravity.  The  maximum  value  of  D  occurs 
when  the  upward  force  due  to  surface  tension  is  equal  in  magnitude 
to  the  downward  force  W.  Thus  we  find,  letting  S  equal  to  L 

DMAX  3  4^V/Lpg  .  (5) 

Maximum  values  of  D  calculated  by  use  of  Eq.  (5)  for  some  common 
liquids  retained  in  a  cell  with  glass  windows  for  which  L  =  5  cm 
are  presented  in  Table  XIX.  However,  our  experience  is  that  liquids 
are  retained  in  the  cell  even  when  values  of  D  are  about  30% 
larger  than  those  calculated  by  use  of  Eq.  (5).  This  is  probably 
due  to  surface  tension  effects  at  the  vertex  where  the  two  windows 
are  in  contac^as  well  as  to  our  ignoring  the  first  term  in  Eq.  (3). 

e.  Some  Experimental  Results 

Table  XX  presents  values  of  aCv)^  that  we  measured  for  20, 

30,  and  40  per  cent  aqueous  solutions  of  zinc  chloride  and  for 
50,  65,  75,  and  85  per  cent  phosphoric  acid  at 
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TABLE  XIX 

Calculated  upper  bounds  for  the  cell  thickness  D  for  some 
common  solvents.  Values  were  obtained  using  Eq .  (5)  with 
L  =  5  cm. 

Solvent  D[un] 

Acetone  240 
Benzene  270 
Carbon  Tetrachloride  140 
Ethyl  Alcohol  230 


Water 


590 


223 


TABLE  XX 

Lambert  absorption  coefficients  a(v).  and  values  for  k(v)  = 
a(v)^/4rv  at  v  =  4444.4  cm-'*'  for  aqueous  solutions  of  ZnCl2 
and  for  phosphoric  acid. 


Sample 

a(V)  l 

k(v)£ 

20%  ZnCl2 

22.06  +  0.08  cm-1 

(3.950  +  0.014)xl0-4 

30%  ZnCl? 

21.73  +  0.07 

(3.891  +  0.012) 

40%  ZnCl2 

21.64  +  0.08 

(3.874  +  0.015) 

5C%  H3P04 

75.29  +  1.56 

(13.48  +  0.28) 

65%  E3?04 

93.99  +  0.39 

(16.83  +  0.07) 

75%  H3P04 

108.40  +  0.22 

(19.41  +  0.04) 

85%  H3P04 

117.28  +  0.34 

(21.00  +  0.06) 
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v  =  4444.4  cm  \  Values  of  k(v)^  =  4Trva(v)^>  the  imaginary  part 

A 

of  the  complex  refractive  index  N(v)  =  n(v)+ik(v) ,  are  also 

presented  in  Table  XX.  We  use  values  of  k(v)£  at  v  =  4444,4  cm"1 

in  a  double  subtractive  Kramers-Kronig  analysis  of  reflectance 

spectra  to  determine  N(v)  for  the  zinc  chloride  solutions— ^  and 

*  8  / 

for  the  phosphoric  acids—  throughout  the  spectral  region 
270  <_  v  <_  5000  cm'1. 

f.  Conclusion 

We  constructed  and  tested  a  wedge  shape  absorption  cell 
which  is  extremely  versatile.  By  use  of  the  cell  one  can  easily 
and  accurately  measure  the  Lambert  absorption  coefficient  of 
highly  absorbent  liquids.  The  provision  of  ’a  variable  apex  angle 
permits  cell  parameters  to  be  optimized  for  any  given  sample. 

The  cell  can  be  used  to  measure  absorption  in  any  spectral 
region  because  several  different  sets  of  windows  are  freely 
interchangeable.  And,  the  cell  is  relatively  inexpensive  to 
construct  and  maintain. 
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4.  Microprocessor  controlled  Interface  fcr  a  Desktop  Calculator 
a.  Abstract 

We  have  designed  and  constructed  a  microprocessor-based  interface/controller 
which  significantly  expands  the  capabilities  of  a  programmable  desktop 
calculator.  With  the  aid  of  the  unit,  a  calculator  such  as  the  Hewlett- 
Packard  Model  9820  equipped  with  an  xy-piotter  and  cassette  tape  deck 
(a)  controls  an  optical  spectrometer,  (b)  reads  and  stores  optical  data 
from  the  spectrometer,  (c)  transmits  data  to  a  central  computer  via  a  ISO 
terminal,  and  (d)  receives  data  from  the  central  computer  and  plots  them 
on  the  xy  plotter.  The  techniques  employed  are  quite  general  and  the 
system,  as  built,  can  be  used  to  automate  a  wide  variety  of  laboratory 
equipment.  The  capability  of  communicating  with  the  central  computer  also 
makes  the  system  a  convenient  substitute  for  a  "CALCOMP"  plotter. 
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b .  Introduction 

i'I)  Control  of  Experimental  Equipment. 

In  the  past  year  or  so  there  has  been  a  marked  increase  in  the 
uumber  of  commercial  instruments  which  perform  a  fixed  number  of 
specified  complicated  experimental  operations.  The  heart  of  the  control 
system  for  these  units  is  a  microcomputer  (MC).  Usually  the  only 
feature  which  makes  these  instruments  more  than  their  predecessors  is 
the  automation  provided  by  the  MC  and  associated  electronics.  Therefore 
automation  of  a  laboratory  procedure  often  involves  replacing  perfectly 
good  laboratory  equipment  with  expensive  new  equipment  for  the  sake 
of  a  small  amount  of  control  electronics.  In  addition  one  may  later 
find  that  the  fixed  number  of  control  algorithms  provided  by  the  new 
instrument  no  longer  are  sufficient. 

To  avoid  these  problems  some  investigators  have  begun  to  construct 
MC-based  control  units  for  their  existing  equipment.— ^  The  controllers 
are  inexpensive  and  can  be  quite  versatile  but  they  may  require  more 
familiarity  with  microprocessor  programming  procedures  than  is  reasonable 
for  these  who  use  the  equipment  only  occasionally.  They  also  are  not 
well-adapted  for  performing  intricate  decision-making  procedures 
sometimes  required  during  experiments. 

The  ideal  control  unit  would  be  one  which  was  easy  for  anyone  to 
operate  and  yet  retained  the  versatility  of  a  user-programmable  MC. 

The  newest  generation  of  desktop,  programmable  calculators  approaches 
this  ideal  but  only  at  a  very  high  cost  ($12,000-315,000). 

| 

I 
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Table  XXI 


List  of  abbreviations  used  in  the  text 


AC 

cc 

Cl 

CPU 

ECZ 

FP1 

FP2 

HP  9820 

LD 

MC 

PAD 

ROM 

RWM 

5DS 

5M 

3MD 

3T 

SR 

CART 


CPU  accumulator  register 
Central  computer 
Communications  Interface 
MC  central  processing  unit 
Experimental  Control  Interface 
Front  Panel  of  Cl 
Front  Panel  of  ECI 

Hewlett-Packard  Model  9820  calculator 

Load-in-memory  strobe  pulse  from  CPU 

Microcomputer 

Peripheral  Address  Decoder 

CC  read-only  memory 

CC  read-write  memory 

Serial  Data  Switcher 

Stepper  motor 

Serial  Mode  Decoder 

Store-in-memory  strobe  pulse  from  CPU 
ECI  shift  registers 
Universal  asynchronous  receiver- 
transmitter 


We  have  chosen  an  alternative,  less  expensive,  approach.  An 
MC  is  utilized  as  an  interface  between  experimental  equipment  and 
one  of  the  previous  generation  of  desktop  calculators  —  a  Hewlett- 
Packard  Model  9820  (HP  9820).  In  this  approach  the  MC  is  programmed 
to  interpret  and  to  carry  out  a  specified  set  of  commands  entered 
on  the  calculator.  The  entered  commands  are  in  simplified  mnemonic 
form  thus  permitting  anyone  familiar  with  the  calculator  to  automate 
experimental  procedures  once  the  MC  has  been  programmed.  Furthermore, 
the  calculator  can  be  programmed  in  the  usual  way  to  determine,  for 
example,  the  sequence  of  scan  increment  sizes  to  scan  a  spectrometer 
linearly  in  wavelength  or  in  wavenumber.  Because  the  desktop  calculator 
retains  its  normal  operating  capabilities,  it  can  also  carry  out 
complicated  mathematical  calculations  on  experimental  data  during  an 
experimental  run.  These  calculations  can  then  be  used  to  modify  the 
experimental  procedure  automatically  during  the  run.  For  example,  it 
is  possible  for  an  operator,  familiar  only  with  the  standard  calculator 
programming  procedures,  to  program  a  routine  which  automatically 
locates  absorption  peaks,  rescans  through  these  peaks  with  enhanced 
resolution,  and  finally  calculates  relative  peak  heights,  widths,  and 
areas. 

(2)  Communication  with  a  Central  Computer 

An  additional  advantage  of  the  interface  we  have  constructed  is 
that  it  can  be  used  to  connect  the  programmable  calculator  to  a 
central  computer  (CC)  via  a  modem  and  telephone  circuit.  This  not 
only  permits  tne  rapid  transmission  of  experimental  data  stored  on 
cassette  tapes  but  also  permits  the  CC  to  send  the  results  of  calcula¬ 
tions  back  ro  the  calculator.  "here  are,  of  course,  manv  instances 


when  a  medium-resolution  plot  of  computer  calculations  is  much 
preferable  to  the  standard  computer  printout.  The  delays  and  other 
aggravations  involved  in  obtaining  a  "CALCOMP"  plot  of  CC  calculations 
are  well  known.  Our  interface  permits  us  to  plot  CC  output  directly 
on  the  plotter  accessory  of  the  calculator.  The  time  required  to  obtain 
a  plot  of  400  datum  points  with  a  resolution  of  better  than  0.252  is 
thirty  seconds  —  a  significant  improvement  over  the  usual  "CALCOMP" 


turnaround  time. 
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c.  system  Design 

An  interconnection  diagram  of  the  various  units  comprising  the  complete 

system  is  shown  in  Fig.  91.  The  additions  described  in  this  article  are  those 

enclosed  within  the  dotted  lines.  The  Communication  Interface  (Cl)  contains 

2/ 

the  microprocessor  chip—  ,  1/2  K  of  read-only  memory,  and  3  K  of  read-write 
memory  as  well  as  the  decoders,  associated  electronics,  and  switches 
referred  to  in  the  next  subsection.  The  Experiment  Control  Interface  (ECI) 
is  under  direct  control  of  the  Cl.  It  facilitates  operation  of  the  spectro¬ 
meter  and  data  acquisition  system  but  can  be  omitted  at  the  expense  of 
additional  software.  The  system  was  expressly  constructed  for  use  with 
the  HP  9820  but  could  be  modified  for  use  with  any  ASCII-coded  programmable 
calculator  equipped  with  parallel  data  input/output  capability. 

In  the  following  subsections  we  describe  first  the  principles  utilized 
in  communicating  between  the  various  units  in  the  system  and  secondly  the 
function  of  the  two  interfaces. 

(1)  Principles  of  Operation 

For  the  purposes  of  this  discussion  the  MC  can  be  represented  as  a 
box  with  three  sets  of  lines  emanating  from  it  —  sixteen  address  lines, 
two  control  lines,  and  eight  data  lines.  As  indicated  in  Fig.  32, 
the  MC  itself  consists  of  three  sections.  The  first  is  the  central  pro¬ 
cessing  unit  (CPU)  which  contains  an  accumulator  (AC)  and  several  general 
purpose  registers  each  capable  of  holding  one  8-bit  binary  data  byte. 

The  second  is  a  random  access  or  read-write  memory  (RWM)  which 
contains  a  large  number  of  8-bit  memory  cells  each  labeled  with  a 
16-bit  binary  address.  The  third  is  a  read-only  memory  (ROM)  which 


Figure  92.  A  block  diagram  Indicating  the  electronic  principles 
employed  to  control  peripherial  devices. 


contains  permanent  data  used  to  program  the  microcomputer.  The  CPU 
is  connected  internally  to  the  ROM  and  RWM  by  the  address  and  data 
lines  discussed  above. 

During  the  execution  of  a  program  the  CPU  transfers  data  from 

the  RWM  into  its  registers,  performs  logical  and  arithmetic  operations 

on  these  data,  and  restores  the  modified  data  to  the  original  or  some 

other  address  in  the  RWM.  The  procedure  for  programming  the  MC  to 

3/ 

perform  desired  operations  is  not  important  here.—  What  is 
important  is  to  note  how  the  CPU  communicates  with  the  RWM  and  how 
this  communication  procedure  can  be  utilized  to  enable  the  CPU  to 
communicate  with  peripheral  devices. 

In  sending  data  from  the  AC  to  one  of  the  addressed  memory  cells 
in  the  RWM  the  CPU  first  enters  the  address  of  that  cell  on  the  16  address 
lines  which  internally  connect  the  CPU  and  RWM.  Then  it  puts  the 
data  in  the  AC  onto  the  data  lines.  Following  this  the  CPU  strobes 
the  ST  (store)  control  line  low  for  one  microsecond.  This  strobing 
of  the  control  line  effects  the  transfer  of  the  data  on  the  8-data 
lines  into  the  specified  RWM  cell.  Loading  data  from  a  specified  WRM 
register  into  the  AC  register  is  accomplished  as  follows:  The  CPU  sends 
an  address  via  the  address  lines  to  the  RWM  and  then  strobes  the  load 
control  line  (LD)  low  for  one  microsecond.  When  I.D  is  low  the  WRM 
puts  the  data  in  the  addressed  cell  onto  the  data  lines.  When  LD  goes 
high  the  CPU  loads  these  data  into  the  AC. 

Because  the  address,  strobe,  and  data  lines  are  externally  accessible 
it  is  possible  to  substitute  external  de\rices  —  peripherals  —  for 
some  of  the  memory  ceils.  Usually  the  substitution  is  for  cells  with 
addresses  lOOOOOOOGOOOCOOO  to  1111111111111111.  Each  peripheral  can 
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nave  as  many  as  eight  input/output  leads  for  each  address  register  it 
replaces.  This  procedure  is  called  memory-mapped  I/O  (Input/output). 

Peripherals  are  interfaced  to  the  microcomputer  by  the  method 
shown  in  Fig.  92.  First  input  and  output  ports  (buffers)  controlled 
oy  the  LD  and  ST  pulses  from  the  CPU  are  interposed  between  the 
microcomputer  parallel  data  lines  and  the  peripheral  device.  These 
prevent  the  peripherals  from  interfering  with  data  appearing  on  these 
lines  during  normal  operation  of  the  MC.  Secondly,  a  peripheral 
address  decoder  (PAD)  consisting  of  standard  decoder  chips  is  connected 
to  the  strobe  and  address  lines  coming  from  the  MC.  The  PAD  routes 
the  LD  and  ST  pulses  to  an  appropriate  pair  of  outputs  corresponding 
to  the  memory  address  placed  by  the  CPU  on  the  address  lines  at  the 
time  the  pulse  appears.  The  control  pulse  is  then  directed  to  an 
input  port  (ST)  or  output  port  (LD)  of  the  peripheral  with  the  same 
address, thereby  activating  the  port  and  effecting  transfer  of  data 
between  the  peripheral  and  the  MC  data  lines. 

As  an  example  of  how  this  interfacing  can  be  used  we  consider 
the  operation  of  a  stepper  motor  (SM).  A  typical  SM  has  four  windings, 
a,  b,  c,  and  d,  each  of  which  is  energized  and  de-energized  in  a  pre¬ 
scribed  sequence  to  advance  the  rotor.  Letting  a  logic  1  (+5v)  represent 
an  energized  winding  and  logic  0  (OV)  represent  a  de-energized 
winding,  clockwise  shaft  rotation  can  be  achieved  by  cyclic  repetition 
of  the  four  patterns  0101,  0110,  1010,  and  1001.  In  these  four  binary 
numbers  the  least  significant  bit  represents  winding  a  and  so  forth. 

The  SM  windings  can  then  be  wired,  through  power  amplifiers,  to  four 

of  the  eight  output  data  lines  from  an  addressed  output  port  and 

latch.  Operation  of  the  SM  then  simply  requires  programming  the  MC 

to  send  the  data  prescribed  above,  at  a  rate  proportional  to  the 
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desired  SM  speed,  to  the  address  of  the  output  port. 

In  practice  most  peripheral  devices  require  more  than  one  address 
for  operation.  For  example  data  communication  between  the  MC  and  the 
CC  is  effected  through  a  universal  asynchronous  receiver-transmitter 
(UART)  chip.  The  chip  receives  data  from  CC  in  series  form  at  a  rate 
of  3C0  bits  sec  Only  after  a  full  8-bits  of  data  have  been  received 
can  the  UART  transmit  these  data  in  parallel  via  the  data  lines  to  the 
MC.  It  is  therefore  necessary  for  the  MC  to  wait  until  it  receives 
an  indication  from  the  UART  that  the  data  are  ready  before  it  executes 
a  LD  command  to  the  UART  data  address.  For  this  reason  the  UART  is 
equipped  with  status  lines.  These  lines  are  connected  to  an  output 
port  with  another  address.  During  the  execution  of  a  subroutine  to 
receive  data  from  the  CC  the  MC  is  programmed  to  load  and  analyze  the 
data  from  the  UART  status  address  repeatedly  until  these  data  indicate 
that  real  data  are  available  at  the  UART  data  output  port. 

Similarly  the  HP  9820  parallel  interface  communicates  with  the 
MC  via  two  addresses,  the  second  serving  for  status  communication. 

A  simple  way  to  control  a  wide  variety  of  devices  is  shown 
connected  to  outputs  (1)  and  (2)  of  the  PAD  in  Fig.  91.  This  is  an 
on-off  switch  consisting  of  a  set-reset  flip-flop  which  is  driven 
to  logic  1  (ON)  by  a  ST  command  from  the  MC  to  its  "ON"  address  and 
to  a  logic  0  (OFF)  by  a  ST  command  to  its  "OFF"  address.  The  advantage 
of  this  device  is  that  it  is  a  simple  matter  to  include  a  mechanical 
switch  circuit  in  a  manner  such  that  the  last  action  taken  —  mechanical- 
or  MC  derived  --  determines  the  state  of  the  output.  Such  an  arrange¬ 
ment  is  very  convenient  in  the  operation  of  electromagnetic  clutches, 
etc.  If  AC  power  must  be  controlled  the  set-reset  flip-flop  can  be 
used  to  drive  a  solid  state  relay. 


(2)  Circuit  Description 


In  this  section  we  describe  the  circuitry  using  only  block  diagrams. 
Detailed  circuit  diagrams  are  available  from  the  authors  on  request. 

(a.)  Communications  Interface 

A  block  diagram  of  the  MC  Cl  is  shown  in  Fig.  93.  Input  lines 
terminating  in  bubbles  represent  control  strobe  lines  which  originate  in 
the  PAD  (upper  left)  and  are  used  to  activate  input  and  output  ports. 

Broad  open  lines  represent  parallel  data  lines. 

All  such  lines  are  buffered  by  input  or 

output  ports  before,  connecting  with  the  MC  bidirectional  parallel  data 
lines  whether  shown  explicitely  or  not.  Broad  dark  lines  represent 
buffered  device  status  lines. 

The  HP  9820  is  connected  to  the  system  through  both  parallel  and 
serial  interfaces.  The  serial  interface  is  used  only  to  transmit  data 
from  the  calculator  to  the  central  computer.  As  discussed  earlier,  data 
from  the  CC  is  transmitted  serially  at  300  bits  sec  \  The  HP  9820 
cannct  receive  serial  data  faster  than  260  bits  sec  ^  hence  the  UART  chip 
is  used  as  an  intermediary.  Serial  data  are  received  by  the  UART.  Once 
a  full  8-bit  byte  has  been  received  it  is  transmitted  in  parallel  form 
to  the  HP  9820. 

The  Serial  Data  Switcher  (SDS)  is  worth  considering  in  greater  detail. 
Commercial  computer  devices  often  interface  with  RS-232  logic  (1  =  -10v, 

0  =  -*-10v) .  For  Completely  automated  operation  it  necessary  that  the  MC 
be  capable  of  mediating  communication  not  only  between  TTL  peripherals 
(1  =*  +5v,  0  ■  0V)  but  also  between  TTL  and  RS-232  peripherals.  This  is 
accomplished  by  converting  all  RS-232  input  signals  to  TTL  by  means  of 
RS-232  line  receiver  chips.  All  internal  switching  is  then  effected  by 
the  microcomputer  in  TTL  via  the  Serial  Mode  Decoder  (SMT).  The  SMD  converts 
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Figure  93.  A  block  diagram  o 
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the  four  least  significant  bits  on  the  data  line  into  their  BCD  equiva¬ 
lent  when  the  SMD  address  is  strobed.  The  TTL  portion  of  the  SDS  con¬ 
sists  of  a  number  of  logic  gates.  These  gates  are  controlled  by  the 
SMD  BCD  output.  Each  BCD  digit  thus  selects  one  of  16  possible  gate 
configurations  each  determining  a  particular  routing  of  serial  data 
from  peripherals  through  the  SDS  and  out  to  other  peripherals.  Finally 
those  data  routed  out  to  RS-232  peripherals  are  converted  back  to  RS-232 
by  means  of  RS-232  line  drivers  before  leaving  the  SDS. 

This  approach  provides  a  general  method  for  doing  any  degree  of 
complex  MC-controlled  switching  using  low-cost  TTL  components.  Switching 
between  RS-232  devices  is  normally  cumbersome  and  expensive  requiring 
reed  relays  or  manual  switches.  A  well-shielded  junction  box,  however, 
is  not  subject  to  noise  levels  requiring  +  10V  logic.  Thus  by  con¬ 
verting  to  TTL  for  switching  one  can  reduce  costs  by  an  order-of- 
magnitude  without  sacrificing  reliability. 

b)  Experimental  Control  Interface 

A  block  diagram  of  the  Experiment  Control  Interface  (ECI)  is 
shown  in  Fig.  The  same  conventions  are  used  in  this  diagram  as  were 
used  in  Fig.  93.  Thus,  each  block  in  the  diagram  to  which  bubbles  are 
attached  represents  a  memory  mapped  I/O  device,  and  so  forth. 

The  ECI  is  unaer  direct  control  of  the  MC.  The  MC,  in  turn,  is 
under  the  control  of  the  HP  9820  calculator.  The  MC  is  programmed  to 
check  repetitively  for  ASCII-coded  parallel  input  (commands)  from  the 
HP  9820.  Once  the  command  is  received  the  MC  analyzes  the  command 
characters  and  then  jumps  to  that  particular  subroutine  stored  in  the 
RWM  required  to  execute  the  command.  Execution  involves  sending,  in 
the  proper  sequence,  one  or  more  ST  commands  to  memory  addresses  of 
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peripherals  in  the  ECI.  Once  the  subroutine  is  executed  the  MC  again 
begins  a  repetitive  check  for  new  commands  from  the  HP  9820,  This 
provides  interlocked  operation  between  the  MC  program  and  the  HP  9820 
program. 

The  ECI  is  currently  being  used  to  control  a  Perkin-Elmer  El  mono¬ 
chromator  and  a  data  acquisition  system.  The  monochromator  grating  is 
scanned  by  an  SM  drive.  The  processed  signal  from  the  monochromotor 
detector  is  digitized  by  a  PAH  Model  260  A/D  converter  under  control 
of  memory  mapped  I/O  signals. 

As  was  seated  in  the  introduction,  the  versatility  of  the  MC-HP  9320 
combination  allows  anyone  familiar  with  the  HP  9820  programming  procedures 
to  automate  extremely  complicated  scan  routines.  This  is  accomplished 
by  utilizing  the  HP  9820  to  perform  complex  calculations  and,  on  the 
basis  of  these  calculations,  to  determine  a  desired  wavelength  and  wave¬ 
number  region  to  be  scanned  as  well  as  the  number  of  scans  and  spacing 
of  data  points.  The  MC  then  need  only  be  programmed  to  execute  a 
very  small  numoer  of  commands  such  as  "scan  up  n  steps"  or  "read  the 
A/D  converter".  Each  such  command  is  written  in  mnemonic  form  and  the 
operator  is  supplied  with  this  short  list  of  "calculator  keyboard 
operations"  in  addition  to  the  calculator's  usual  keyboard  functions. 

In  what  follows  we  describe  the  operation  of  the  ECI  by  listing 
the  HP  9820  mnemonic  commands,  following  each  with  a  short  command 
description  and  finally  listing  the  sequential  operations  carried  out 
by  the  system  to  effect  the  command. 

9/ 

SL’xxxx:—  Scan  up  xxxx  steps  where  xxxx  ranges  from  0000  to  9999. 
a)  The  number  xxxx  is  series- loaded  from  the  MC  into  the  ECI  shift 
registers  (SR). 
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b)  This. number  is  parallel- loaded  into  the  ECI  latches  for  front 
panel  display  and  into  the  ECI  A  Step  Counter. 

c)  The  MC  activates  the  SM  scan  via  the  ECI  Scan  Control  Logic. 

The  pulses  running  the  SM  are  also  directed  to  the  A  Step  Counter 
and  to  the  Step  Counter.  Each  pulse  decrements  the  number  in  the 
A  Step  Counter  by  one  digit.  When  this  number  reaches  zero  a 
pulse  is  generated  within  the  A  Step  Counter.  This  pulse  is 
directed  to  the  Scan  Control  Logic  where  it  serves  to  stop  the 
scan.  The  MC  monitors  the  Scan  Control  Logic  unit  to  determine 

'  when  the  Scan  has  stopped,—^ 

d)  Each  pulse  received  by  the  Step  Counter  increments  the  number 
stored  there  by  one  digit.  At  the  end  of  the  scan  the  MC  parallel 
loads  this  number,  via  the  2;1  multiplexer  and  SR  into  the  Latches 
for  front  panel  display. 

SDxxxx:  Scan  down  xxxx  steps. 

The  sequential  operations  are  the  same  as  SUxxxx  except  for 

obvious  changes, 

RA:  Read  the  A/D  converter, 

a)  The  MC  commands  the  A/D  converter  to  digitize  the  signal  voltage 
appearing  at  its  input. 

b)  These  digitized  data  are  parallel- loaded  via  the  2:1  multi¬ 
plexer  into  the  SR. 

c)  From  the  SR  the  data  are  loaded  serially  into  the  MC. 

d)  The  MC  waits  for  an  input  (read)  command  from  the  HP  9820  after 
which  it  sends  the  data  in  parallel  to  the  HP  9820  for  analysis 


or  storage. 
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RS:  Read  Step  Counter  and  transmit  the  number  to  the  HP  9820. 

a)  The  number  in  the  Step  Counter  is  parallel-loaded  into  the  SR 

via  the  2:1  Multiplexer. 

b)  From  the  SR  the  number  is  serial-loaded  into  the  MC. 

c)  The  MC  then  waits  for  an  input  (reac$  command  from  the  HP  9820 

after  which  it  sends  the  number  in  parallel  to  the  HP  9820. 

P  xxxx:  Pause  for  xxx.x  seconds  where  xxx.x  ranges  from  000.1-999.9. 

The  MC  is  equipped  by  the  manufacturer  with  a  pause  command  in  the 
ROM.  The  subroutine  for  thi3  command  simply  translates  the  P  xxx.x  into 
the  equivalent  assembler  command. 

The  following  commands  are  effected  by  a  single  ST  command  to  the 
appropriate  peripheral  addresses 

CS  Clear  Step  Counter  and  Step  Counter  Display. 

CD  Clear  A  Step  Counter  and  A  Step  Counter  Display. 

A  Turn  on  External  Alarm  to  alert  operator. 

CA  Turn  off  External  Alarm. 

Cc)  Operational  Problems 

We  conclude  this  section  with  a  discussion  of  some  unexpected 
problems  which  were  encountered  in  the  operation  of  this  interface  system 
and  of  how  these  problems  were  overcome. 

(1)  Echo  during  direct  derial  communication  between  the  HP  9820  and 
the  CC. 

The  HP  9820  Serial  Interface  cannot  be  used  without  the  SDS  to 
receive  signals  directly  from  the  CC.  This  is  because  all  received 
signals  are  echoed  back  out  to  the  CC  via  the  HP  serial  output  line 


and  interfere  with  the  CC  operation.  Even  though  the  HP  9320 
is  not  used  to  receive  data  directly  from  the  CC  it  _is  used  to 
receive  "handshake"  signals.  By  routing  HP  9820  serial  input/output 
through  the  SDS  this  interference  can  be  avoided.  The  SDS  logic 
circuit  simply  blocks  the  serial  output  channel  when  it  is  placed 
in  the  HP  9820  serial  input  mode  by  the  MC. 

1)  Interference  of  "handshake"  operations  between  the  HP  9820  and  CC. 

In  the  transmission  of  data  from  the  HP  9820  to  the  CC  it  is 
necessary  both  for  the  CC  to  be  notified  when  one  line  of  data  has 
been  entered  cn  the  TSO  terminal  and  for  HP  9820  to  be  notified 
when  the  entered  data  have  been  received.  We  utilize  the  standard 
"handshake"  signals  between  TSO  and  CC  for  this  purpose.  Thus  when 
the  HP  9820  has  entered  a  full  line  of  data  on  the  Decvriter  terminal 
it  terminates  the  data  with  the  ASCII  equivalent  of  CNTL/D  entered 
on  the  keyboard.  As  usual,  when  the  CC  has  received  these  data 
it  unlocks  the  TSO  keyboard,  returns  the  carriage,  and  advances 
the  paper  with  the  ASCII  equivalent  of  CR/LF  (carriage  return/ line 
feed)  entered  on  the  keyboard.  The  HP  9820  is  programmed  to  look 
for  the  ASCII  characters  CR/LF.  When  these  are  received  it 
begins  entering  another  line  of  data  on  the  keyboard,  and  so  forth. 

A  problem  arises  because  when  the  HP  9820  sends  serial  data 
it  first  enters  them  in  parallel  on  internal  shift  registers. 

It  then  proceeds  to  execute  the  next  programmed  command  while  these 
data  are  being  shifted  out  serially.  Normally  the  next  programmed 
command  would  be  to  receive  the  CR/LF  signal  from  the  CC;  but  this 
involves  having  the  SDS  block  the  HP  serial  output  channel  before 
the  serial  output  data  have  left  the  shift  registers.  This  problem 
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can  be  avoided  by  programming  the  HP  9820  to  execute  several 

I 

1  "mark-time"  operations  between  the  send  data  and  receive  data 

commands. 
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(d)  Performance 

The  microprocessor-based  interface/controller  has  been  in  operation 
for  over  two  years.  It  has  proven  to  be  extremely  versatile  and 
virtually  trouble-free.  The  construction  process  itself  was  not  trivial 
but  an  undergraduate  with  a  solid  background  in  digital  electronics  should 
be  able  to  build  a  similar  unit  as  a  summer  project.  The  cost  of  the 
components  required  would  be  about  $800, 
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We  used  a  SC/MP  microprocessor  manufactured  by  National  Semiconductor 
Corp,  Santa  Clara,  CA.  Complete  microprocessor-plus-memory 
boards  are  now  available  from  a  number  of  manufacturers. 

See  for  example,  W.  J.  Weller,  Assembly  Lend  Programming  for  Small 
Computer,  (Lexington  Books,  Lexington,  MA. ,  1975).  H 

In  order  to  transmit  the  ASCII  characters  SUxxxx  via  the  EP  9820 
parallel  interface  the  actual  keyboard  operations  are  FMT  "SUxxxx"; 
WRT  8  where  8  is  the  Select-Code  for  the  HP  Parallel  Interface. 

The  procedure  described  here  is  primarily  "hardware  controlled”  and » 
differs  from  the  "software  controlled"  stepper  motor  scan  procedure 
described  earlier  in  this  article.  However,  we  found  during  the 
course  of  construction  of  this  unit  that  software-dominated  control 
is  usually  preferable  in  terms  of  total  time  required  for  con¬ 


struction  and  programming 


